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Abstract 

It is the purpose of this research to examine the engineering process required to 
adequately predict the total resistance of unique hullforms. A comparison of the significant 
physical effects present in monohull and advanced hullforms testing is made. The historically 
accepted method of scaling using Froude’s Hypothesis is studied and its underlying 
assumptions are considered in detail. An overall design philosophy is created for achieving 
resistance minimization through a range of operating speeds. This discussion is applied to a 
specific hullform called SLICE which is an adaptation of the Small Waterplane Area Twin 
Hull (SWATH) concept. 

To further develop the body of available SLICE resistance data, three models were 
built and tested at the US Naval Academy's 380 feet towing tank in .Annapolis Maryland. The 
results of these tests is reported on and a video demonstration of their performance is 
provided as part of this report. 

Finally, a numerical method approach to the solution of Mitchell’s integral is presented 
as a possible means of calculating SLICE wave making resistance. This is considered as an 
intermediate step tool in design optimization, useful after preliminary model tests are 
completed and before final designs are decided for powered model testing. 
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Chapter 1 Introduction 



The majority of the world's commerce is moved by monohull displacement 
ships. Ironically, in this age of supersonic air travel, guaranteed overnight mail 
delivery' and fibre optic data transfer at the speed of light, modem displacement ships 
are capable of speeds no faster than those built fifty years ago. The apparent lack of 
progress in high speed ocean travel stems from the physical realities of the interaction 
between the displacement monohull form and the water in which it travels. 

The total resistance curve for a ship is generally plotted as the Total Drag 
Coefficient. C, versus Froude Number. JC These are expressed as: 



^v/.rTgLj 



(i-D 



fl r .(ipK 2 S)C, (1-2) 

where L is some characteristic length, usually length of the waterline of the body, p is 
the density of the water. V is the ship’s velocity and S is the Wetted Surface Area. 
Total resistance is generally recognized as the additive effects of friction, wave 
making and spray generation of both the main body and its appendages. The extent to 
which each of these components effects the total resistance is. however, the subject of 
ongoing debate even though it has been over one hundred years since Froude gave us 
this theory. For a monohull displacement ship, the wave making resistance portion of 
the total resistance curve has the characteristic shape shown in Figure (1-1). Since 
resistance is obtained by multiplying the coefficient shown in Figure (1-1) by V 2 , the 
total resistance curve becomes nearly vertical and practically insurmountable by 
conventional hull forms at the speeds shown. 
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Figure 1-1: Characteristic Wave Making Resistance Curve 
(from LMSC presentation! 



For Naval Architects, the search is ongoing in advanced huLl form design for a 
practical method to "get over the hump". Most advanced hull forms do not resemble 
the classic displacement monohull shape on which the vast majority of Naval 
Architecture experience is based. The purpose of this research is to examine the 
underlying assumptions regarding commonly accepted full scale resistance prediction 
methods in an effort to understand the important features of design optimization and 
draw some conclusions regarding resistance minimization methods . The Lockheed 
Missile and Space Company in conjunction with Pacific Marine (LMSGPM) has 
proposed a hull form which is a modification of Small Waterplane Area Twin Hull 
(SWATH) concepts. Their proposal, called SLICE, will be presented to the United 
States Government as a 177 Lton vessel in an Advanced Technology Demonstration 
(ATD) in June 1996. A sketch of this hullform is shown in Figure 1-2. The SLICE 
vessel will be the subject hullform for the research conducted as part of this paper. 
SLICE is protected by Lockheed under patent law. 






I 

I 



Figure 1-2: Artist's Concept of SLICE from an unpublished LMSC, 

Pacific Marine. ONR Article. 




1.1 The SLICE Concept 

High speed is not an inherent quality of SWATH ships. Although there may 
be some reduction of the wave making resistance due to smaller waterplane area, there 
is also an increase in the total wetted surface which increases frictional resistance. In 
SWATH ships, high speed is a design feature. When compared to an equivalent 
displacement monohull. SWATH design must be considered as a trade-off between a 
definite increase in frictional resistance and the potential for a decrease in wave 
making resistance. The design question becomes how to ensure that the decrease in 
wave making resistance is greater than the gain in frictional resistance. 

LM SC's SLICE concept differs from traditional SWATH design in that the 
demi-hull is replaced by two shorter pods. Each pod is attached to the main box 
structure by a single strut. It is through this design feature that LMSC intends to 
minimize the impact of wave making resistance. By designing the vessel with an 
unusually short characteristic length, the hump of the wave making curve will in 




theory be driven to a ver\ low speed regime. Since resistance is a function of V 2 , any 
left shift of the curve could have a substantial minimizing effect on the total 
resistance associated with the hump. It is further theorized that in the high speed 
regime, total resistance will be dominated by frictional effects and therefore design 
efforts focused on minimizing the wetted surface of the hullform will yield the best 
result. LMSC has chosen the length of the submerged pod as the characteristic length 
critical to design considerations. It should be emphasized that this hullform does not 
eliminate the hump effect characteristic of the wave making resistance curve. It does 
not specifically cause wave cancellation effects nor does it attempt to reduce viscous 
drag effects. Rather, the apparent thnist of the LMSC theory is to stagger the peak 
effects of wave and skin resistance, thereby reducing the overall resistance value at 
specific design speeds. 

1.2 LMSC Model Test Synopsis 

LMSC has conducted two different sets of model tests at the US Navy’s 
Surface Warfare Center, Carderock Division. Bethesda, Maryland (NSWC). In the 
first senes of tests, called the HM&E tests, several different configurations of struts, 
pods and foils were used to make a determination of the overall best configuration to 
use as the baseline for the production vessel. These tests were conducted in the Spring 
1994. 

The second set of tests, also conducted at NSWC. were used to obtain detailed 
resistance and control surface performance data for the baseline production vessel. 
These tests, conducted in August 1994 and continued in December 1994 were 
observed by several MIT students. The model that was tested was a 1 /8th scale of the 
production vessel. It’s construction was fiberglass bodies mounted to an aluminum 
box frame. For all the tests conducted, the model was fixed in pitch and heave and 
data was collected for forces and moments about the three principle axes. There were 
two major groupings of data collection runs performed. The first group was called the 
dynamic control tests and was primarily used to judge the required effectiveness of 
control surfaces to hold the ship at a constant draft with an even trim. Each run had a 
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slightly different model configuration with canards, stabilizers and rudders tested in 
different locations and with varying angles of attack. Each configuration was tested a 
single time unless there was some specific problem noted in the conduct of the test. 
The second major grouping of tests consisted of rather major modifications to the 
model configuration in an attempt to improve the design. Haunches on top of the struts 
were removed to allow for deeper strut length, wedges located at the top of the struts 
were added and removed to experiment with better surface piercing shapes, the overall 
draft of the model was adjusted to force a deeper submergence of the rather bluff pods 
and finally, a new model was tested in which the longitudinal position of the struts 
relative to the pods was altered to modify the flow over the pods. 

Within these two major test categories, the test procedure remained consistent. 
The model remained fixed in pitch and heave for all tests. Two different speed 
regimes were tested, a lower speed run consisted of twelve speed increments ranging 
from 3.2- 8.8 feet per second (fps) corresponding to full scale speeds of 9-25 knots. A 
high speed run used seven increments from 5.7-12.4 fps which corresponds to full 
scale speeds of 10-35 knots. All speed increments were fit into a single pass along a 
tow tank length of approximately 1200 feet. All of the dynamic control tests were 
conducted only in the lower speed regime, the major model modification tests were 
conducted in both speed regimes. 

1.3 Research Motivation 

The US Navy has potential applications for relatively small (200-700 Lton) 
high endurance craft that embody the seakeeping attributes of SWATH. The Navy’s 
interest in the SLICE program is demonstrated by the funding the Navy has committed 
to the ATD that LMSC will perform in June 1996. It was LMSC/PM's entrepreneurial 
spirit that initially brought the SLICE concept to the attention of the Navy’ but this 
same entrepreneurial spirit has the potential to limit the usefulness of the SLICE 
design as well. LMSC/PM's allocated timetable from concept development to final 
product delivery’ is just over two years, very’ short by Navy shipbuilding standards. As 
a result LMSC must focus on a relatively narrow set of engineering problems that 



specifically address the needs of the 177 Lton ATD vessel. Broader based questions, 
such at an in-depth understanding of the effects of body interactions, how these 
interactions vary with modifications to the basic dimensions, the practicality of scaling 
the SLICE concept into a different size regime (1500 Lton. 5000 Lton. etc.) and a 
relative comparison of the merits of a SLICE ship as compared to a SWATH or 
Monohull design for a given mission are all facets of the development of this concept 
that are not specifically addressed by LMSC in the ATD. 

This thesis will attempt to add to the developing body of SLICE knowledge by 
focusing on the resistance and powering aspects of the SLICE vessel in an 
academically rigorous maimer. In so doing, three benefits will be realized: 

1) Independent research into this hull form will give both the US Navy and 
LMSC a second hydrodynamic viewpoint of the 177 Lton vessel prior to the ATD. 

2) Development of a sufficient body of data will aid optimization of the 177 
Lton hull form. 

3) A broader understanding of the hydrodynamic features of the hull form 
will facilitate deciding the merits of pursuing designs in other displacements more 
practical to the US Navy. 

1.4 Approach 

To achieve the objectives outlined above, there are several questions raised by 
the LMSC efforts that require answers. These questions and the method that will be 
pursued for answering them is outlined in this section. 

THEORY 

LMSC theorizes that the appropriate characteristic length for SLICE is the 
length of a single pod. In conventional designs. is based on the length of the 
waterline (LWL). Representative characteristic lengths for SLICE would seem to be 
bounded on the low end by the length of a single strut and on the high end by the 
length between the leading edge of the forward strut to the trailing edge of the aft 
strut. LMSC's choice of an underwater body for the characteristic length is not an 



obvious one. Understanding the appropriate characteristic length for the SLICE 
configuration is necessary to understand the magnitude of savings that can be expected 
over competing hull forms. 

The interaction between bodies is the second area for further study. There are 
two specific situations that require detailed analysis. The first is the effect a strut and 
its attached pod have on one another, the second is the effect that a strut/pod 
combination have on any of the other three strut/pod combinations that form the 
underwater body. Understanding these effects is necessary for design optimization and 
also deciding the practical upper limit of displacement for a SLICE ship. 

Finally, the LMSC proposed design places the forward struts inboard of the aft 
struts as illustrated in Figure 2. This arrangement may cause some significant 
outward lateral forces on the struts that must be considered in the structural design. 

TESTING 

There were several procedural questions raised by the model tests conducted by 
LMSC in August. Keeping the model fixed in pitch and heave was the most 
significant issue. Although it was intentionally done to facilitate the dynamic control 
analysis, it also left open significant questions regarding how the hull form behaves at 
different speeds. Considering the recognized sensitivity that SWATH type hulls have 
in tons per inch immersion (TP1) and the moment to trim one inch (MT1) this may be 
significant. 

The short duration of each test run speed increment, the lack of repetition of 
tests, and the truncation of most data collection runs at the full scale equivalent of 25 
knots for a ship required to achieve 30 knots raises questions regarding the 
thoroughness of the results. Although this is understandable considering the expense of 
model testing and the limited time available, more testing must be completed before 
the SLICE body of data can be utilized to make design decisions. 

OPTIMIZING 

There are several striking features incorporated in the LMSC design. The most 
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striking is the bluffness of the leading sections of the pods and struts, especially for a 
high speed craft. It would appear that minimizing the wetted surface to displacement 
ratio was LjMSC's driving concern in hullform selection. Although this approach 
minimizes frictional resistance, it also introduces significant form drag. With a pod 
length to diameter ratio (L/'D) of 4.2, SLICE cannot be thought of as a slender body. 
This introduces another concern in evaluating SLICE performance. The vast majority 
of traditional resistance analysis has slender body theory incorporated to some extent. 
Any evaluation process undertaken must be carefully scrutinized to ensure that the 
SLICE hullform does not violate this imbedded assumption. A design approach which 
more thoroughly addresses the entire speed spectrum is preferable even with a specific 
mission speed identified. All speeds from zero through the mission speed must be 
attainable so some design consideration must be given to the entire speed spectrum. 
With eight component bodies ( 4 struts and 4 pods), wave cancellation effects through 
body interaction should be a design focus. 

RESEARCH METHOD 

In an attempt to answer these questions the following work was accomplished: 

1) A literature review that encompassed one hundred years of ship resistance 
research was performed. This was accomplished to ensure that the full body of 
assumptions made in this field of engineering is understood. 

2) A study of model testing techniques was made to develop an appreciation 
for proper test methods. 

3) Senes 58 resistance data was studied for potential pod optimization. 

4) A literature survey was conducted to ascertain the state of the science in 
spray drag predictions. 

5) Three models were built and tested to expand the overall data bank. 

6) A numerical calculation method for wave making resistance based on the 
solution to Mitchell’s integral was attempted in an effort to expedite the hullform 
optimization process. 
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Chapter 2 Theory 



The difficult) in predicting full scale resistance is in understanding the extent 
to which different physical phenomena contribute to total resistance. Although effects 
such as boundary layer development, wave making, spray generation and flow 
separation are known to be present as a body moves through the water, methods for 
quantitatively evaluating the contribution of each of these effects to total resistance 
and the extent to which these phenomenon may interact continues to be the subject of 
ongoing debate. 

This chapter studies the two available approaches to resistance prediction, 
model testing and analysis The objective is to identify those methods that will be best 
suited for SLICE resistance predictions. Since the vast majority of available resistance 
data was developed from displacement monohull analysis, care must be taken to 
ensure that underlying assumptions to these methods are not violated. 

2.1 Model Testing 

The classical approach to model testing is well documented. Principles of 
Naval Architecture 1 (PNA), and Rawson and Tupper are probably the two most 
recognized authorities outlining this process. The textbook approach is definitively 
stated and easy to follow Nonetheless, the SLICE experience has shown that specific 
model testing methodology and subsequent data reduction can be hotly contested. The 
difficulty lies in deciding if the model tested accurately replicates the physical 
phenomena that will be present in the full scale ship and what to do with the model 
scale data once it is obtained. 

As background, the model testing/scaling process outlined by William Froude 
is presented. In this discussion the subscript m refers to the model and the subscript s 



1 Comstock. J P . Principles of Naval Architecture . SNAME, New York, 1967. Chapter 7 
Rawson, K J & Tupper. E C Basic Ship Theory . 3rd Edition. Vol 2, New York. 1983 
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refers to the full scale ship 

1) Build a model that is geometrically similar to the full scale ship. This is 
accomplished by choosing a scale factor to apply to all linear dimensions of the ship. 
The scale factor X is then defined as: 




^ 'model 



( 2 - 1 ) 



2) Test the model at equivalent speeds. Through dimensional analysis, it is 
easily shown that the full scale and model scale speeds are related through equation 
( 1-1 ). Specifically: 



^ship V model X V^ 



( 2 - 2 ) 



3) Convert the total model resistance measured R lm to a coefficient of total 
resistance C, m through relationship (1-2). 

4) Calculate the coefficient of frictional resistance C tm . The method for 
accomplishing this will be discussed shortly. 

5) A quantity called the coefficient of residuary resistance is obtained by 
subtracting C tm from C lnr The reason C tn , is not subjected to the scaling law expressed 
in (2-2) will be discussed shortly. 

6) The full scale coefficient of residuary resistance C rs is then expressed at the 
same Froude Number as C^,. Specifically. 

C =C (2-3) 

rs rm 

7) The full scale coefficient of frictional resistance C ls is calculated in the 
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same manner as in step 4 and is added to C rs to obtain the full scale coefficient of 
total resistance C ls . 

8) Full scale total resistance R IS is then found using (1-2). 

This process can be thought of as a system of one equation and one unknown. 
The unknown is C^; the equation is C u = + C ts . This method of scaling data is 

useful only if the following conditions exist: 

1) A geometrically similar model can be built and tested in compliance 
with the equalities set forth in (2-1) and (2-2). 

2) The model is tested in a condition that can be considered equivalent 
to the normal operating condition of the full scale ship. 

3) Frictional resistance can be reasonably calculated at both the model 
and full scale size. 

4) All other physical phenomenon present must either obey equality 
(2-3) and therefore be lumped into the residuary component or be of such little 
concern that ignoring them is inconsequential to the final outcome. 

If the last condition can not be satisfied for a particular phenomenon, then its 
effect must be calculated at both the model and full scale size and algebraically treated 
the same as the friction resistance component. 

Froude's experiments were conducted on slender bodies, generally thought of 
as those with L/D > 10. Residuary resistance was defined as the combined effect of 
wave making resistance and eddy resistance. At low speeds, eddies appear to be a by 
product of the boundary layer and therefore part of frictional resistance. 
Experimentally, Froude demonstrated that eddy resistance did not change 
proportional ly with skin friction as parallel midbodv was added to a given hull form. 
He concluded that eddy resistance was actually related to the form of the hull. Since 
wave making and eddy resistance are both related to form and are pressure related 
phenomenon, they historically have remained combined as the residuary resistance 



component' For slender bodies the only significant physicaJ phenomena are friction, 
waves and eddy formation. The process outlined above was developed with slender 
bodies and only these three phenomena in mind. 

2.2 Frictional Resistance 

The problem of calculating the frictional resistance of a ship was first studied 
by considering flat plates. The essential issue focuses on the development of the 
boundary layer along the length of the plate and its transition from a laminar to a 
turbulent regime. Frictional resistance is usually expressed as a coefficient of friction 
in the form of (1-2) and is generally plotted as a function of the Reynolds Number, '.H,: 



8 = 



V*L 



v 



(2-4) 



where L is some characteristic length, usually length of the waterline of the body, o 
is the kinematic viscosity of the water and V is the ship's velocity. The development 
of a laminar flow equation was the work of Blasius (1904). while the turbulent 
solution was formulated independently by Prandtl and von Karman in the 1920's. In 
the range of Reynolds Number between 5x10' and 9xl0 6 , neither the turbulent nor the 
laminar expression yields accurate results. This region of uncertainty is known as the 
transition region. Experimentally it has been shown that the transition from laminar to 
turbulent flow occurs gradually across the length of a plate as its speed increases. This 
gives physical meaning to the rather large band of uncertainty between fully laminar 
and fully turbulent flow. The curves of laminar, transition and turbulent plate flow are 
shown in Figure (2-1). 

Model testing through the 1940's attempted to draw parallels between monohull 



3 Gilmer.T C and Johnson. B. Introduction to Naval Architecture . Naval Institute Press, Annapolis. 
1982, p 217 
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torms and the known Hat plate solution. The most accepted flat plate solution to 
frictional resistance adapted to ship form use is the Schoenherr Formula. This flat 
plate formulation is frequently seen in the literature labeled as the ATTC Line. In an 
attempt to account tor the effect of hull form on boundary layer formation the 
International Towing Tank Conference (ITTC) adopted an interim solution that w-as 
derived from a correlation of previously obtained model tests. 



;009 

. ITT 




Figure 2-1: Laminar and Turbulent Flow Skin Friction Lines 4 



Principles of N’avai Architecture , op. cit , p295 
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This formulation, known as the ITTC 1957 line is given as: 




0.075 

flog, 0 3i-2) 2 



(2-5) 



It should be stressed that this line was developed from a correlation of model data, it 
does not physically represent the behavior of a flat plate and it was not obtained by 
extrapolating flat plate data to ship like forms. Forty years later, the ITTC has not 
replaced this interim solution with a better alternative. 

At the same time that the ITTC 195" formulation was published. Hughes 
proposed a method for calculating flat plate frictional resistance given bv: 




0.066 

(log I0 3i-2.03) 2 



( 2 - 6 ) 



To the basic Hughes line, a form factor is then added. The purpose of this form factor 
is to help account for the differences in the laminar to turbulent transition seen by a 
flat plate and a fuller body. The attractiveness of Hugh's method is that it attempts to 
keep friction and eddy resistance linked. Physical observation of a body moving 
through the water would seem to make this reasonable. The method proposed by 
Hughes has not received the same acceptance within the Naval Architecture 
community as the ITTC 1957 line. The earliest criticism was the low value that 
equation (2-6) yields when compared to flat plate testing. Additionally, the use of a 
form factor is also questionable as outlined by Oosterveld\ quoted below': 

" It has now become clear that in many cases the wave 
pattern resistance is less than the residuary resistance estimated 



' Oosterveld. M W C. ed . "Report of the Resistance Committee". 15th International Towing Tank 
Conference . The Hague, 1*578. p 21 



by assuming a form factor independent of Froude number and 
Reynolds number. This will certainly be the case whenever bow 
wave breaking occurs. One possible reason for this is that the 
viscous resistance may vary significantly with Froude number 
(which effects the flow at the edge of the boundary layer) as 
well as with Reynolds number, so that it is not well predicted on 
a constant form factor basis. However it is conceivable that this 
variation arises at least in pan from an interaction of the 
boundary layer and wake with the wave system, so that energy 
losses which originate with wave making, and are primarily 
Froude number dependent, may manifest themselves as an 
apparent increase in the viscous wake losses." 

The frictional lines discussed in this section are shown in Figure (2-2). 




5 « lO* 10* 



1 » J J i 1 



Figure 2-2: Standard Skin Friction Lines 



This discussion of frictional resistance has bearing on the SLICE project 



6 Principles of Naval Architecture , op.cn . p. 298 



respects First there is the issue of turbulent flow. The general practice in model 
testing is to artificially stimulate turbulence. This is done so that the model and the 
full scale ship operate in the same regime. Without turbulence tripping, most models 
would operate at least partially in the transition regime shown in Figure (2-1), making 
data interpretation difficult since C, is highly variable in this regime. Hughes' 
describes standard methods of stimulation (usually pins, wires or sand strips), 
methods of accounting for the resistance caused by these devices and the uncertainty 
involved in insuring that turbulence is actually achieved. The values of 9? at which the 
SLICE models and full scale ship operate are shown in Table (2-1). A comparison of 
these values to Figure (2-1) would appear to demand that some form of turbulence 
tripping be used in testing. 
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Table 2-1: Operating Reynolds Numbers for SLICE 
based on Strut and Pod Length 

Figure (2-1) however shows the laminar to turbulent transition for a flat plate. SLICE, 
with a nearly hemispherical nose shape and an overall L/D * 4.5 cannot be likened to 
a flat plate. The laminar to turbulent transition for a sphere is known to occur in the 



Hughes, G and Allan, J , "Turbulence Stimulation on Ship Models. SNAME Transactions, 1951, 
p28 1-3 14 
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region 'Ji = 4.5 x 10 J - 5 x 1CP or two orders of magnitude lower than that of a flat 
plate. The actual transition for SLICE probably lies somewhere between these two 
extreme boundaries. When viewed in this light, turbulence stimulation is probably not 
worth the uncertainty involved in mathematically removing the resistance effects of the 
trip device. For the purposes of testing in this project, at the model scale, turbulent 
flow will be assumed even though turbulence stimulation is not provided 

The second issue is to decide the most appropriate method to compute the 
frictional component of SLICE resistance. The ATTC line should be excluded since it 
is derived strictly from flat plate results. The Hughes line with form factor corrections 
may initially appear attractive since it incorporates calculated adjustments for form. 

The fullness of the SLICE body and the definite appearance of a breaking bow wave 
may cause the inaccurate accounting discussed by Oosterveld. The ITTC 1957 line is 
the best choice. It is based on actual ship and model results. Even though there may 
be uncertainty in describing the details of laminar to turbulent transition, and the 
impact the fullness of a ship shape has on that transition, the ITTC 1957 line must be 
thought of as correctly including the phenomenon since it is based solely on test 
results. This is a fallout of the model to ship correlation method used to derive the 
line. It is the only line discussed that was developed in this fashion Further work in 
this project will utilize the ITTC 1957 line as the basis for calculating frictional 
resistance. 

2.3 Wave Making Resistance 

Wave making resistance has two parts. In an obvious wav it accounts for the 
energy expended creating the visible wave system seen as a body moves through the 
water. It also accounts for the energy expended overcoming drag due to the pressure 
distribution created by unseparated inviscid flow around a body. The literature is 
vague on the second part of this definition of wave making resistance. The pressure 
distribution around the body is directly related to the form of the body but this effect 
is not the commonly termed "form drag". "Form drag" is a Reynolds scaled viscous 
effect that is best thought of as the correlation between flat plate frictional resistance 
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prediction and actual ship friction results. The form part of wave making resistance 
should not be confused with eddy resistance either which accounts for the effects of 
separated flow. Since both parts of the definition of wave making resistance are 
pressure effects, this phenomenon will follow Froude scaling. 

One of the reasons that model testing is performed is because the ability to 
calculate wave making resistance remains an elusive goal. This section reviews the 
classical work that has been accomplished in this area. The objective is to gain some 
insights into the physical phenomenon and how it applies to the SLICE hullform. 

The foundation for a mathematical approach to wave resistance was set in place at 
nearly the same time that Froude was developing his model testing methods. J.FL 
Mitchell gave an analytic solution for wave resistance of the form s 



R I sec 3 0l f f aexp{(-^-)sec 2 d(v-ixcosd)}dxd\’ 2 dd (2-7) 

2 V 2 J JJ V 2 



Where p is the density of water. V is velocity and a is the potential function that 
represents the body of interest. The combined complexity of writing a to comply with 
the necessary boundary conditions at the free surface as well as with the relatively 
complex form of a monohull and the difficulty in actually solving this integral without 
the advantage of numerical computing techniques, were serious setbacks to the 
popularity of this approach. The most significant pursuit of Mitchell's work was 
undertaken by T. Havelock and W.C. Wiglev. Wiglev was the first to successfully 
decompose the waves generated by a hullform into components and thus gain an 
understanding of how they are formed. This was first achieved with what is now 
called the Wiglev hull, a wedge followed by parallel midbody followed by an identical 
wedge pointed opposite to the bow. Eventually Wiglev refined the hull to a smooth 



* Newman. J.N., Marine Hydrodynamics . The MIT Press. Cambridge. 1989. p 282 
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form. He discovered that there are actually four waves produced by this shape. The 
first is produced by the leading edge of the hull, the next two are produced by the two 
shoulders formed by the intersection of the fore and aft wedge with the parallel 
midbodv and the final one is produced by the trailing edge of the aft wedge. This 
decomposition held true for hullforms with continuous lines as well. Figure (2-3) is 
reproduced from Wiglev's writings 9 which graphically demonstrates this additive 
effect. 

This discussion is useful to the SLICE program because it gives insights to an 
approach at hullfonn optimization. Since the amplitude of the wave system generated 
is related to the form of the body, maintaining lines that are reasonably fine should 
help to minimize the hump of the wave making curs e. As a approximation without 
actual calculations, the shape of the struts and pods of the SLICE form will be 
adjusted in an effort to find a reasonable compromise between frictional and wave 
making resistance. This is different from the LMSC approach which appears to have 
concentrated on the operating range of ^ > 0.9 with effort focused on minimizing 
frictional resistance. 

The SLICE hullform can be decomposed into eight different bodies, four struts 
and four pods. Each of these component bodies has a form that is reasonably similar 
to the hullforms tested by Wigley. If each of these bodies can be thought of as 
producing four waves that combine into a single wave system for the body, then an 
approach to optimizing wave making resistance can be created by adjusting body 
spacing once form issues are decided. A three tier system of looking at this problem is 
proposed: 

1) First consider the interaction between a strut and its attached pod. This step 
can be thought of as creating an oversized bulbous bow for the strut. 

2) Next consider the effect of longitudinally aligned strut pod 
combinations. This will probably have a diminished effectiveness compared to step 1 

w Wigley, W C S, "A Comparison of Experiment and Calculated Wave-Profiles and Wave- 
Resistances for a form having Parabolic Waterlines", Proceedings of the Roval Society . Series A, 
London, 1934. pp 144-159. 
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Figure 2-3: Calculated Wave Profile Analysis 



since the amplitude of the forward system will not equal that of the aft system once 
they interact. 

3) Finally consider the interactions of the transversely and diagonally aligned 
strut pod combinations. Since these are not in the principle direction of motion, it is 
expected that the effectiveness of these interactions will again be diminished from the 
previous step. 

2.4 Eddy Resistance 

As previously noted, the literature is noticeably vague in defining both form 
and eddy resistance. They are frequently used together but are not used 
interchangeably Articles written about resistance prediction usually account for eddy 
resistance as an input to residuary resistance. The research for this thesis uncovered no 
direct usage of form resistance in prediction schemes. As discussed in Chapter 2-2, 
this thesis considers form drag to be the increase in frictional resistance felt by a ship 
shape compared to an equivalent flat plate. The 1TTC 1957 line directly accounts for 
form drag since it is based on experimental data from monohull forms. Other schemes 
such as the ATTC line and Hughes Line require some form compensation for accurate 
use since they are based on flat plate resistance. 

Eddy resistance on the other hand has physical origins. In this thesis it is taken 
to mean the result of separation effects from discontinuities in the hull form. These 
can result from actual breaks in the ships lines such as the intersection of the strut and 
pod at both the leading and trailing edge, intersection of appendages and the main 
body, or radii of curvature which are too hard for the flow to follow. As discussed in 
chapter 2-1, eddy resistance is a pressure effect and therefore should be scaled as part 
of the residuary resistance component. 

Following the classic approach, the minor swirling seen along the sides of a 
ship are associated with viscous effects of well behaved flow and are accounted for in 
the frictional component' 0 . Significant swirling that can be associated with flow 
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Principles of Naval Architecture , op cit p 312. 



separation, such as that seen at the transom, describes the type of effect accounted for 
by eddy resistance. There is some obvious room for interpretation of cause and effects 
in these definitions of form and eddy drag. This highlights the importance of observing 
physical behavior through model testing. 

2.5 Spray Resistance 

The formation of spray and the impact it has on total drag has historically been 
assumed insignificant to monohull testing. Simple observation of full and model scale 
monohulls operating at speed validates this assumption Although spray is generated, 
the region of generation is highly localized and the quantity of spray generated is 
usually small. With the increased interest in advanced hullforms such as hydrofoils 
and planing craft, neglecting spray was recognized as a poor assumption but the body 
of available knowledge regarding the mechanism of spray generation and the 
appropriate method of scaling spray results is relatively small. 

When spray becomes a major concern, the one equation, one unknown analogy 
used in section 2-1 becomes one equation and two unknowns. The typical solution to 
this problem has been to test at very high speeds {£F > 3.0) to eliminate wave making. 
This allows spray drag to be analyzed separately but is only a partial answer. A gap 
exists at the lower, more practical values of cF. In 1971, Chapman 11 conducted a 
significant test program on strut forms with parabolic leading and trailing edges. He 
concluded the following: 

1 . Spray drag is partially dependent on both strut form and thickness. Blunt 
leading edges produce more spray. 

2. The spray mass flow rate of spray depends on the thickness to chord ratio 
(t/c) and on Blunt bodies send the spray up at higher angles causing greater 
wetting. 

3. Skin friction due to wetting is the primary source of spray drag. 



11 Chapman. R.B. , Spray Drag of Surface Piercing Struts . Naval Undersea Research and 
Development Center Report AD 730710. 1971 
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4. The momentum in the spray sheet is sufficient to create a large pressure 
drag on any object it strikes. 

Chapman’s work was expanded by Zhu Bing-quan in 1 986 12 . A total of twelve 
struts were tested, all with parabolic trailing and leading edges. The objective was to 
correlate strut performance to the parameters c (chord length), t ( maximum thickness) 
and x ( longitudinal position of the maximum thickness). The relationship derived is: 



R=pV 2 x 2 C 



(2-8) 



and C 5 is given as: 



C =0.068- -0.555(-) 2 +0.696(-) 3 +0.083(-) 2 ^? r (2-9) 

t t t t 

The research also showed that the spray drag was minimized when lx = 0.313. 

This has applications to the SLICE analysis. The strut form selected by LMSC 
is contrary to that recommended by the published literature for minimizing spray drag. 
The apparent design goal of minimizing frictional resistance fails to consider spray 
effects. In deciding the extent to which the spray effects will scale from model size to 
full scale, the dependence of (2-9) on ^implies that spray drag will follow Froude 
scaling. The conclusion reached by Chapman and Zhu Bing-quan that the initiation of 
spray is highly form dependent also supports this assertion. 

The following is recommended when dealing with spray effects. Choose strut 
shapes that are known to minimize spray generation. Unless there is a significant 
structural or interior arrangements reason to do otherwise, struts should have parabolic 
leading and trailing edges and lx should be reasonably close to 0.313. For the ATD, 



** Zhu Bing-quan, Ge Wei-zhen. "Experimental Study on Spray Performance of Surface Piercing 
Struts. Journal of Hydrodynamics , ser 2 no 4. 1990, pp 91-98 
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this will help minimize spray effects and allow the research efforts to focus on the 
area of principle interest, the wave making resistance. Any spray that is observed in 
testing should be treated as part of the residuary resistance. If the tentative conclusion 
that spray follows Froude scaling is wrong, then this process will over predict the total 
resistance. Considering the uniqueness of the SLICE hullform, a conservative error of 
this nature is well justified for the first full scale ship. Design work in this project will 
follow these conclusions. 

2.6 Appendage Resistance 

The effect of appendages on total drag has historically been assumed small in 
monolnill testing. This assumption is based on the fact that the overall size of the 
appendages is usually one or two orders of magnitude less than the main hull. For 
comparative testing of hulls, it is reasonable to ignore the appendages. For obtaining 
exact resistance predictions, appendages are usually included in the final testing and 
stripping tests have been devised to help account for their effects. In stripping tests, 
the model is towed first with all the appendages on. In subsequent tests, the 
appendages are removed one at a time. In this manner the impact of each appendage 
can be assessed. This approach assumes that the interaction effects of the appendages 
with each other and with the main hull are negligible. 

For SLICE, appendages occupy a larger percentage of total volume compared 
to monohulls. The flow around SLICE is more complex and therefore interactions may 
be more pronounced. A full stripping test is warranted in the case of SLICE. Since 
this project is focused a comparison of alternate hull forms, appendages will be 
excluded, leaving the stripping test to future work. 

2.7 The SLICE Data Scaling Method 

Based on the above discussion, the traditional method of data scaling outlined 
in section 2-1 is recommended for the SLICE project and will be used in this research. 
Calculating the frictional resistance using the ITTC 1957 line will yield a moderately 
conservative answer. Although the ITTC 1957 line was derived from monohull tests, it 



is the only method that naturally accounts for the Froude Number dependence of 
boundary layer development. All other phenomena (wave making, eddy formation and 
spray development) will be treated as part of residuary resistance and scaled 
accordingly. 

It is important to recognize is that this technique of data scaling is considered 
reliable due to the enormous volume of model testing that has been performed on 
numerous series of monohull models. The validation of this process for SLICE will 
only come when the full scale ship is built and tested. This is the real value of the 
ATD. Building and properly testing the full scale ship will significantly contribute to 
the Naval Architecture community's understanding of the robustness of current testing 
and scaling techniques. 

2.8 A Numerical Approach 

Model testing is a time consuming, detail intensive process that requires 
significant capital expenditure. For a new hull form, some degree of model testing is 
mandatory in the early stages of design so observations of physical behavior can be 
made. This is critical to validating assumptions, but from a cost savings perspective, 
the ability to numerically model the ship's behavior is attractive. Small changes to the 
hull can be made and tested without the expense of building a new model. If the 
numerical model can be kept mathematically simple, comparative results can be 
obtained quickly and cheaply There must be a degree of confidence however that the 
computer model fully represents actual resistance. 

The physical phenomena of interest are friction, wave, spray and eddy 
resistance. Friction and spray resistance can be calculated using (2-5), (2-8) and (2-9). 
If there is no significant separation or vortex phenomena, it is reasonable to neglect 
eddy resistance (only observation can answer this) and a prediction for full scale 
resistance can be obtained if a method of solving (2-7) can be found. This is 
significant to the SLICE discussion since it may provide a simple and fundamental 
tool for predicting and minimizing the wave making hump magnitude. 
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A solution to Mitchell’s integral was outlined by Lunde 11 in 1951. The process 
was further applied to a ship called the Wave Cancellation Multi-hull Ship Concept by 
Wilson and Hsu ' They express the total wave making resistance as the summation of 
body and interaction terms as shown: 




( 2 - 10 ) 



where N represents the number of bodies the overall hull can be divided into and: 
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13 Lunde. J K . "On the Linearized Theory of Wave Resistance for Displacement Ships in Steady 
and Accelerated Motion". SNAME Transactions . Yol 59, 1951. pp 25-76 

1-1 Wilson. M B . and Hsu. C C , "Wave Cancellation Multihull Ship Concept", Presented at the 
Intersociety High Performance Marine Vehicle Conference. Arlington, Va . 24-27 June 1992. 
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In these expressions, the like product terms are the body terms and the cross product 
terms are the interaction terms. 

The underlying assumptions to Lunde’s approach are : 

1. The fluid is incompressible, irrotational and inviscid. 

2. The wave height is small compared to wave length. 

3. The ship motion can be reasonably represented by a fine distribution of 
sources and sinks along the centerline of the body. 

Prior observation of SLICE models in the towing tank indicates that conditions 
2 and 3 may not be satisfied ideally. Each of these phenomenon must be watched 
closely in subsequent tests and considered in the interpretation of results. Havelock 15 
demonstrated good correlation using this technique with hull forms having L/B as 
small as 10.6. Data presented by Newman 16 indicates that there is better correlation 
between model test data and thin ship theory at higher values of ZF. SLICE is a full 
form with L.D = 4. 5-5.0 but operates at &"=■ 0.9-1. 2. Although SLICE is not an exact 
fit for thin ship theory, the approach may be reasonable. 



L ' Havelock. T H . "Wave Resistance Theory and its Application to Ship Problems", SNAME 
Transactions . 1951, pp 13-24 

16 Newman. J.N. . op.cit . p282 
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Chapter 3 Model Testing 



As pan of this research, three models were designed, built and tested. Appendix A 
provides detailed lessons learned from this testing experience which may be useful to 
other students attempting to undertake similar model testing projects. 

3.1 Model Test Objectives 

There were three objectives underlying the model testing phase of this project: 

1. Expand the quantity of available model scale resistance data for the 
baseline SLICE configuration. 

2. Visually observe model performance in sufficient detail to 
understand the physical behavior of this hullform As discussed in chapter 2. this is 
important for justifying the scaling technique selected. 

3. Test operationally equivalent hullforms to help assess baseline 
SLICE advantages and disadvantages. 

The models were tested at the United States Naval Academy's 380' tow tank 
facility during the week of 6-10 March 1995. 

3.2 Model Descriptions 

The three models tested as part of this research were designed to meet the 
same operational constraints. This was done to ensure that comparisons made at the 
end of the testing could be made on the basis of similar ships. All models when scaled 
to full scale have a displacement of 177 Lton. a pod diameter of 7-8 feet and a strut 
beam of 3.25 feet. These dimensions are the working dimensions of the ATD vessel, 
driven by requirements other than hydrodynamic performance. The pod diameter is 
controlled by engine size requirements and the strut beam is controlled by access 
trunks into the pods. With these requirements imposed on all three models, they are 
considered equivalent hulls. 

The three models tested were designated M-l. M-2 and M-3. Detailed tables 
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of hull offsets and design calculations for each model are contained in Appendices B. 
C and D respectively. Individual descriptions of each model follow. 

MODEL M-l BASELINE: 

This model utilizes LMSC's offsets as tested in December 1994 with the 
exception that the fillets at the strut/pod interface and the haunch structures are 
excluded. The fillets were excluded to help simplify model building. This modification 
is not considered detrimental to the testing outcome. The combination of the pod, strut 
and fillet in the LMSC model creates a effective L D of 3.9, the omission of the fillets 
increases this ratio to a more favorable 4.2. This modification should result in slightly 
lower resistance values. The bow wedges installed in the December LMSC model are 
replicated in this model. The purpose of this model is twofold. First, it provides a 
means of correlating the data obtained in these tests with the data obtained by LMSC 
in December 1994 This is considered important since the tests were conducted in two 
different facilities, under different conditions with different scale models. The ability 
to correlate data between Model M-l and the LMSC model will be valuable in 
demonstrating repeatability of results. The second use of Model M-l is to expand the 
available range of data for the baseline SLICE. Model tests were conducted at model 
equivalent speeds through 40 knots to ensure that the entire operating range was 
properly represented by data. 

Table 3-1 provides the principle dimensions of Model M-l. Profile and top 
view sketches of the forward and aft strut/pod combinations are in Figure 3-1. In this 
drawing the alignment of the struts with the corresponding pods is correctly 
represented. The spacing between the forward and aft combination is not preserved 
however in the interest of fitting the drawing onto one page. This model was 
constructed with 29 lbf of reserve buoyancy. The struts are manufactured from balsa 
wood with a fiberglass shell, the pods are molded fiberglass. 
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Scale Factor 


17.26:1 


Displacement 


75 Lbf 


Beam 


36.15 inches 


Draft 


8.34 inches 


Pod Diameter 


5.5625 inches 


Fwd Pod Length 


22.9 inches 


Aft Pod Length 


25 inches 


Fwd Pod Separation (CL to CL) 


23.4 inches 


Aft Pod Separation t CL to CL) 


30.25 inches 


Strut Length 


16.68 inches 


Stmt Max Beam 


2.26 inches 



Table 3-1: Model M-l Baseline 
Principle Dimensions 




Figure 3-1: Model M-l Forward and Aft Strut/Pod 

(Note: In this drawing the fwd-aft pod spacing is not preserved) 
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MODEL M-2 RATIONAL DESIGN: 



The purpose of this model is to explore the effects on total resistance of using 
finer offsets for the struts and pods . To avoid changing too many variables in a single 
test, this model maintains all the gross parameters of the Model M-l Baseline (ie: 
Strut/Pod alignment, fore/aft spacing and centerline to centerline spacing). The 
expectation for this design is that by comparing this model's results with that of M-l, 
some conclusions can be drawn about the relative magnitude of both friction and wave 
making resistance through the range of speeds tested. The anticipated result from this 
model is that the amplitude of the wave making hump will be diminished. 

The offsets for the pods are adapted from Series 58 data 1 ’ This series was 
selected since it represents a concerted effort to optimize submarine hull forms. The 
model selected from the senes was Model 4155. This model was selected based on its 
length to diameter (L/D) ratio of 5 which is similar but finer than the baseline SLICE 
L/D of 4 2. It should be noted that within the Series 58 results. Model 4155 falls 
outside the range of optimum L D which was found to be between 6 and 8. Model 
4155 was selected despite this drawback in order to satisfy the operational requirement 
of creating a full scale 177 Lton vessel with pod diameters of 7-8 feet. These criteria 
could not be satisfied with a series 58 hull with a L/'D of 6-8. 

A trait of the Series 58 hull forms is that they have no parallel midbody. 
Although hydrodynatnically this is supenor. it would significantly complicate the 
building of the full scale ship. The Model M-2 therefore incorporates a modification to 
Model 4155 by replacing some of the midships shape with parallel midbodv. 

The strut offsets utilize a symmetric parabolic leading and trailing edge. The 
general parameters of this strut are outlined in Chapter 2-5 but are also constrained by 
the design requirement that the full scale maximum beam must be at least 3.25 feet. 
The calculations used to develop the stnit and pod offsets and a table of the actual 

1 Gertler. M. "Resistance Experiments on a Systematic Senes of Streamlined Bodies of 
Revolution for Application to the Design of High Speed Submarines", David Taylor Model Basin Report 
C-297. unclassified and approved for public release 18 Dec 1*572, April 1950 
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offsets used are contained in Appendix C. 

Model M-2 represents a partial optimization. The offsets for the struts and pods 
are based on research conducted in optimizing these types of shapes. The optimized 
values are then adjusted to conform to design requirements. Model M-2 make no 
attempt to optimize the placement of the strut pod combinations relative to one 
another. Table 3-2 provides the principle dimensions of Model M-2. Profile and top 
view sketches of a the forward and aft strut/pod combinations are at Figure 3-2. This 
model was constructed with 30 lbf of reserve buoyancy. The struts are manufactured 
from balsa wood with a fiberglass shell, the pods are molded fiberglass 

MODEL M-3 OPTIMIZED HIGH SPEED SWATH 

The purpose of this model is to compare SLICE performance to a more 
conventional SWATH. The offsets used for the high speed SWATH were reported by 
McGregor 18 . The demihull diameter was forced to meet the 7-8 foot full scale 
requirement, the strut maximum beam was forced to meet the 3.25 foot full scale 
requirement and the overall displacement was forced to meet the 177 Lton full scale 
requirement established earlier. In order to meet these requirements, the demihull L/D 
was established at 12.86. This corresponds to a full scale demihull length of 90 feet, 
about 15 feet shorter than the Length Overall (LOA) of the proposed ATD SLICE. 

The consequence of this is that although Model M-3 is hydrodynamically equivalent to 
Model M-l. depending on the extent to which the box structure is cantilevered forward 
and aft of the struts, the SWATH will end up with as much as 600 square feet less 
arrangable deck space. This is an 11% decrease from the SLICE ATD vessel. 

A tandem strut arrangement was selected in order to minimize the displacement 
occupied by the struts. This was necessary in order to maximize the demihull 
displacement and maintain a reasonable demihull length. As with model M-2, the 



18 McGregor. R.C. and Chun. H.H., "On the Potennal of SWATH Ships for Very High Speed 
Operation", Fast 91 . pp 491-?06 
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strut offsets are based on the guidance provided in Chapter 2-5. In his SWATH 
design, McGregor provided guidance for the optimum placement of the struts in 
relation to the pods. This placement was utilized in Model M-3. The calculations used 
to develop the strut and pod offsets and a table of offsets are contained in Appendix 
D. Table 3-3 provides the principle dimensions of Model M-3. A profile and top view 
sketch of a the demihull is at Figure 3-3. The struts of this model are manufactured 
from balsa wood with a fiberglass shell, the demihulls are hollowed PVC rod, shaped 
to the proper offsets. There is 10 lbf of reserve buoyancy. 

The reason the three models were not all built to the same scale factor stems 
from the actual building process. There is a slight variation in the full scale diameter 
of the three underwater hull diameters. In full scale. M-l is 8 feet. M-2 is 7.74 feet 
and M-3 is 7 feet. PVC pipe was used to mold parallel midbody sections. This 
material is relatively inexpensive and comes in a standard pipe outer diameter 
dimension of 5.5625 inches. This dimension became the driver for determining the 
model scale factor. Since the three models have different full scale diameters, the scale 
factors vary accordingly. 
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Scale Factor 


16.68:1 


Displacement 


83.3 Ibf 


Beam 


37.41 inches 


Draft 


8.63 inches 


Pod Diameter 


5.5625 inches 


Fwd Pod Length 


27.8 inches 


Aft Pod Length 


27.8 inches 


Fwd Pod Separation (CL to CL) 


24.17 inches 


Aft Pod Separation (CL to CL) 


36.26 inches 


Strut Length 


17.25 inches 


Strut Max Beam 


2.34 inches 



Table 3-2: Model M-2 Modified Baseline 
Principle Dimensions 




Figure 3-2: Model M-2 Forward and Aft Strut/Pod 

(Note. !n this drawing the fwd-aft pod spacing is not preserved) 
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Scale Factor 


15.10.1 


Displacement 


115.1 lbf 1 


Beam 


39.93 inches 


Draft 


1 1.13 inches 


Demihull Diameter 


5.5625 inches 


Demihull Length 


71.5 inches 


Strut Length 


19 10 inches 


Strut Max Beam 


2.58 inches 



Table 3-3: Model M-3 McGregor's High Speed SWATH 
Principle Dimensions 




Figure 3-3: Model M-3 Demihull and Strut 
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3.3 Test Description 

A copy of the test plan developed for the Naval Academy Tow Tank tests is 
included as Appendix E. This enclosure explains the planned tests in detail. 

Two basic tests were planned for the models. The first measures resistance and 
lift with the models fixed in pitch and heave. This test replicates the test configuration 
used by LMSC thereby allowing for correlation with the LMSC data. Data collection 
runs consisted of the model held at its design draft with no trim. At the slow end of 
the speed spectrum, one speed change was allowed during data collection. At higher 
speeds only one speed was used per test run. Data was collected at full scale 
equivalent speeds from 3-40 knots with the following model configurations: 

1) Model M-l Baseline in the LMSC configuration. 

2) Model M-l Baseline with the forward and aft struts in line. 

3) Model M-l Baseline in the LMSC configuration at varying drafts. 

4) Model M-2 Modified Baseline in the LMSC configuration. 

5) Model M-3 McGregor's SWATH. 

The second test, quasi free in heave, was not conducted. The intent of this test 
was to replicate the effects of control surfaces through a counterweight system with 
the hope that information could be gained regarding the necessary sensitivity of the 
control system. Conceptually, the counterweight would provide the same lift as the 
control surfaces without added model building complexity. The model would then be 
free to respond to small perturbations and its stability qualitatively assessed. This test 
was considered a second priority to the project and time in the test facility ran out 
before it could be accomplished. It is considered a worthwhile area of further study. 

The test plan indicates that a fourth model, the Midfoil, was tested. This model 
was built and tested at the request of ONR but is not considered part of this research. 
It is not reported in this thesis. 
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3.4 Physical Observations 

The importance of observing the behavior of the models can not be overstated. 
A videotape of the model tests is included as part of this thesis. Shown in the tape is a 
port beam perspective, both above and underwater as well as a stem view of all the 
models tested. The video is organized by speed increment, showing each model first at 
the peak of the wave making hump (10 and 12 knot full scale equivalent) then in the 
post hump valley (16 knot full scale equivalent) and finally at three high speed 
increments (25, 30 and 35 knot full scale equivalents). The following qualitative 
observations are made regarding the performance of the models. 

PROPELLER SUBMERGENCE 

The ATD vessel will have propellers mounted on the aft end of the forward 
pods. After the August LMSC model tests, there was significant concern that the wave 
system created in the hump was large enough to expose the propeller. This evaluation 
was made by observing the wake from above the free surface but was difficult to 
assess properly due to spray and distortion. With M-l, the underwater camera used in 
this project reveals that there is better submergence at the propeller than apparent 
when viewed from above. 

HUMP SPEED PERFORMANCE 

It was expected that M-2 with its finer lines would exhibit less wave making 
than the very bluff M-l in the hump speed range. In fact, the surprising result was 
that the wakes produced were nearly identical. The wake of M-3 at the same speed 
was much smaller than either M-l or M-2. The implications of this observation are 
unclear without examining actual data. The similarity' between the M-l and M-2 wakes 
would suggest that the form of the bodies, both strut and pod. is not a controlling 
factor in hump speed resistance. The vast difference between M-l /M-2 and M-3 
suggests that the underwater body plays a role in wave cancellation. 

The correct conclusion may simply be that this is a demonstration of the flow 
differences of bluff bodies (M-l L/D = 4.2, M-2 L/D = 5.0) and more slender bodies 
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(M-3 L/D = 12.9) or there may be some more complicated mechanism of wave 
cancellation that needs to be explored. 

HIGH SPEED PERFORMANCE 

In the full scale range of 25-40 knots, there is a significant difference in the 
performance of M-l and M-2. At these speeds, M-l creates a very large vortex at the 
aft intersection of the strut and the pod. Above the surface, this is manifest as a 
rooster tail. As the speed increases, the inception point of this vortex travels further aft 
along the pod tail and at the full scale equivalent of 40 knots almost reaches the end 
of the pod. This effect is absent in M-2 and M-3. 

This effect does not appear to be caused by the relatively small scale of the 
models. Reviewing footage of the LMSC tests, the same phenomenon is apparent. It is 
more pronounced in the August tests than the December tests, but it is difficult to tell 
if that can be attributed to the design changes LMSC made. There is not enough 
quality video footage of the December tests to draw this conclusion. It is possible that 
the August and December tests did not allow a full steady state condition to be 
established. In both the August and December tests, speed increments were maintained 
on average for 10 seconds measured from the start of a speed change to the start of 
the subsequent speed change. As observ ed in the March tests, this amount of time is 
barely adequate to establish steady state conditions. When accelerated from rest, the 
March data shows that forces did not stabilize until 8-10 seconds after the steady state 
speed was achieved. The LMSC testing may have stabilized sooner since the 
acceleration was only between speed increments and not from rest. Nonetheless, the 
data collection period was extremely short. The inception speed of the rooster tail 
appears to be comparable in both the LMSC tests and this project's tests. This 
indicates that the effect follows Froude scaling. 

After the August LMSC tests, much attention was focused on scaling spray 
effects. This was in response to the extreme wake created at high speed. LMSC's 
position was that spray effects could be ignored in full scale predictions and designed 
out using spray rails. The underwater observation of M-l demonstrates that much of 
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what was called spray in August was actually produced by an entirely different 
mechanism. The vortex that creates the rooster tail is a separation effect created as the 
flow unsuccessfully attempts to conform to the hullform. When M-l, M-2 and M-3 are 
compared at the same speed, it becomes obvious that only a very little of the splashing 
seen above the surface in M-l can be attributed to spray generated by the struts. The 
majority of the splashing effect is more accurately accounted for as a by-product of 
eddy resistance. 

WAKE CHARACTERISTICS 

The wake produced by the SLICE hullform is compact and has a greater 
amplitude than that of the SWATH. The SWATH produces a relatively flat but 
significantly longer wake. From an energy in the water standpoint, comparison 
becomes difficult since the entire volume of the wake must be considered. 

3.5 Test Data 

Raw data collected during the testing is presented in Appendix F. Test 1 
results (fixed in pitch and heave measurements of drag and lift) are presented 
graphically in. Drag ( positive x) is taken to be opposite the direction of forward 
motion. Sinkage (positive z) is down from the free surface. Several speed settings 
were repeated during the testing. Repeated speeds are recorded in these figures and are 
designated Dual Point (DP) No direct comparison of performance using Figures 
(3-4) - (3-6)should be done since each model is built to a slightly different scale 
factor . 

Figures (3-7) and (3-8) present model scale data at the same scale factor. This 
required scaling M-2 and M-3 data to the M-l scale. Data scaling was accomplished 
using the process outlined in Chapter 2-7. Scaling between models was accomplished 
using the composite scale factors shown: 
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(3-1) 



^ = Knodel2 _ 16.68 

2 " Ko*en~ 17.26 



=0.966 



^ _ ^model 3 _ 15.10 
3 ~ Uii" 17 -26 



= 0 . 875 



(3-2) 



The calculations used in this process are in Appendix G. 

Figure (3-9) shows the data collected for Model M-l tested in both the 
Baseline configuration and an in-line configuration of the forward and aft strut/pod 
groups. The results of testing Model M-l at varying depths was inconclusive. 
Insufficient data was collected to develop a trend. 



46 



40 

35 

30 

25 

20 

15 

10 

5 

0 

-5 



□ □ Sinkage Force 

G G Drag Force 



DP indicates a repeated data point. Both 




5 10 15 20 




Velocity (fps) 



Figure 3-4: Model M-l Fixed in Pitch and Heave 
Unsealed Data 
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Figure 3-5: Model M-2 Fixed in Pitch and Heave 
Unsealed Data 
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Figure 3-6: Model M-3 Fixed in Pitch and Heave 
Unsealed Data 
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data at 17.26:1 scale 




Velocity (fps) 



Figure 3-7: Models M-l/M-2 and M-3 Fixed in Pitch and Heave 
Data Presented in 17.26:1 Scale 
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data at 17.26:1 scale 




Velocity (fps) 



Figure 3-8: Models M-l/M-2 and M-3 Fixed in Pitch and Heave 
Data Presented in 17.26:1 Scale 
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Figure 3-9: Models M-l Fixed in Pitch and Heave 
Tested in the Baseline and In-Line Configurations 
Data Presented in 17.26:1 Scale 
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Several comments can be made about this data. 



REPEATABILITY OF RESULTS 

A large number of data points were repeated through the course of the testing. 
Without exception, the difference in recorded data is statistically insignificant. Most 
repeated data was recorded in successive runs down the tank. A notable exception to 
this are runs 70 and 71 in Appendix F. In this case, testing of Model M-l was 
completed. Model M-2 and Model M-3 were tested, then Model M-l was remounted 
and realigned on the carriage and five data points collected. There was excellent 
correlation between these points and data collected two days earlier. 

HUMP SPEED MODEL PERFORMANCE 

The resistance values for M-l and M-2 through the hump, Figure 3-7, support 
the observation made earlier that the wakes are very similar. The similarity in 
resistance values was an unanticipated result. M-2 represents the extent to which L/D 
can be extended and still retain the features of SLICE. Since this variation was 
insufficient to impact the amplitude of the hump, the only other alternative for 
diminishing the hump magnitude is altering the relative position of struts and pods. 

The fact that the hump for M-l, M-2 and M-3 all occur at essentially the same 
speed indicates that the Strut Length is the appropriate length to characterize SLICE. 
This conclusion is reached by considering M-3 in contrast to M-l and M-2. The hump 
is expected to occur at 27= 0.45. For the demihull of M-3 this corresponds to 
V= 6.25 fps. for the struts it corresponds to V= 3.2 fps. The hump actually occurs at 
V=4.5, 

demonstrating dominance by the struts. The pods of M-l and M-2 are less than half 
the length of M-3's demihulls and yet the hump for all three is at the same speed. The 
struts for all three models were essentially the same length. The underwater bodies 
may be contributing to wave making through cancellation effects but the position of 
the hump is driven by the surface piercing bodies. 
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HIGH SPEED MODEL PERFORMANCE 

In the range 6.5-10 fps (full scale equivalent: 16-25 knots), M-l total resistance 
is up to 20% greater than that of M-2 even though M-l has less wetted surface and 
therefore less frictional resistance than M-2. This is the first indication that model 
form must be traded off with friction to optimize post hump resistance. 

For speeds greater than 10 fps. the rapid increase in M-l resistance can be 
directly correlated with the appearance of the trailing edge vortex noted in the physical 
observations. 

M-l LMSC COMPARISON 

One test objective for M-l was to compare M-l results with the December 
LMSC model. Figure 3-10 presents the M-l data as collected as well as a sample of 
LMSC data from both August and December scaled to 17.26:1. Model M-2 results are 
also presented scaled accordingly. The August LMSC curve represents LMSC data 
runs 106-121, The December LMSC curve represents LMSC data runs 86-97. These 
models were tested with appendages and haunch structures in place. 

M-l correlation with the December tests is good below 8.25 fps (full scale: 20 
knots), becomes marginal from 8.25-10 fps (full scale:20-25 knots), and does not exist 
above this value. M-l correlation is more satisfactory with the August test results 
The lack of high speed correlation requires comment. 

Of concern is the vortex created at 10 fps (full scale; 25 knots). The 
phenomenon is documented as existing in the August and December tests as well as 
the tests conducted on M-l. The question is deciding if the lack of correlation 
indicates that effect does not scale, or if there is some other mechanism that is causing 
the difference. The possibility that the LMSC data may have been obtained before 
steady state was achieved has already been discussed. Although this is a fault of the 
LMSC test program, it is difficult to attribute the entire difference in test results on 
this shortcoming alone. 



Drag Force (Ibf) 




Velocity in FPS 



Figure 3-10: Models \I-1 Performance Compared With LMSC Results 
Data Presented in 17.26:1 Scale 
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In the tests of M-l. there is direct impingement of the rooster tail on the model 
frame. During the testing, this was considered acceptable for two reasons: 1) 
Impingement also occurred in the LMSC model tests and will occur in the full scale 
ship depending on the constniction of the box structure. 2) During the testing, the 
impingement was considered a minor effect since the profile of the frame structure 
was small (60 in : total). Reviewing video footage of these tests, the direction of 
rooster tail flow was perpendicular to the frame, thereby maximizing the force of 
impingement. In the LMSC tests impingement was not as direct. Appendix H 
provides a calculation that demonstrates that impingement could account for up to 
11 pounds of the measured force at 15 fps. 

There are enough differences in test configurations to preclude a definitive 
explanation for the lack of correlation. The rooster tail impingement may contribute to 
the lack of correlation. The conclusion drawn at this point is that the trailing vortex is 
real and will scale but may not directly contribute to overall resistance to the extent 
indicated by the Model M-l results in Figure (3-7). Failure to achieve steady state in 
the LMSC test may also be a significant factor in this difference. Re-testing M-l with 
the frame isolated will help eliminate this ambiguity but ultimately full scale test 
results are required 

3.6 Full Scale Predictions and Discussion 

With the understanding that the M-l results may over predict full scale 
resistance, figures (3-1 1) and (3-12) present full scale data for all three models. Figure 
(3-13) provides power requirements based on these values. Data scaling was 
accomplished using the process outlined in Chapter 2-7. The calculations used in this 
process are at Appendix G. Based on figure (3-1 1), the high speed SWATH offers the 
best resistance performance through the entire speed range. 

Figures (3-14)- (3-19) present M-l. M-2 and M-3 data in model and full scale, 
decomposed into the frictional, residuary and correlation factor components. In the 
post hump region of the resistance curve, residuary (ie wave making) resistance 
accounts for at least 60% of the total resistance value. This further supports the 
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concept that the hull form must be designed as a compromise between frictional and 
wave making resistance, not a single component optimization. 
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Figure 3-11: Models M-l/M-2 and M-3 Fixed in Pitch and Heave 
Full Scale Resistance Data 
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Figure 3-12: Models M-l/M-2 and M-3 Fixed in Pitch and Heave 
Full Scale Resistance Data 
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Figure 3-13: Models M-l/M-2 and M-3 Fixed in Pitch and Heave 

Powering Data 
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Figure 3-14: Model M-l Resistance Components as a Percent of Total Resistance 

Derived from I'nscaled Test Data 
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Figure 3-15: Model M-l Resistance Components as a Percent of Total Resistance 

Based on Full Scale Results 
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Figure 3-16: Model M-2 Resistance Components as a Percent of Total Resistance 

Derived from Unsealed Test Data 
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Figure 3-17: Model M-2 Resistance Components as a Percent of Total Resistance 

Based on Full Scale Results 
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Figure 3-18: Model M-3 Resistance Components as a Percent of Total Resistance 

Derived from Unsealed Test Data 
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Figure 3*19: Model M-3 Resistance Components as a Percent of Total Resistance 

Based on Full Scale Results 
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Chapter 4 Numerical Solution 



4.1 Objective of a Numerical Approach 

The ability to predict a ship's resistance using numerical methods has many 
advantages. If an appropriate method can be found, the expense and time of model 
testing may be reduced. Current state of the science is approaching realistic treatment 
of the free surface (ie: non linearized theory) but these methods are still too expensive 
to use as an iterative design tool. Simpler less expensive methods may not be accurate 
enough to exactly predict resistance but can be extremely useful for comparative 
analysis of small changes to the same hullform. These methods require tow tank 
model testing as verification to ensure that all physical phenomena are properly 
represented. 

fhis thesis attempts to validate the SLICE applicability of a numerical method 
solution of Mitchell's integral as an approximation for wave resistance. As cited in 
Chapter 2. the approach, first outlined by Lunde. was used by Wilson in the Wave 
Cancellation Multihull Ship concept. It is not expected that this thin ship 
approximation will provide exact resistance values. The goal of this chapter is to 
develop an understanding of the body interactions in the SLICE hullform to assist 
resistance optimization through body placement. Of specific interest is reducing the 
magnitude of the wave making hump. The thin ship approach is useful since it allows 
wave making resistance to be segregated into all component body and interaction 
terms. This is important while studying the effect of small geometry changes on total 
resistance. 

Model M-2 is selected as the subject for this stud} - . Physical observation during 
the model tests demonstrate that this model has a reasonable chance of satisfying thin 
ship flow requirements. 
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4.2 Solution Technique 

Restated from Chapter 2.8. the equation to be solved is: 



r N N N N 

R w = 1671 p>c* f 

v i 1 1 J 1 1 J 



(4-1) 



where: 



TL 

P=J r JoU^)cos[KCtro^+3w'ne)sec : 0]e ,wc ‘ 6 fiwt (4-la) 

0 0 



TL 

Q=J J‘o(x r :)sm[K(xcosQ+ysinG)se(r6]e^ c ^dxdz (4-lb) 

00 



0 = ~X“(/C*^)) ( 4 ‘ lc ) 

2 71 bx 

k=-£- (4-ld) 

V 2 

The coordinate system that will be used has its origin in the undisturbed free surface 
with f-X in the direction of ship motion. - i -Z vertically up from the free surface and 
+Y in accordance with the right hand rule. 

DEVELOPMENT OF THE BODY AND POTENTIAL FUNCTIONS 

The first requirement is to write a function that describes the strut and pod 
forms. This function is written with the bodies at the global origin. Later, these bodies 
will be translated to their actual relative positions through the X, Y and Z terms in 
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Equation (4-la) and (4-lb). 

For M-2. the strut form was developed using an identical parabolic nose and 
tail shape with parallel midbody in between. The equations used to describe the strut 
shapes are: 

Nose Offsets - — * [ 1 (4-2) 

2 



MidbodyOffsets = ^ 
z 



(4-3) 



TailOffsets=— *[1 

2 i- x 



(4-4) 



where for this specific strut. 

B (the maximum beam) = 0.195 ft 
L n (nose length) = L, (tail length) = 0.51 ft 
Lp m b (parallel midbody length) = 0.42 ft 



To reduce the complexity of the strut bod} - potential function as applied in (4- la) and 
(4-lb), the offsets calculated by (4-2), (4-3) and (4-4) were curve fit to a single high 
order polynomial given by: 



Sfra =0.4893 ** 6 +2.1 126 ** 5 +3. 1328 ** 4 + 1.7248 3 -0.0786 -0.3890 **+0.0008 

(4-5) 
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The source distribution representing the strut is written as: 



o 






— * ( 2.936 ** 3 + 10.563 ** 4 + 12 . 53 1 ** 3 + 5. 174 ** 2 - 0. 1572 **- 0 . 389 ) 
2n 



The pod shapes are created in a similar fashion. Now however, the nose is 
represented by an ellipse and the tail by a parabola. The pod offsets are given by the 
following: 






PodNose = — * [1 -(— — ) ] 

2 L n 



(4-7) 



PodMidbody=^ 



(4-8) 



D -L_,, n *u 

PodTail=-*[ l-( 2— S5) ] 



(4-9) 



where lor this specific pod. 

D (the maximum diameter) = 0.4635 ft 

L 0 (nose length) = 0.67 ft 

L, (tail length) = 1.1 ft 

Lpmb (parallel midbodv lenath) = 0.60 ft 

IU. = 2.00 
- 2.40 



Again a single high order polynomial was generated to represent the pod and is given 
by: 
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Pod= -0.04 *x 7 -0.32 ** 6 - 1. 16 -2.37 *x 4 -2.96 *jc 3 -2.27 *;r -0.99 **+0.038 



The source distribution representing the pod is written as: 



= -— * ( -0.253 -1.93 ** 5 -5.82 ** 4 -9.49 ** 3 -8.88 ** 2 -4.54 **-0.99 

m 2n 

The pods are treated as a line source in (4-1 a) and (4- lb) by fixing the z term 
in the exponential as a constant equal to the strut depth plus the maximum radius of 
the pod. This follows the example provided by Chapman 19 for dealing with submerged 
bodies of revolution. 

PAIRING OF TERMS 

Struts and pods are numbered in the traditional naval manner. The scheme is 
presented in Figure (4-1). There are a total of 36 possible combinations of the eight 
bodies that form the SLICE hull. Eight combinations are "like" or body terms and 28 
combinations are different or "cross product" terms. These terms are presented below 
grouped by relative position in the hull. The numbers indicate which bodies are 
interacting in the term calculated: 

1) Body Terms: 

1-1. 2-2. 3-3. 4-4. 5-5, 6-6. 7-7. 8-8 

2) Same Comer Cross Terms: 



19 Chapman. R.B., Hydrodynamic Drag of Semisubmerged Ships, American Society of Mechanical Engineers, 
paper number "2-Wa, 1972. pp 879-884. 



71 



1-5. 2-6. 3-7. 4-8 

3) Longitudinal Cross Terms: 

1-3. 1-7. 5-3. 5-7. 2-4. 2-8. 6-4. 6-8 

4) Diagonal Cross Terms: 

1-4. 1-8. 5-4. 5-8. 2-3. 2-7. 6-3. 6-7 

5) Transverse Cross Terms: 

1-2. 1-6. 5-2. 5-6. 3-4. 3-8. 7-4. 7-8 

This physical grouping is retained through the calculations to simplify' comparisons. 




Figure 4-1: SLICE Body Numbering Scheme 
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CODING 



MATLAB was used for the required computations. This choice was made to 
utilize the installed graphics capabilities of the program. Plotting results as 
computations were made was useful both in validating geometry assumptions and 
debugging. A logic flow chart of the overall computation scheme is at Appendix I. 

The routines used are also presented in this appendix. The pod and strut geometries 
were developed in the files PODGEO.m and STRUTGEO.m, The main wave making 
resistance routine is contained in the file Q.m. The solution to equations (4-1), (4-la) 
and (4-lb) are solved in called subroutines for each of the thirty six paired terms in 
files named for the bodies that are paired. Samples of a strut body, pod body and 
strut'pod cross term are contained in the files ONEONE.m. SIXSLX.m and 
ONESIX.m. All other body and cross product subroutines are patterned after these 
files. This apparently cumbersome organization had one significant advantage. Once 
the files were written and assembles, bodies could be toggled on and off in the 
computation simply by adding a "°o" at the head of the appropriate subroutine call 
line. This allows the overall SLICE geometry to be reduced into smaller pieces rather 
easily and helped demonstrate the complex interactions that occur in the full 
multibodied hull. 

4.3 Results 

Only partial results were obtained for the numerical solution during the 
allocated period of study for this thesis. Complete validation of results was not 
obtained and the results contained in this section require further analysis before design 
decisions can be made. Sufficient progress was made to draw some preliminary 
conclusions. Model M-2 was computed in the speed range of 3-9 fps (full scale 11-22 
knots) using the values shown in Table (4-1). The results presented in Figure (4-2) 
and Table (4-2) are based on using fifty subdivisions in the x integration, 12 divisions 
in the z integration and 10 divisions in the 0 integration. Testing was conducted 
varying the number of divisions and this was found to be a minimum acceptable limit 
with the integrity of the result maintained. 
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Strut Length 


1 .44 ft 


Strut Beam (max) 


0 195 ft 


Strut Height 


0.36 ft 


Fwd Strut Set Back (x axis reference) 


0.0 ft 


Fwd Strut Offset (y axis reference) 


± 0.979 ft 


, Aft Strut Set Back (x axis reference) 


-3.48 ft 


Aft Strut Offset (\ axis reference) 


± 1 .396 ft 

1 


Strut Submergence iz axis reference) 


-0 1 8 ft 


Pod Length 


2.321 ft 


Pod Diameter (max) 


0.4635 ft 


Fwd Pod Setback (x axis reference) 


0.0 ft 


Fwd Pod Offset (y axis reference) 


± 0.979 ft 


Aft Pod Set Back (x axis reference) 


-3.16 ft 


Aft Pod Offset (y axis reference) 


± 1.396 ft 


Pod Submergence <z axis reference) 


-0.59 ft 



Table 4-1: M-2 Gross Parameters Used For Calculation 

Figure 4-2 presents only the wavemaking portion of resistance as predicted by 
the numerical scheme. The negative overall result is attributed to the overall sign 
convention that is used in the calculation which has x positive in the direction of 
motion. An unresolved problem with this conclusion is that the body terms (1-1 
through 8-8) should have been negative on their own. This however was not the case. 
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K (lb!) 




Figure 4-2: Wave Making Resistance of SLICE calculated from Thin Ship Theory 
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Speed 


l-i 


2-2 


3-3 


4-4 


5-5 


6-6 


7-7 


8-8 


(fps) 


(lbf) 


(lbf) 


(lbf) 


(lbf) 


(lbf) 


(lbf) 


(lbf) 


(lbf) 


3 


0.39E-3 


0.39E-3 


0.39E-3 


0.39E-3 


0.0004 


0.0004 


0.0004 


0.0004 


4 


0.56E-3 


0.56E-3 


0.56E-3 


0.56E-3 


0.0020 


0.0020 


0.0020 


0.0020 


s 


0.62E-3 


0.62E-3 


0.62E-3 


0.62E-3 


0.0038 


0.0038 


0.0038 


0.0038 


6 


0.65E-3 


0.65E-3 


0.65E-3 


0.65E-3 


0.0052 


0.0052 


0.0052 


0.0052 




0.65E-3 


0.65E-3 


0.65E-3 


0.65E-3 


0.0061 


0.0061 


0.0061 


0.0061 


8 


0.64E-3 


0.64E-3 


0.64E-3 


0.64E-3 


0.0067 


0.006~ 


0.0067 


0.0067 


9 


0.63E-3 


0.63E-3 


0.63E-3 


0.63E-3 


0 . 00^0 


0.0070 


0.0070 


0.0070 



Table 4-2a: Wave Making Resistance Component Terms 

BODY TERMS 



Speed 

(fps) 


1-5 

(lbf) 


2-6 

(lbf) 


3-7 

(lbf) 


4-8 

(lbf) 


" 1 

3 


-0.556 


-0.556 


-0.181 


-0.181 


4 


-1.413 


-1.413 


-1.037 


-1.037 


5 


-2.043 


-2.043 


-1.773 


-1.773 


6 


-2.406 


-2.406 


-2.220 


-2.220 


- 


-2.598 


-2.598 


-2.466 


-2.466 ' 


8 


-2.~00 


-2.700 


-2.601 


-2.601 


9 


-2.751 


-2.751 


-2.675 


-2.675 


Tab 


e 4-2b: Wave Making Resistance Component Terms 



SAME CORNER CROSS TERMS 
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Speed 

(fps) 


1-3 

(lbf) 


1-7 

(lbf) 


5-3 

(lbf) 


5-7 

(lbf) 


2-4 

(lbf) 


2-8 

(lbf) 


6-4 

(lbf) 


6-8 

(lbf) 


3 


0 276 


-0 249 


- 0.245 


0.001 


0.421 


- 0.241 


- 0.436 


0.0 


4 


0.11 


- 0 . 54 " 


- 0.103 


0.002 


0.141 


- 1.054 


- 0.488 


0.004 


5 


0.356 


- 0.135 


- 0.528 


0.0 


0.350 


0.383 


- 0.371 


- 0.003 


6 


- 0.264 


1.128 


0.922 


- 0.007 


- 0.422 


2.035 


1.709 


- 0.010 


- 


- 0.412 


1.325 


1.248 


- 0.008 


- 0.882 


1.879 


2.059 


- 0.008 


8 


- 0.388 


1.142 


1.264 


- 0.006 


- 0.903 


1.106 


1.527 


- 0.004 


9 


- 0.103 


0.682 


0.811 


- 0.004 


- 0.626 


0 302 


0.799 


0.0 


Table 4-2c: W 
LO> 


ave Making Resistance Component ' 
sGITUDLNAL CROSS TERMS 


ferms 


Speed 

( fps ) 


1-4 

(lbf) . 


1-8 

(lbf) 


5-4 

(lbf) 


5-8 

(lbf) 


2-3 

(lbf) 


2-7 

(lbf) 


6-3 

(lbf) 


6-7 

(lbf) 


3 


0.04 E -3 


1 0.119 


- 0.018 


0.0 


0 . 02 E -3 


-0 043 


- 0.018 


0.0 


4 


0 . 33 E -3 


- 0.250 


- 0.574 


0.001 


-.11 E -3 


- 0.045 


- 0.57 


0.001 


s 


- 58 E -3 


0.986 

1 


1.090 


- 0.005 


- 09 E -3 


- 0.016 


1.089 


0.001 


6 


-. 29 E -3 


0.401 


0.815 


- 0.003 


-. 30 E -3 


0.870 


0.815 


- 0.004 


- 


0 . 08 E -3 


- 0 . 44 " 


- 0.241 


0.002 


-. 09 E -3 


0.540 


- 0.242 


- 0.004 


8 


0 . 75 E -3 


- 1.308 


- 1.240 


0.006 


-.31 E -3 


0.986 


- 1.240 


- 0.005 


9 


0 . 58 E -3 


- 1.711 


- 1.580 


0.009 


0.01 E -3 


0.072 


- 1.380 


- 0.002 


Table 4-2d: Wave Male 


ing Resistance Component r 


>rms 



DIAGON AL CROSS TERMS 
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Speed 


1-2 


1-6 


5-2 


5-6 


3-4 


3-8 


7-4 


7-8 


(fps) 


(lbf) 


(lbf) 


(lbf) 


(lbf) 


(lbf) 


(lbf) 


(lbf) 


| 

dbf) ; 


3 


0 0 


-0 079 


-0 079 


0 0 


0 0 


0.06 

; 


-0 115 

' 


0.0 


4 


-.03E-3 


-0.213 


-0.213 


0.001 


0.24E-3 


0.112 


-0.403 


0.001 


5 


0.25E-3 


-0.492 


-0.492 


0.003 


-.23E-3 


0.131 


-0.383 


0.001 


6 


0.13E-3 


-0.484 


-0.484 


0.003 


0.28E-3 


-0.406 


-0.489 


0.003 


- 


0.52E-3 


-0.840 


-0.840 


0.005 


0.20E-3 


-0.048 


-0.541 


0.003 


8 


-.25E-3 


-0.404 


-0 404 


0.005 


0.34E-3 


-0.491 


-0.944 


0.004 


9 


0.07E-3 


-0.4~2 


-0.4^2 


0.005 


0.16E-3 


-0.612 


-0.417 


0.004 



Table 4-2e: 



Wave Making Resistance Component 



erms 



TRANSVERSE CROSS TERMS 



Understanding that the code used is only partially validated, the following 
observations are made: 



STRUTS 

The contribution made by the struts is small. This is true both of the body 
terms and the cross terms between any two struts. This is reasonable for the strut 
geometry used. There is symmetry about midships, the nose and tail sections are 
slender and the overall LD is 7.4. To verify that the low values obtained were not a 
code problem, an alternate geometry was selected that did not possess midships 
symmetry but retained the same LD. A significant increase in resistance values was 
noted. These results follow the conclusions drawn by Chapman and Zhu Bing-quan on 
strut optimization cited in Chapter 2.5. 

CROSS TERMS 

As anticipated, the Same Comer Cross terms have the largest impact on wave 
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making resistance. This is followed in importance by the Longitudinal Cross Terms. 
The Diagonal and Transverse Cross Terms had more significant impact than was 
anticipated. Although smaller than the other terms, it would not be correct to eliminate 
the Transverse terms from optimizing iterations. Also clear from a review of Table 
(4-2) is the complexity of the body interactions. At a given speed, the transverse term 
for a specific body may have an additive effect while the longitudinal term for that 
same body has a subtracting term etc. Additionally, due to the lack of fore aft 
symmetry in the baseline SLICE layout similar body groupings do not behave the 
same at a given speed. Note the difference between 3-8 and 7-4 in the Transverse 
Terms Table as an example of this. This decomposition of SLICE wavemaking 
resistance demonstrates that optimization work must be conducted with the entire 
hullform present. Attempting to sub-optimize a component body (such as a single strut 
pod combination) without observing the impact on the entire hullform may lead to 
erroneous conclusions. 
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Chapter 5 Conclusions 



5.1 Restatement of Thesis Objectives 

This thesis attempted to add to the developing body of SLICE knowledge by 
focusing on the resistance and powering aspects of the SLICE vessel in an 
academically rigorous manner. Three specific objectives were identified: 

1) Independent research into the SLICE hullform to give both the US Navy 
and LMSC a second hydrodynamic viewpoint of the 1 77 Lton vessel prior to the 
ATD. 

2) Development of a larger body of data to aid in optimizing the 177 Lton 
hull form. 

3) A broader understanding of the hydrodynamic features of the hull to 
facilitate deciding the merits of pursuing designs in other displacements more practical 
to the US Navy. 

In order to meet these objectives, several questions required answers: 

1) LMSC theorizes that the appropriate characteristic length for SLICE is the 
length of a single pod. Is this valid? 

2) The specific interaction between bodies has not been defined by LMSC. 
What is the relative importance of form; strut pod alignment; and the effects that a 
strut/pod combination have on any of the other three strul pod combinations that form 
the underwater body? Understanding these effects is necessary for design optimization. 

3) The LMSC proposed design places the forward struts inboard of the aft 
struts. Is the outward lateral force (y axis lift) created by an effective angle of attack 
on the aft struts significant enough to create structural concerns? 

4) The LMSC model tests were all conducted with models fixed in pitch and 
heave. Is this an accurate enough representation of reality to make predictions? 

5) The LMSC tests maintained speed increments for an extremely short 
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duration (usually less than 10 seconds) and frequently truncated test runs at a model 
scale equivalent of 25 knots full scale. Was this sufficient to characterize a ship with 
an intended speed of 30+ knots? 

6) Is the apparent focus on skin friction reduction, as reflected in the body 
forms selected, reasonable? 

7) Does SLICE offer the stated hydrodynamic advantage over a SWATH ship. 
5.2 Conclusions 

CHARACTERISTIC LENGTH 

The appropriate characteristic length for SLICE is the strut length, not the pod 
length. This is illustrated in the resistance curves for M-l.M-2 and M-3. These models 
had significantly different underwater bodies but very similar strut lengths. On a 
resistance vs speed plot, the position of the wave making hump for all three hullforms 
is nearly identical. In order to align these with 4.5. strut length must be chosen as 
the characteristic length used in the definition 

THE IMPORTANCE OF BODY FORM AND BODY INTERACTION 

Model M-2 was developed with the expectation that finer lines would help 
improve the hump speed performance of the baseline SLICE hullform. The offsets 
selected for the pods increased the constructive LD of the pods from 4.1 in the 
baseline to 5.0 in M-2. This was the maximum improvement that could be made to 
L/D and still maintain the gross characteristics of the Baseline SLICE. This change in 
LY) w as not sufficient to impact the hump speed performance. If the SLICE hullform 
is maintained, improvements in hump speed resistance must come from body 
interaction. Underwater photography showed that propeller submergence is better 
through the hump speed than previously thought. 

The finer lines of M-2 did make a significant improvement in the resistance 
measured in the high speed range. Underwater photography revealed that the bluff 
lines of M-l causes the shedding of a significant vortex at the trailing edge of the 
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struts starting at model scale equivalents of 25 knots. These vortices were incorrectly 
interpreted as spray phenomenon in previous testing. In the full scale ATD ship, the 
flow modification caused by the propellers may cause these vortices to be ingested 
into the propellers with a significant loss of thrust. 

MODEL TEST CONDITIONS 

Attempting to test the SLICE hullform in a condition other than fixed in pitch 
and heave introduces a significant level of complexity to the model design. The added 
complexity does not necessarily make the results more useful. The difficulty lies in 
developing a system that connects the model to the tow carriage through the line of 
force. With four underwater bodies all towed from a central point, this is extremely 
complex to manufacture. Any connection point somewhere other than in the line of 
force is easier to manufacture but creates variable applied moment. This effect cannot 
simply be subtracted out since it is impossible to separate it from the effects of real 
hydrodynamic forces. The remaining option then is to lock the model in pitch and 
heave. For the purposes of conducting comparative testing, this simplification is 
justified and was a reasonable simplification made by LMSC. 

The tests conducted as part of this thesis showed that there was a significant 
transient time associated with acceleration from rest. In the higher speed ranges, forces 
required 4-6 seconds to achieve steady state after the steady state speed was achieved. 
The LMSC tests incorporated more speeds in a single run with fewer accelerations 
from rest. LMSC allowed for significantly less data collection time at each speed 
increment (less than 10 seconds compared with greater than 20 seconds). This may not 
have been sufficient time to achieve steady state and may contribute to the differences 
in LMSC's test results and those conducted as part of this thesis. 

The numerical analysis conducted in Chapter 4 highlighted the extreme 
complexity of the flow around SLICE. Relationships presumed negligible such as the 
transverse and diagonal body interactions were shown to contribute noticibly to the 
overall resistance value. This has an implication in the LMSC test results. When the 
August control surface tests were conducted, only one side of the model (starboard) 
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was adjusted for various angles of attack. Based on the results of the wave making 
decomposition however, drawing conclusions about whole body behavior based on 
partial both' measurements may lead to incorrect results. 

HULLFORM OPTIMIZATION 

The data decomposition presented at the end of Chapter 3 demonstrates that it 
is incorrect to optimize the hull for high speed operation by focusing on frictional 
resistance. Powering through the wave making hump does not imply that wave making 
resistance becomes negligible. It is true that the coefficient of wave making resistance 
C w becomes a small number but resistance continues to increase as velocity squared. 
The design perspective should be that Froude scalable effects will account for at least 
half the total resistance and the hull should be optimized accordingly. 

SLICE APPLICABILITY TO THE US NAVY 

Figure (5-1), obtained from a set of presentation graphics used by Dr. R. 
Compton at the U.S. Naval Academy, depicts a predecessor to SLICE. The literature 
survey conducted as part of this thesis did not uncover any design work accomplished 
on this hull. 

The comparison of either M-l or M-2 with M-3 indicate that if there is an 
advantage to the SLICE hullform over conventional SWATH hullforms. it must be 
found somewhere other than in powering performance. M-3 demonstrates that it is 
possible to create a SWATH to the design constraints of SLICE and still meet or 
exceeded SLICE resistance performance through the entire speed range. 

SLICE offers a possible advantage in the ratio of Deck Area to Displacement. 
The development of M-3 showed that the 177 Lton SWATH could have up to 600 fit 2 
less deck area than the equivalent SLICE. The extent to which the crossbox is 
cantilevered fore and aft of the struts becomes the controlling factor. For low density 
payloads, such as ferry traffic, this is significant. For higher density payloads typical 
of military’ craft the advantage is diminished. 
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Figure 5-1: A SLICE Predecessor 
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The following recommendations are made for future study of SLICE 
technology: 

1) The question of lateral lift on the aft struts created by the flow direction 
remains unanswered. This question should be addressed before the ATD vessel is 
placed in service to ensure a structural issue, either from an integrity or fatigue 
standpoint, is not overlooked. 

2) This thesis presented a reasoned approach to scaling model data to full 
scale. The ATD will be the validating test of assumptions and decisions made in 
presenting the case for using traditional methods of data scaling. If conducted 
properly, the ATD will be invaluable in demonstrating the robustness of Froude 
hypothesis. During the ATD, significant effort should be invested in photographically 
documenting flow behavior both above and below the surface. Without this 
observation of physical phenomena, specific conclusions about scaling issues will not 
be possible. 

3) The numerical approach used in this thesis requires further development 
before tuning of the SLICE hullform can be fully demonstrated. Thin ship theory 
appears to be reasonable for comparative analysis purposes and this approach should 
be further developed. After further validation, SLICE component positioning can be 
optimized and then a model of the final outcome should be built and tested in a 
similar scale to M-2 so that results can be compared. This could be the subject of a 
reasonable follow on thesis. 

4) The SLICE ATD has the potential to answer many questions that have been 
raised regarding the classical naval architecture approach to resistance predictions. 

With careful planning, the ATD will be much more than a demonstration of a single 
ship’s performance. It has the potential to contribute to the Naval Architecture 
community's understanding of the robustness of current testing and scaling techniques. 
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Appendix A 



Model Testing Lessons Learned 
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The Science and Art of Model Testing 



The purpose of this appendix is to outline the thought process that was used in 
designing the experiments and the models that were used in this thesis. The tone of 
this section will be conversational and it represents a compilation of lessons learned as 
the research for this paper progressed. It is hoped that this appendix may be of use to 
other students who may choose to undertake a similar type of research effort. 

When the work for this paper began, it was decided that in support of a Naval 
Construction program, there would be a w r ealth of experience to be gained from 
designing and building a ship (albeit a small one) that actually had to work. With this 
decided. I embarked on a challenging, frustrating, rewarding, extremely time 
consuming, detail intensive process. 

The title to this appendix is not meant to be trite. There are clearly 
scientific/engineering principles that may not be violated if reasonable test results are 
expected. Geometric similitude must be maintained, tolerances must be decided on, 
load and material strength estimates must be reasonably predicted, scaling must be 
calculated properly and data reduction must be interpreted accurately. Beyond this 
however, I discovered an entire world of additional issues regarding the planning of 
the experiment actual construction of the models and successful time and resource 
management during the critical and always too short time allocated in the test facility. 
You cannot gain an appreciation for this facet of experimentation by taking a course 
or reading a book. 

RECOMMENDATION #1: Prior to doing anything on your own experiment, 
observe similar projects for several days. In the case of this research, I had the 
opportunity to observe the SLICE testing conducted at NSWC in August and 
December. For eight days that I was able to watch other people devise solutions and 
cope with frustrations without impacting my own work. Many ideas for my own 
project were bom during these eight days and most importantly, I met the experts in 
the field of model testing. Later, these associations were invaluable to me. 
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If this period of observation does not scare you away from proceeding with 
your own experiment the next question you must honestly answer is: How handy are 
you and what resources do you have at your disposal? When I began building the 
models. I considered myself reasonably proficient with tools but I had never 
undertaken a project of this magnitude before. 



RECOMMENDATION #2: As a means of determining your handiness, I 
suggest you consider the type tolerances you will have to build to if the model is to be 
considered a reasonable geosim. For argument lets say that in full scale, the greatest 
building error that will be tolerated is 6 inches (pretty big error!). If the scale that you 
are building to is 16:1. inaccuracies in your model of greater than 3 8 inch exceed the 
building tolerance and you have introduced a level of uncertainty as to the actual 
geometric similarity of your model to the real ship. If your personality is such that you 
fail to recognize the difference between 3 8 and 1 2 inch when measuring things, then 
model building may not be the thing for you. Scaling errors to full scale to determine 
reasonableness was a process that slowly evolved during my own project 



At this stage, some decisions need to be made regarding the specific test 
objectives. To a large measure, this will determine the complexity of the model to be 
built. The speeds to be modeled mast be matched with tow tank capabilities. For this 
project, full scale speeds of 40 knots were desired. Based on the maximum carriage 
speed at the MIT tow tank, this would have required building the models in a 100:1 
scale, a size so small that no useful data would be collected. The alternatives were to 
either slow the maximum speed or find another facility. The United States Naval 
Academy was literally the closest facility' with the necessary capability'. In addition to 
these concerns, there were also the questions regarding what the model should be 
capable of doing. Early on in the project it was decided that the models would have 
no control surfaces. This removed a significant level of complexity from the model 
design. It was also decided that the models should be capable of being modified for 
self propulsion at some point in the future. Additionally, since the general test scheme 
called for a test free in heave, the models had to be built with sufficient reserve 
buoyancy to support themselves, the mounting frame, the towing carriage heave post 
and force blocks and an estimated reserve for propulsion equipment. This simple 
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exercise placed many constraints on the model: 

1) The models would have to be particularly rugged to enhance the chance of 
survival during the 400+ mile trip to the test facility in the back of a van. I eliminated 
the conventional foam construction in deference to this requirement. 

2) By understanding what the models had to do. I could make a weight budget 
to build to. It turned out that the heave post and force block weighed 25 pounds 
collectively. Properly scaled I could estimate the displacement of the models (between 
75. 83 and 115 lbf for models M-l, M-2 and M-3 respectively). I then budgeted 10% 
for reserve buoyancy and allocated 15 lbf for frame weight. The remainder could be 
used in materials to actually build the models. This weight budget further drove the 
building plan. The material selected had to be strong, extremely light and simple 
enough for an amateur to work with. 



RECOMMENDATION #3: Design your models only for the capability you 
plan to test. Build only what you need. Clearly establish exactly what the models will 
be used for before planning their construction. Eliminate multiple functions where 
possible since multiple functions will most likely increase the model's level of 
complexity. In the case of this project, significant effort was expended to meet an 
established weight budget. The weight budget was based on the assumption that the 
model would be tested in a free to heave condition but this test was never performed. 
For the fixed in heave tests, the models would have been acceptable even if they were 
negatively buoyant. 



At this stage there is almost enough preliminary groundwork set to begin 
building. With any luck you will be able to enlist the help of others in the actual 
building process. The models built for this project represent the collective work of 9 
individuals with a total of roughly 700 manhours invested between December and 
March. Admittedly a portion of this time was expended inefficiently, devising 
solutions to emergent problems, obtaining supplies etc. but these are the facts of model 
building and this sort of time expenditure should be planned for. 



RECOMMENDATION #4 If you are not literate in a CAD system, develop 
that proficiency well in advance of the model building. As early as possible, make full 



89 



scale working drawings. This is critical if other people are helping you. It is the only 
way to know that you helpers understand what to build. If the drawings are done soon 
enough in advance, they may help point out errors in the design. I did not develop a 
good set of working drawings early enough in the process and it hindered progress on 
several occasions. Once the drawing process was brought under control, I found it 
hand}’ to cut up the rather disposable drawings for templates. 

Originally PYC was the material of choice for the model construction. It is 
easy to shape and bond and is a relatively inexpensive alternative. Additionally, water 
absorption is not a concern. There are two drawbacks though. PVC is relatively dense. 
At 81 lb/in 3 , it can be made to float if rather thin walled sections are used. This is 
useful for long stretches of parallel midbody but becomes a significant drawback for 
the shaped pieces which must be turned on a lathe. Hollowing these sections out to an 
appropriate wall thickness was not practical. Additionally, the time to turn a single 
nose cone (inside and outside) was fund to be about 10 hours after considerable 
practice. A total of 16 nose and tail cones were required for two models. 160 hours 
was considered excessive to finish this single piece of the project. The second 
drawback of PVC is also one of its advantages. PVC is relatively easy to cut because 
it is soft. Although there w as little danger of it deforming in a gross sense, it is very 
easy to gouge and scratch. Once damaged in this fashion, it is difficult to repair 
because although it is easy to sand and it bonds well to itself with the proper adhesive, 
other bonding and filling agents do not adhere well. In short the probability of 
achieving a satisfactory repair once damaged is small. In light of the trip the models 
would take to the test facility and the potential for damage during this trip, this was 
considered a major drawback. 



RECOMMENDATION #5 Treat your model building project like an 
actual shipbuilding project in the preliminary phase. Account for weight and buoyancy 
to the greatest accuracy possible. I developed a weight budget for Structure (the 
frame), Propulsion (the tow pad), the Payload (the heave post), a building margin of 
10% and a design reserve buoyancy of 10% (allocated because there is some 
intention to eventually convert my models to self propulsion so there may be some 
growth). By calculating the displacement I then knew how much each component 
could weigh and I did weigh them frequently through the building process. As work 
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proceeded, building margin was allowed to be utilized at a rate commensurate with the 
stage of building. I incorporated this method early enough in the project to abandon 
PVC as a building material before too much work was done, but not soon enough to 
avoid buying a quantity of material that was not useful. 

RECOMMENDATION #6 Make frequent tests of methods you think you 
will use. well in advance of when you need them in the building project. If you plan 
to use a particular machine, check it out to ensure all the pieces are there and it 
functions properly. If you plan to use a particular adhesive to bond two different 
materials, try it out early on with two scrap pieces, even if the adhesive claims to be 
formulated to do the job. I had significant difficulty with this on several occasions. 
Testing gave me the time to either call the manufacturer for assistance or to find a 
different product. The two times I did not perform tests on intended procedures 
resulted in near disaster. The one time I did not perform a bonding test was with PVC 
and PVC cement items I had worked with in the past and was assured by the PVC 
supplier would work. For a reason I never determined, the cement did not bond. This 
occurred late on a Saturday immediately proceeding the Sunday I was to depart for 
Annapolis. The search for a new bonding agent was frustrating since it was not 
supposed to happen, difficult because of the hour, and nearly cost me a significant 
piece of the test program (the high speed SWATH model did not require significant 
reserve buoyancy and so was built from PVC to utilize the material purchased.). I 
intended to paint the models bright yellow with marine paint. Yellow is the color of 
choice for model testing because it enhances contrast and is useful for photography. 
The marine paint was difficult to apply, did not cover well and was slow to dry'. This 
provided for extreme frustration at the already frantic end of the building phase. In the 
end analysis, any spray enamel with a 10-15 minute dry time would have performed 
equally w ell in the tow tank tests, and I would have known that if I tested the paint 
systems in February rather than waiting until the end of the project. 



For the SLICE hullforms. molded fiberglass was found to be the best 
construction method. With other materials such as shaped foam. PVC rod turned on a 
lathe or stacked wood, construction is a 100% hands on process. Only one piece can 
be manufactured at a time. Each SLICE model has the following components: 4 struts, 
4 nose cones. 4 tail cones and 4 parallel midbodv sections. There were two SLICE 
models, therefore requiring the manufacture of 32 components. Hand crafting each 
component would not only have been prohibitively time consuming but also would 
have introduced slight variations between components that should have been identical 
(the four nose cones for example). Molded fiberglass construction offered the 
following advantages: 
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1. A mold insures that the four components have the same shape. 

2. Once the molds are made, component fabrication time is primarily 
epoxy curing time. This frees up time for the builder to do something else. 

3. By laying in the fiberglass cloth in the correct fashion, excellent 
strength is obtained for very little weight. 

4. Using the proper additives in the epoxy provides an outside surface 
that requires no finish work once the piece is removed from the mold. 

The West Epoxy System was the material of choice for this project. It is 
relatively expensive but is readily available at Marine Supply stores, very easy to use, 
and is associated with an excellent consumer product hotline that is responsive and 
helpful. Most of the products referred to in the building process are West System 
names. The West system can be used with either a fast setting hardener or a regular 
hardener. The fast hardener cuts the working time and also the curing time of the 
mixture in half. Additives which are called for in some of the steps also reduce 
working time. Times cited below are based on using the fast hardener. The West 
System sells a set of proportioning pumps which helps to ensure that the correct mix 
of resin and hardener is achieved. The pumps are a worthwhile investment. 

Note: The epoxy will not adhere well to plastic so working on a plastic drop 
cloth will help cleanup. Mixing the epoxy in plastic pails can be reused. Just let the 
epoxy harden and the tap the pail. 

The following outlines the steps used constructing the pods: 

1. A CAD drawing was made of the pod. Five or six tangent lines were drawn 
along the curvature. .Angles with the vertical and insets from the leading edge were 
calculated. This is shown in figure A-l. 

2. A male plug was cut out of PVC rod on a lathe. This was done by marking 
the insets, and then cutting at the appropriate angle from the inset mark to the leading 
edge. The transitions between angles was fared with a file and the entire plug polished 
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with consecutive applications of 200, 320 and 400 grit wet sandpaper followed by an 
abrasive cleanser and finally brass polish. All these steps were performed on the lathe 
to ensure roundness. 

3. A female mold was then made from the male plug. This was done with the 
West System. First the Plug is coated with a layer of wax. This is necessary to ensure 
the epoxy does not adhere to the plug. The first epoxy coat consists of the appropriate 
mix resin and hardener with sufficient microfibers added to give the mixture the 
consistency of paste. Since the additives decrease the working time. I suggest you 
pump the resin into the mixing pail, add the additive and mix it thoroughly and then 
as a last step, add the correct amount of hardener. This method will save a few 
minutes of working life. This mixture is painted onto the plug and serves as very 
smooth gel-coat layer. Brushes can be conserved by washing them with acetone after 
use but I recommend always doing this first layer with a new brush. This application 
is allowed to begin to harden but not cure (approx 2-3 hours). Even with steep angles, 
this mixture should be thick enough so there is very little running of the coat while it 
sets. A second layer is added which is also a paste like consistency. This time, the 
additive is a 50°o-50% mix of microfibers and colloidal silica, which will be less 
smooth but has greater strength. This also is set aside to begin to harden but not cure 
(2-3 hours). Finally, the real strength layer is added. For the molds, this consisted of 
squares of fiber matting soaked in mixed epoxy (no additives in the mixture this time). 
The matting has a random fiber lay and tends to be fairly (0.125" thick). Do not let 
the matting get too wet or else it will tend to slide when placed on the plug. The most 
effective way to do this step is with your hands. Latex gloves are an absolute 
requirement. This step is just like working with paper mache. When the matting has 
begun to harden (2-3 hours) a layer of 12 ounce woven fiberglass cloth is added as a 
final step. Pre-cutting the cloth into smaller patches will help since your gloved hands 
will be covered in epoxy. Since there are no additives in the epoxy in the last two 
steps, use the epoxy sparingly, just enough to barely wet the cloth. .Allow the entire 
mold to cure for the appropriate time ( 8 hours with fast hardener). The mold may 
have to be beaten pretty hard to pull it free. In the case of my shapes. I found that 
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removing the plug actually drew a significant vacuum because of the closeness of the 
fit. The mold is strong enough to be beaten pretty hard with a rubber mallet. A small 
hole drilled into the mold and later patched may also help. 

4. Manufacturing pieces is now essentially a repeat of step 3 with a few 
exceptions. I chose to leave out the fiberglass matting in the final pieces. There was 
sufficient strength with the two initial coats and the 12 ounce cloth and I wanted to 
reduce weight. Take Note: The epoxy cures using an exothermic process. In the 
steps used up till now. this was insignificant. Since the female mold is a cup, as you 
manufacture the component pieces, any excess resin will pool in the bottom of the 
mold. On two occasions I cracked my molds because of the heat released by the 
excess pool. I recommend that you set the mold into a pot of cold water at least 
during the final step (glass cloth and epoxy). This seems to be adequate to control the 
temperature. Take .Vote: You need to provide a method of pulling the piece out of 
the mold. I did this by setting two 12" pieces of nylon string into the wet epoxy of the 
final coat and covering the end with a square of fiberglass cloth. This gave a secure 
method of pulling the piece out. 

5. .After the components are pulled from the mold, they should be washed 
with acetone to remove any residual wax. If this is not done later epoxying during 
assembly may not bond well. 

6. When all the components are manufactured, they can easily be assembled by 
epoxying them together. I cut rings out of balsa wood that w ere fitted into the joints to 
increase the gluing surface. The timing of this assembly is dependent on how the 
struts get attached. For the SLICE models, the struts are bolted to the pods. In order to 
facilitate this, the tails were joined to the parallel midbody sections. The struts were 
then attached to the sub assembly and finally the nose cones sealed the assembly. 

The following process was used to build the struts. 

1. A full scale CAD drawing of the strut cross section was made. This was cut 
out and glued to a piece of thin plywood. The cross section shape was cut out of the 
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wood so that a "slip over template was created. 

2. Blocks of balsa wood were cut to rough dimensions of the strut and glued 
together to create the correct strut length. 

3. Once dry. the balsa wood was sanded to shape so that the template made a 
loose fit over the balsa wood. 

4. The balsa wood was sealed using a coating of epoxy painted onto it. This 
was left to cure. 

5. The strut was sanded to remove large defects. (The balsa wood will cause a 
fair amount of bubbling in the epoxy)- It is then covered in 12 ounce fiberglass cloth. 
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Figure A-l: Example of Nose Cone Template 
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Appendix B 

Model M-l Baseline 
Table of Offsets 
and 

Design Worksheet 
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Fowara Pod 



X fp . x i n 
0 : 
0.695 

I. 391 - 
2.086 I 
2.781 
3.477 : 
4.172 | 
4.867 

I I. 82 

1 1 3.2 ! 1 . 
14.602! 
15.992! 
17.3S3! 
18 7*3 
20.164! 
21.555! 
22.945! 
24 3361 
25.72"! 



Model M-1 Table of Offsets 

in Inches Aft Pod in Inches Struts in Inches 



Y :b,> 


1'p 


am x ; 


3P/ 


^ opO 


0 




0 




0 


! 534 1 




0.695 j 




1.534 


2.05 | 


1 


1.391 




2.05 


2.364| 




2 086 




2.364 


2.568| 


j 


2.781 




2.568 


,2.692 




3 477 | 




2.692 


2.762 




4 172 | 




2.762 


2.781 


1 


4 S67 1 


2.781 


2. "SI | 


11 !25| 


2.781 


2.6631 


12.5161 


l 


2.663 


2.44l| 


j I3.906| 


1 


2.441 


2. 1 55 j 


15 297 


i 


2.155 


: 84 9 j 


16.687 


l 


1.849 


1 54 1 j 


18.078 




1.541 


^ 1 

’ i J 




19469 




1 233 


0 925 j 


20.859 




0.925 


^617 




22.25 




0.617 


,0 308 




23.641 




0.308 


0 




25.03 1 




0 



0 1 


0.695 | 


1.391 j 


2.086 | 


ro 

■"j 

00 


3.456 1 


9.06 


9.734 


10.43 | 


1 1.125 


1 1.82 


12.516| 


13.2 1 1 1 


13 906 1 


14 6021 


I5.297| 


15.992 


16.687 



0 


0.658 


0 91 


1.037 


1.1 


1.13 


1.13 


1.105 


ii.048| 


0.978) 


|0. 892 1 


,0.788) 


0 672) 


0.542 


0.406 


0.271 


j0 136 


i 0 
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These are the full scale SLICE offsets taken from the December 1994 LMSC 
variant used to generate Model M-l offsets. 



Pod X Coord (ft) 


Pod Y Coord (ft) 


Strut X Coord (ft) 


Strut Y Coord (ft) 


1 * 


0 


0 


0 


| 1 


2.21 


• 


0.95 




2.95 


2 


1.31 


> 


3.40 


3 


1.49 


4 


3.69 


4 


1.58 


5 


3.87 


4 97 


1.625 


6 


3.97 


13.03 


1.625 


7 


4.00 


14 


1.59 


16 


4.00 


15 


1.51 


18 


3.81 


16 


1.41 


20 


3.51 


17 


1.28 


22 


3.10 


18 


1.13 


24 


2.66 


19 


0.97 


26 


2.22 


20 


0.78 


28 


1.77 


21 


0.58 


30 


1.33 


22 


0.39 


32 


0.89 


23 


0.20 


34 


0.44 


24 


0 


36 


0 
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Model M-1 Baseline LMSC SLICE 



This is based on the drawings provided by LMSC for the December baseline ship. 

la. Units Defined 

ft = I L lb = 1 M sec^lT Base Units 

run -6076-l't Nautical Mile 

, . run .. 

kt - 1 — Knot 

h r 

Iton = 2240 lb Long Ton 



1b. Foward Pod Offsets in Full Scale 

i = u j s 

\ 4'p - READ( nilplx) - READ(mlpty) 



Aft Pod Offsets in Full Scale 

\’ n p - READ(mlpox) y = READ(mlpny) 

Strut Offsets in Full Scale 



J !; 

\ s -READ(misx) y s = READ(mlsy) 

’ j ’j 



Model Scale Factor 



8 ft 

5.5625 in 



7=17 258 



100 



Model Offsets 






^ fp V. X 1 p j x fp. L 12 



Y fpv> fp; = - v fp. ;;1 ' 12 FowardPod 



X \ - \ 

a PV a P ' a P; 



. - 1 



Y ap\- ap) = ap. > ' ‘ 1 2 Aft Pod 



X sv x si :x s I 2 



Y :v = v ;. ‘ i2 

S\- 

J 



Struts 
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1c. Model Scale Weight Balance 

Fixed Variables 



/. = I" 258 



P.SU = : 9912 lb ^£1 

a 


density cf Salt Water at 56° F An assumed average 
operating temDerature. 


-> 





v ^ 3243 I () ^ ^ nematic ..scosity of Salt Water at S6° F 

sec 



■\-s s — *» slC 

I’tiv 9j6/ ib — — 

it 


density cf -'esn Water at 55 w r ~n assumed average tow 
:anK temcerature Use tms if prec.ction is for a moaei test. 


v^ v - 0804 10' 5 — 

sec 


n, nematic V scosity of Fresn Water at 68° F 


g = 3: tfsec 2 


Acceleration due to gravity 


* V = 0 kt. 1 kt.. 35 kt 


Modify this to reflect the velocity range of interest. 



Indicator Settings 



Water 




s\V 


“1" noicates the snip is operating in Sait Water 



"0“ indicates the snip is operating in Fresn Water 



p - if SW= 1 


’ P s\\ ’ f\v 


v -u:SW=l 


’ v SU ’ ’ f\V 
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Hull Displacement Definition 

... _ 1 70-lton This defines the the reauireded (,-j) displacement of the vessel 

A d : : — 

y *' note 1 70 Iton is usea here since that is the displacement of the n filletless‘ t hull. 

1 • 

V d "A d — V d = 1-lV The requireded ((j) displaced volume 

p g 

Bouyant Force 

F b = V j p g F b = 7 4 .084 * lb 

Weight Account 
FRAME 

W .'-IS lb 

Force Block - Heave Post 
w bl —5 lb 



Other Model Supplies 



W o.h -- 21b 



Desired Reserve Bouyancy 



brd --° 06 F 
brd =4 445 ' 



Allowed Weight of a single strut pod 



W 

.sp 



F b - ’.V)-- W b| - •V olll - "brd 
4 



W =6 16* lb 
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Appendix C 

Model M-2 Modified Baseline 
Table of Offsets 
and 

Design Worksheet 
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Model M-2 Table of Offsets 



Pods (both fwd and aft have the same offsets) 
in inches 



Struts 
in inches 



\mi \). in 1 ym(y) in 



0 I 


0 8 


1.1451 


1.41 |j 


1.632 i 


1.82 | 


1.9831 


2.1241 


2.2461 


2.351. 


2.441! 


2.5 1 Si 


2.582 


2.635 


2.678 


2.713 


2 74 1 


2.7581 


2.77 1 i 


2 . 7 7 9 i 


2.7811 


2 78 


2.771 


2.763 


2.75 


2.735 



0 


0 556 


1.1 12 


1 669 


2 225 


2.781 | 


j. j j7 


3.S94 


7 45 : 


\£> 

o 

o 

'f \ 


5.562 1 


6.1 19 j 


6.675 


7 231 1 


7.788 | 


8.344 j 


8 9 


9 456 


10012 


1 0 569 


11 1251 


1 1 68 1 


j 12 23X| 


j i: 7941 


13.35 1 


13.906 



\m(\). in 



14.463 


15.019 


15.575 


16.131 


16687 


17.244 


17 8 


18.356 


! 18.9 121 


19 469! 


120.025 1 


20.5811 


21.137 


2 1 .694 


22.25 


22.806 


23.3621 


23.919 


24 475 


25.03 1 j 


25.588 j 


26 144j 


26.7 1 


27 256| 


j 27 . 8 1 2 



ynuy) in 



2.716 


2.6941 


2.6691 


2.6421 


2.61 1 


2.575 j 


2.53 5 i 


2.4911 


'2 4- i 1 


:.3S5i 


(2.321! 


'2.2511 


2.1T2| 


2.0831 


1.9S4| 


l.873| 


1-5 j 


1 614 j 


1.4631 


1 296f 


1 1131 


0.9091 


0 6841 


|0 4291 


° ! 



\scale(\l ) - 1 



0 


0.72 


1.44 


2.16 


2.87 


3.59 


4.31 


5 03 1 


5.-5 1 


6.4“ 


7.19 


7 91 


862 


9.34 


1006 


1 0.78 


1 1 5 


12.22 


12.94 


13 65 


14 37 j 


15 09| 


15 811 


, 

16.531 


1 7.25 j 



yscale( \ 1 ) 
ft 

0 

0.287 
0.529 
0.729 
0.889 
1.009 
1.093 
1.144 
I 166 
1.168 
1.168 
I 168 
1.168 
1.168 
1.168 
1.168 
1.166 
1.144 
1.093 
1.009 
0.889 
0.729 
0.529 1 
0.287 

0 
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Model M-2 The Better Slice 



This variant uses the same gross dimmensions as the baseline SLICE but varies the 
body shapes to follow more traditional forms. The struts use parabolic leading and 
trailing edges and the pods are based on series 58 shapes. 

la. Units Defined 

fts|L lbsIM secs IT Base Units 

nm =6076 ft Nautical Mile 

kt = 1 — Knot 

hr 

Iron = 2240-lb Long Ton 

lb. Fixed Variables 



Psw r 1 99l21b-~- 


Density of Salt Water at 56° F. An assumed average 
operating temperature 


A 2 

v sw = 1 3343 10' 5 — 
sec 


Kinematic Viscosrty of Salt Water at 56° F 


2 

Pft v = 1.9367 lb — 
ft 4 


Density of Fresh Water at 68° F. An assumed average tow 
tank temperature Use this if prediction is for a model test 


ft 2 

Vfw = I 0804 I0' 5 — 
sec 


Kinematic Viscosrty of Fresh Water at 68° F. 


g = 32. 174*ft*sec 2 


Acceleration due to gravity 


V = 0 kt. 1 kt.. 35 kt 


Modify this to reflect the velocity range of interest. 
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1c. Indicator Settings 

Water 

$$$ SW =1 *1* indicates the ship rs operating in Salt Water 

"0" indicates the snip is operating in Fresh Water 

P =if(sw=i, Psw , Pfvv ) 

V =if(SW=l.v sw , VfVv ) 

Body Segmentation 

$$$ seg = I - ft 

2a. Hull Component and Required Displacement Definition 

*** = ! 77 Iton This defines the the requireded (<j) displacement of the vessel 

V d = Vj = 6189‘ft 3 The requireded (<j) displaced volume 

pg 

Strut Height 

H =6 ft Ammount of Strut in the Water Defined per the Baseline 



Strut Beam 



*** D | 53 25 It Maximum Strut Beam Required Defind per the Baseline 
Considered a required parameter. 

Overall Strut Length 

... L | 524 ft Overall Strut Length 



Length to 



Beam Ratio 



LDI =- 



I 



D 



LDI =7.385 Length to Maximum Diameter Ratio 
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Length Breakdown 

*** Lp m kj =7 ft Length of Parallel Midbody 

*** L n( =0.5 (L | - L pmb( i L n , = 8.5- ft Strut Nose Length (2.4 x d , s optimum) 

L , | = L | - L pmb | - L n | L t | = 8 5 * ft Strut Tail Length (3.6xDis optimum) 

Nose and Tail Shape Exponent 

*** n n , = 2.25 *** n t ,=2.25 

Calculation of Nose and Midbody Shape 

xl Oft. -beg..- L n |-Lp mb | defines the x values of the subdivisions 



D 1 , 

vnl(xl) 1 I - 

i l 



Lnl- 



'nl 



1 nl 



defines the y values of the nose based 
on a parabola 



D l 

vpmbl(xl) : — definesthey values of the midsection 

2 

Alfxl ) = iff xl i -L n j ,yni(xl ),ypmbl(xl ) 



Calculation of Tail Shape 

xl = 0 ft, - seg.. - L n j - Lp m0 | - L t |^ defines the x values of the subdivisions 



D 1 I 

vtl(xl) =— • I 

2 I 



- L 



nl 



I I" 11 ' 

^pmbl 



defines the y values based 
on a parabola 



body I ( x 1 ) = iff | x 1 1 — L n | -+- L p m b l > A 1 (x I ),ytl(xl)) Function that defines the 

body geometry 
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4 



Strut Section 




- 2 - 



- 4 1 

*24 “20 “16 “12 - 8-4 0 

— Body 1 
Body 1 



Wetted Surface and Volume Calculation 



Perimeter 

Nose 



N p =05 w 'D 1 2 .16L nl 2 ---Lln| H ' I ' nl 



D * 2 41 , - ,'D , 2 - 16 L 21 



'nl 



8-L 



nl 



D 



Tail 



N p = 17.55-ft 



1 2 2 D 1 

T p -0.5 -^D , 2 - 1 6 L t , 2 — Ini 



41, 



D | 2 - 16 L ( , 2 



8-L 



tl 



D 



T p = 17.55- ft 

Midsection 

^ p - ^ T p m ^, | 

.S p = 14-ft 
Total Perimeter 

P=N P-V T P 

P = 49. ] * ft 

Wetted Surface of the Strut 

ws s = P-H WS s = 294.602-ft 2 for one strut 
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Areas 



nose area 



N a =f >-nl D t 
tail area 


N a = 18.417* ft 2 




T a =f L.i-D , 
mid section area 


T a = 18.417-ft 2 




S a =D l L pmbl 
Total Area 


S a = 22.75*ft 2 




Area r N a + S a -T a 
Volume per Strut 


Area = 59.583-ft 2 




vol = AreaH 
Total Strut Volume 


vol = 357.5-ft 3 




VOL =4 vol 
Pod Volume Required 


VOL = 1430* ft 3 




V pods ^ V d - V 0L 


V pods= 47 ^ 


for four pods 


\, . V pods 

V P°d 4 


V p od = ll90*rt 3 


per pod 



Pod Definition based on the Model 41 55 of Series 58 

). | = 4.3 the scale factor used to get from the data on page 49 of senes 58 (Lpod= 9 

ft) to the size needed per baseline Lpod= 38.7 ft) is 4.3 

L p = 9 ft X , D p = 1.8 ft J. | D p = 7.74*ft 

i =0.1. .53 

x. = READ(betx) y. = READ(betv) non dimmensional 

xfs(x) =x-L p yfs(y) =yD p 

WS 58 = 39.75 ft 2 
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full scale dimmensions 



for one pod 



WS fs '.VS 58 , 2 WS fs = 734 977-ft : 

Vol < s - 14 89 fr* As printed in the Series 5S data 

Vol = Vol 58 /. , 3 Vol fs = 1 1 83 .859* ft 3 

VOLUME - 4 Vol fs - 4 vol 

VOLUNtE = 6 1 65 .437* t7 J V J = 6188.727'tV Compares calculatea volume 

to needed volume. 

A - VOLUME p g A = 176.334 -lion 



Mow Scale to Moael Size 
5 5625 in 




l -0.1 25 

\mi\) : \ i s( \ ) 



7 ^ = 0.06 

j = 26. 27 50 
n =>1s(v )•?.-> 



— = 16.698 
7 o 



Better Slice Strut Offsets }.*> = o 06 



\ i i 2 A n 

xscaiei \ 1 ) - “ 

11 



yscale( \ 1 ) - body 1 ( \ 1 )• 1 2 7 
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Weight / Volume Balance 



Density of PVC 



D r -6 


1875 in 


L r = 7.9375 in 


diameter and length of i 


D r 2 




weighed 


V r = * 


T L '- 


V r = 0 1 38-rt' 


volume of 1 foot of rod 


P pvc 


1 ! .7 lb 
V r 


P pvc = H4 70S* :b-3l •' 


density of the PVC 


■’ !\v - = 


h2 Hb-i 


1 


density of FW at 68°F. 



Displaced Volume 

V mode! -'VOLUME;.,' 



^ model 1 



Bouyant Force 

1 - V o • •’ 

' o mooe; 1 1\\ ^ 

Weight Account 



FRAME 



Fu =s:. 523-lb 



W r = 15 lb 

Force Block Heave Post 
W hl ,25 1b 
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Parallel Midbody X 4 pods 



w 5.5625 2 5.5625 - 0.5 

’ pvcpmb ~ ' “ in ~ r ‘ “ in: L pmb 4 

V pvcpmb “0.065*n 

pmo ~ ^ pvcpmb' ? pvc ^ pmb “ ~ 



Pod Nose ana Tail x 4 poas 

pour.t ~ ^ pods '■ 2 ■ 4 ^ ; 

V ;mu:n =0645 ^ 

v :muicl t (2.25 ini 2 l 25 in 2 
’ n bore 2 “^(2.0-inj 1 m 2 
^ :ibore3 " 1 ^ ^ ^ in ' 2 in - 
'2 «hore4 r 1 ( 1 25 in )“ 1 .0 m 2 

'* :;bore5 "" 1(075 ^ 1 0 *« 2 
r tbore I - * ( 2 25 m ) 2 in 2 

V tbore2 ~ 1 ( 2 0- in j 2 3 in 2 

V iborej =n i 1.5-inj 2 2 in 2 
v :bnre4 : * ( 1 () ' m f 2 m 2 
'2 ihopjf “(05 in r ! m 2 



— — — m volume of the nose 

*) 

and tail 

^ nborei ” ^ * in 
^ nbore2 ~~ ^ 55 Mn 

V nbore3= 28 ^ 4 ^ 
V nbore4= 9817 -' 11 '' 

V nborei = 7 534 ,,n ' 

V ibore I = 63 617 -in' 

^ lbore2 - 73 398 " in 

V tbore3 =28274 -> n? 

V thore4 = l2 566 ’ in '' 

37 iborei “ ' 4 ' * " in 



11 ? 



V 



bore 



V nborel ~ ^ nbore2 V nbore3 ^ nbore4 v nbore5 ^ V tborel 
~ ^ :bore2 ~ ^ tbore3 “ ^ tbore4 ^ ^ tbore5 



^ bore 



= 287.947 *in : ' 



^ pvcpod " ^ podnt ^ bore 
V pvcpod= 0479 * ft3 

^ poant ~ ^ pvcpod'P pvc ^ podnt ~ 4 ^ ’ ^ 



4 Struts 



NOlM'C 



>boi\ 



V 



pvcstr 



W 



en 4}.-," 6 m 


V . = 0 42“'- fr 

s 


- ( 10 in 1 S in 6 in 


V Shore 


i 4 ( 2.75 in 0.7 in o m)-2 


V sbore2 =005 > lV '' 


: V - V , - V u 

s shore more2 


V pvcstr “ ] 3 ~ 4 ’^ 


“V o 

pvcstr 1 pvc 





W = 10 499Mb 
struts 1 11 



Other Model Supplies 



Checking for a Balance 

F brd - 10 1b Desired Reserve Bouyancy 

w model t w bl ^ pmo podnl struts oth 
W model = ]03 - 527 * lb A ^ual Model Weight 

F brn - F h - W mode j Actual Reserve Bouyancy (Should = R bd ) 
F br a =‘19 005 Mb 

1 ; bra F bm = -29 005 Mb Must be 0 or greater. 
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Areas to Help the Balance 



reserve = areas with excess weight included that can be tnmmed down by the 
indicated ammount 

r f - 2 lb frame = 0 lb force block r olh = 1 lb other mode! items 

reserve =r f +r w -r olh 

fixed = W f-r W j^j - W ol ^ - reseive Weights that are unchangable 

excess = reserve bouyancy that is desired but not needed 
excess - ' > lb 

* hr * u rd “ 

Calculate the upper limits for strut'nose and tail cone weights 
u modelmax " ^ b “ r br 



x . r _ u modelmax ~ Jixe ^ 


The Max Allowed weight for a 


v ' s P a i 

4 


single strut/pod combination 


W spa = 7 63 1 "lb 




^ stmts ^ podnt ^ pmb 


= 14 1 32* lb Calculated Weight of a single 

strut/pod combination 


4 4 4 


Summation of Allowed Weights 




'.V f - i ( = 1 3-lb 


Frame Weight 


W b |-. b |=25-lb 


Force Block and Heave Post 


£ 

-r 

II 

i 

c 

2: 


Miscalleanous Weight Allowed 


W S pa= 7 63Nb 


The Max Allowed weight for a 
single strut/pod combination 



Model Dimmensions 



" m 16.68 




Lstrutjn 


Lstrut m = 1 439- ft 


WSstrut m = WS s ;. m 2 


WSstrut m = 1.059* ft 2 


Lp°d m ' L p /• m 


Lpod m = 2.32’lt 


WSpod m 5 ws & ;. m 2 


WSpod ;11 = 2. 64 2 -ft 2 
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Appendix D 

Model M-3 Optimized SWATH 
Table of Offsets 
and 

Design Worksheet 
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Model M-3 Table of Offsets 



Demihull Offsets 
in inches 



Strut Offsets 
in inches 



xh!2 ft' 1 • ;20FF(xhl2; :V 



0 




0 


0.794 




0.667 


1.589 




1 212 


2.383 




1.641 


3.178 




1.959 


3.972 


! 2 173 


4 766 




[2.294 


5.561 




2.335 


6.355 




|2.303 


7 15 


12.285 


7 944 


j 2.26X| 


8.738 


1 2.25 1 


9 533 


2.232 1 


10 327 


2.2141 


1 1.1221 




i2.196| 


1 1.9161 


| '2.178 


12.71 


2.16 


13 505 


12.142 


14 299 


,2.124 


15 094 


2.107 


15 888 


2.0X9 


16.682 


1 


2.071 


17.477 


2.053 


18.271 


2.035 


19.066 


2.017 


19 86 


■2.313 


20.654 


1 


2.48 


21.449 




2.647 


22.243 


i 


2.814 


23.038 




2.78 


23.832 




2.78 


24.626 




2.78 


25.421 




2.78 


26.215 




2.78 


27.01 


i 


2.78 


27.804 


i 


2.78 


28.598 


r 


2.78 


29.393 


2 78 



\l ft' 1 - 12 



0 


|0.794 


j 1.589 


|2.3X3 


,3.178 


13.972 

j 


j4 766 


‘5.561 


\6 355 


“ 15 


j7.944 


fs 738 


,9.533 


10.327 


1 1.122 


1 1.916 


12.71 


13.505 


14 299 


15.094 


15.888 


16.682 


17.477 


18 271 


19.066 



body 1 ( n 1 ) ft' 1 12 



0 


0.26 


0.488 


0.685 


0.852 


0.99 


1.099 


1 182 


1 24 


1 275 


1 29 


1 291 


1 291 


1.291 


1 29 


1 275 


1 24 


1 182 


1.099 


0.99 


0.852 


0 685 


0 488 


0 26 


- 1.14710' 15 



Model M-3 The High Speed Swath 



This variant uses the same gross dimmensions as the baseline SLICE but uses a 
shape characteristic of a high speed swath. Tandem strut design was selected to 
minimize strut displacement. 

la. Units Defined 

Its 1 L lbsIM secs IT Base Units 

nm - 6076 ft Nautical Mile 

kt = 1 — Knot 

hr 

Hon - 2240 lb Long Ton 

lb. Fixed Variables 

-0 0662 



Psu - : 9912 !b ^ 
ft 


Density cf Salt Water at 56° F An assumed average 
operating temperature 


. , , . , . 5 It 2 

v s\\ ■ • 10 — 

sec 


Kinematic Viscosity of Salt Water at 55° F 


= : 9367 lb — 

it 4 


Density of Fresh Water at 65° F An assumed average tow 
tank temperature. Use this if prediction is for a model test. 


- ! 0804 10' 5 — 
sec 


Kinematic Viscosity of Fresn Water at 68° F 


g = 32 174* 11* sec 2 


Acceleration due to gravity 


V = O ku 1 kt .. 35 kt 


Modify this to reflect the velocity range of interest. 



1c. Indicator Settings 



Water 

$$$ SW -U V" indicates the snip is operating in Salt Water 

"0" indicates the snic is oDerating in Fresh Water 

p -if, SW=l .p sw ,Pj\ v . 
v :ifSW S |.._ sw .v fu , 

Body Segmentation 

SSS 'cc = ! ft /. 

2a. Hull Component and Required Displacement Definition 

i \ =177 lion /.' This defines the the reauireded (^) displacement of the vessel 

V j =A c j*— V j = 2-iV > The requireded ( d ) displaced volume 

P £ 

2a. i. Strut Creation 
Strut Height 

H = x ft /. Desirea to allow demi hull suomergence of at least one diameter. 

1 1 =6.355 *in 

Strut Eeam 



Bj 3 25 ft/. Maximum Strut Beam Required Defind per the Baseline 
Considered a reauireo parameter. 

B j = 2.5X2 mii 

Overall Strut Length 

... L | = 24 ft /. Overall Strut Length 

L j = 19.066 *in 

Length to Beam Ratio 
L | 

LB I = LB1 =7.3S5 Length to Maximum Beam Ratio 
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Length Breakdown 

*** Lpmbl -^TtX Length of Parallel Midbody 

L n i r 0 5 L j - L j L^j = 0. 695* 1 1 Strut Nose Length ( 2 . 4 x Dis optimum) 

L - L | — L p m ^ | - L n j L t j = 0.695* h Strut Tail Length (3.6 x dis optimum) 



Strut Nose and Tail Shape Exponent 



r ‘ ns 5 125 



n ts .2.25 



Calculation of Nose and Midbody Shaoe 
\1 - w li. L n j - L j 



defines the x values of the subdivisions 



yn 1 ( \ 1 ) 




\l 



n ns i 



L 



nl 



defines the y values of the nose based 
on a parabola 



ypmb 1 ( \ 1 ) 




defines the y values of the midsection 



Al(.\l) = if \1 _L n | , \n 1 ( x 1 ),ypmbl(\l ) 



Calculation of Tail Shape 



\1 t) ft- -sec L n j Lp nlD j L t j 



defines the x values of the subdivisions 



m I ( \ 1 ) 



^ni *'pmbl 



defines the y values based 
on a parabola 



bodv 1 c \ 1 j - if !\1 _L n j - L„ m ^ j , Al(xl ),yil(\l ) Function that defines the 

body geometry 
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Strut .Section 




-16 “ I 4 -1.2 “I -0 8 -0.6 -0.4 “0.2 0 



2.a.ii. Wetted Surface and Volume Calculation For the Struts 

Perimeter 

Nose 



N ,, i 5 B , “ - 16 L —In. 



3 f , 4 L n 1 ~ B | ‘ - 1 6 L n | ' 



s L 



nl 



N p = i4::.,t 



Tail 



i)5 B r - 16 L 



U * 



B 1 

S L, 



4 • L ♦ i — B i — I6L 



In i 



tl 



T P = ' J 3l-ti 
Midsection 

p " " w rmb 1 

S p ='> -'5"-rt 

Total Perimeter 

p : N - S - T 
p p p 

P = 3.24-li 

Wetted Surface of the Strut 

WS s - P H WS s = I 716-l't 2 for one strut 
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Areas 



nose area 







tail area 



T a V L .I B I 



mid section area 
r B , L 



pmb 1 



Total Area 

\rea N , S 

Volume per Strut 
\ui - Area H 

Total Strut Volume 

VOL s = 4 vol 



N a = 0 I -It 2 



T a = 0.1 *1\ 2 
S a = 0.04>ft 2 



\rea =0 242* !'t 2 



■. oi =D i :s-rf 



VOL s =0.513-11'' 



2.a.iii Demi-Hull Creation 



Demi Hull Volume Required 



dh =Vj-VOL s V Jh = !.f.-' 



v dh 



= o.7*fr 



for two hulls 
per demihull 



12 ? 



Demi Hull Shape 



D | -7 ft /. Max diameter of the Hull D ^ = 5 56 1 *m 



D hi 


= 0 84 D h 


D j^j =4.671 *in 


Section #1-2 interface 


D h2 


' 73 D h 


D h2 = 4 059 Mil 


Section #2-3 interface 


D h3 


= D h 


D h 3 - 5.561 -in 


Section #3-4 interface 


D h4 


= D h 


D j l4 = 5.561 mi 


Section #4-5 interface 


D hS 


■'’73 D l j 


D = 4 059*m 


Section #5-6 interface 


»>!«, 


N4 ;•) 


D h( , =4(ri* :n 


Section #6-7 interface 



=90 ft/. Overall Length of the null. L] 1 =5 95S-n 



L h 1 =,)0XL h 


L h , =5.72un 


Section #1 


b h 2 1 ^ ^ h 


L = 1 3 5S4 *m 


Section #2 


L h3 r:105L h 


L h3 =' 575-.n 


Section #3 


L h4 = ll25L h 


L i^ = 17 S74 -in 


Section #4 


^ h5 " ' } ^ b )i 


L ^3 = 10 ^24 nn 


Section #5 


L h 6 :;,09L h 


L h 6 = n 435 ’ in 


Section #6 


L h7 ' !) 19 L h 


L h 7 ^ 1 4.5S4 -in 


Section #7 


l-h 

1.20 : — L2D = 

19 h 


I2.X57 




Segment 1 & 2 






Calculation of Nose Cone Shape 




n nl " 2 - 25 




Nose cone shape exponent 


\ltl2 = O ft,- seg.. L - 


L h 2 i 


defines the x values of 
the subdivisions 
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the subdivisions 



D 



vh 1 ( \h 1 2) =■ 



L j - \h 1 2! 



n nl 



defines the y values of the 
nose based on an ellipse 



vh2(\h!2) 




0 h2 ~ D h I . r 
; L h2 



\jii 2 - 



defines the y values of 
seg sectjon 2 



Al(\hl2) - if' xli 1 2 i <L h] ,\hlf \h!2),yh2(\h 12) 



offsets for the first two body 
segments. 



Adding in Segment 3 



\h 1 2 - 0 ft, - seg L j | * L \ ^ - L , 3 defines the x values of the 

subdivisions 

D ip D h 3 - D 10 

\h3f\hl2; \hl2-L .. - ^ to - see : defines the y values 

- ~ :.3 of section 3 



Bl( \hl 2) - if \h!2 <L h] - L h2 ,Al(\hl2),\h3t \hl2) 



offsets for the first three 
body segments. 



Adding in Segment 4 

defines the x values of the 
subdivisions 

defines the y values of 
section 4 



\ii 12 u It. >eg I- hi ~ L !i2 " L Ii3 ' L h4, 
H hi 

\h4(\hl2) - 



Cl(\hl2) = if \h!2 <L j 1 j - L j.p - L ,B l(\h; 2'. , yh4( \h 12) 

offsets for the first three body 
segments. 

Adding in Segment 5 

\h 1 2 - 0 fu- seg I- hi “ L jp - L i n 3 * L ^ - L 1 ,^ defines the x values 

of the subdivisions 

^ h4 ^ li5 " ^ h4 - 

yh5(,\hl2) - — xh!2 ^ L hl - L h2 - L h3 * L h4 - seg 

- L h5 

defines the y values of 
section 5 



Dl(.\hm =if \hl2 <L h) -L h2 -L h3 -L h4 .CK>Jii2).>h5(\lil2)j 



offsets for the first three body 
segments. 



Adding in Segment 6 



\hl2 -Oft.->eg.. L - L h 2 - L ^3 - - Lj^ - L defines the x values of the 

subdivisions 



D 



vh6( \h!2'i 



h5 



n Ii6 ^ hi - , . . - j 

— ' JlI2 - L hl - L h2~ t '!i3~ L h4- rL h5/- se S 

- L lid • : 



defines the y values of 
section 6 



E 1 ( \h 1 2 1 



\ii i — <L 1^ - L ip - L i,- - L | i: j — L i,j,Di( \ii i 2 ) , vh6( \h 1 2 ) 

offsets for the first three body 
segments. 



Adding in Segment 7 



\h!2 - 0 It. - seg L h i - L h2 -L h3 - L h4 - L, )5 - L h6 - L h7 



defines the x values 
of the subdivisions 



\h7(\hi:, 




D 



h6 



2 L, 



\hl 2 



111 - L Ii2 ~ L h3 ~ L Ii4 ~ L hi L h6 " I 



defines the y values of 
section 7 



OFF(\hir, : :! \h 1 2 : <L hl - L h2 - L h3 - L h4 - L h5 * L h 6 .EI(\hl2),yh7(\hl2)) 

offsets for the first three 
body segments. 



0.4 

o 2 

OFF( \lil2) 

0 

-0 2 

“6 "5 -4 -3 -2 "I 0 

\hl2 
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Wetted Surface of TWO demihulls 



•-L h 

WS h -2. 2 n ()FF ( xJi 1 2 ) d\hl 2. WS h =13.2-ft 2 

. 0- If 



Volume of TWO demihulls 

*-L h 

VOL h =2 aOFF(xhl2) a dxhl2; VOL h = l.3-ft 3 

-0 ft 



Volume Balance 

Available 

V() Lava,i = VOL h - v OL a 



VOL avail = I.S-lV' 



' avail v ^ ^ avail P £ 



'avail 1 



Required 



V d =lSMV 



Excess Volume 



avail ’ d U * !I 



(negative implies insufficient volume) 



Deck Area Lost compared to SLICE 

C - 3 • It This is the combined ammount the bow/stern are 
cantelivered over the ends of the hulls. 

- h 

105 ti C : 55 ti = 660* tr (negative implies a gain in deck area) 

t. 
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Model Building Balance 



Bouyant Force 



F b avail 

Weight Account 

FRAME 

W f = 18 lb 

Force Block + Heave Post 

w bl - 5Ib 

Other Model Supplies 
W o.h --- lb 

Desired Reserve Bouyancy 


F b = 1 12.496-lb 


F brd = 005 F b 
F b,- d = 5 625-lb 





Allowed Weight of a aeminull + 2 Strut combo 



u, . F b- W f- W b|- 

sp 2 


^ olh r brd 


W =30 936- lb 





1 



Appendix E 



Project Test Plan 
and 

Test Journal 
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SLICE Project Test Plan for the Week of 6-10 March 1995 



Anticipated Number of Models to Test: 4 

Desired Data to Obtain: Resistance (x) for all 4 models. Dynamic Sinkage (z) for 2 

models. 

Model Descriptions: Model M-l is the LMSC baseline design. 

Model M-2 anticipates improvements to the LMSC baseline by 
improving the shape of the bodies but retaining the gross 
geometry of the baseline ship. 

Model M-3 is an optimized SWATH based on the results 
presented by McGregor 
Model M-4 is the Mid Foil Concept. 

Desired Testing Schedule: 



6 March 95 Set-Up M-l; M-l Test 1 20runs 

7 March 95 M-l Test 2; Set-Up M-2; M-2 Test 1 30 runs 

8 March 95 M-2 Test 1 (finish) ; Set-Up M-3; M-3 Test 1 30runs 

9 March 95 M-3 Test 1 (finish); M-3 Test 2 30 runs 

10 March 95 Set-Up M-4; M-4 Test 1 20 runs 



Test Description: 

Test 1: Resistance Locked in Heave and Pitch. 

Objective: To measure forces in the X and Z directions over the 
operating range of speeds. 

Discussion: Previous testing has indicated that a significant Z force is 
generated and that the magnitude of the force is geometry dependent. This method of 
testing replicates the LMSC methodology and therefore allows for a basis of model 
comparison with LMSC results. The LMSC justification for this methodology is that 
control surfaces will be employed to fix the operational draft. 

Description: With the model rigidly fixed to the carriage, the following 
sequence of speeds will be used to collect data for both X and Z forces : Model Scale 
Equivalent (MSE) of 3,6,9,12,16,20,25,30,35,40 kts for basic data. Each speed will 
accomplished with one run. Then MSE of 20,25,30,35,40 kts to demonstrate 
repeatability. If time in the day permits then additional runs will be made at MSE of 
10,11,13,14 and 22 kts to fill in hump data and the high speed knuckle. This 
represents a total of 20 runs for the test completion. 



Test 2: Free to Heave. 



Objective: To quantify the difference in resistance results for a model 
rigidly fixed and one that is allowed to oscillate about a desired waterline. 

Discussion: Although control surfaces may hold the ship to a fixed 
depth in a gross sense of the word, it is unlikely that the system's sensitivity will 
exactly replicate a rigid mounting. This test will use a counterbalance to replicate the 
control surface but allow the model to oscillate about this point. 

Description: This test will use the Z force results from Test 1 to 
determine the counterbalance weight required. The model will be attached to the 
carriage free to heave but locked in all other respects. The necessary counterbalance 
weight be added and the model towed at the given speed. (Note, due to the small TPI 
associated with SWATH vessels, a method for adding the weight as the carriage 
conies to speed may be necessary'. This is a Discussion Point for the 3 Feb meeting.) 
The test ng will require a positive stop in the negative z direction with a method of 
determining that the stop has been reached. This stop must be set to no more than 3 
inches above the design waterline of the model. • 

The following sequence of speeds will be used to collect data for X 
force only with note being taken of the counterbalance weight added. Model Scale 
Equivalent (MSE) of 9,12.16,20.25,30,35,40 kts for basic data. Each speed will 
accomplished with one run. Then MSE of 20,25,30,35,40 kts to demonstrate 
repeatability. If time in the day permits then additional runs will be made at MSE of 
3, 6,10,1 1,13.14 and 22 kts to fill in hump data and the high speed knuckle. This 
represents a total of 20 runs for the test completion. 
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Appendix F 



Model Test Results 
Raw Data 
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Index of Test Runs Conducted 6-10 March 199S 



Run# Model/Test 


Full Scale 


Test 


Drag 


Lift 






Equivalent \, 
(Kts) 


Velocity 
(FPS) 


(Ibf) 


(Lbf) 




1 


M-1 in Baseline Config. 


10 


4.08 


6.46 


1.27 


2 


M-1 in Baseline Config. 


3 


1.21 


0.151 


0.191 




M-1 in Baseline Config. 


6 i 


2.41 


0.689 


0.562 




M-1 in Baseline Config. 


9 


3.64 


4.24 


1.09 


3 


M-1 in Baseline Config. 


12 


4.85 


9.08 


2.42 




M-1 in Baseline Config. 


13.5 


5.47 


6.16 


3.33 


61 


M-1 in Baseline Config. 


16 


6.47 


5.61 


5.18 


7 


M-1 in Baseline Config 


20 


8.08 


7.19 


7.73 


R7 


M-1 in Baseline Config 


20 


8.081 


7.17 


77 


81 


M-1 in Baseline Config. 


25 


10.11 


10.66 


10.53 


9! 


M-1 in Baseline Config. 


30 


12.13 


19.54 


15.52 


R9 


M-1 in Baseline Config. 


30 


12.13 


19.39 


15.46 


10 


M-1 in Baseline Config 


35 


14.14 


28.82 


21.02 


11 


M-1 in Baseline Config 


-0 


16.18 


35.32 


26.42 


12 


M-1 in Baseline Config 


‘0 


4.05 


621 


1.02 




M-1 in Baseline Config. 


11 


4.45 


8.91 


1.64 


14i 


M-1 in Baseline Config. 


13 


5.25 


6.86 


2.74 




M-1 in Baseline Config. 


14 


5.66 


5.65 


3.57 


16 


M-1 in Baseline Config. 


15 


6.02 


5.48 


4.25 


17 


M-1 in Baseline Config 


40. 


16.18 


35.38 


26 65 


18 i 


M-1 with Spray Rail Aft 


30 


12.13 


19.45 


15.48 


191 


M-1 with Spray Rail Aft 


35 


14.14 


27.95 


20.61 


20 


M-1 with Spray Rail Aft 


40 


16.18 


36.87 


26.87 


24i 


M-2 in Baseline Config 


10 


4.11 


6.83 


0.443 


25! 


M-2 in Baseline Config. 


3 


1.22 


0.172 


0.247 




M-2 in Baseline Config. 


6 


2.46 


0.667 


0.845 


27 


M-2 in Baseline Config. 


9 


3.68 


4.52 


0.416 




M-2 in Baseline Config. 


12 


4 93 


9.92 


2.37 


28 


M-2 in Baseline Config. 


6 


2.46 


0.679 


0.724 


29 


M-2 in Baseline Config 


12 


4.93 


9.76 


2.51 


30 


M-2 in Baseline Config 


16 


6.58 


6.28 


5.95 


31 


M-2 in Baseline Config 


16 : 


6.58 


6.3 


5.89 


321 


M-2 in Baseline Config 


20 


8.23 


7.25 


9.35 


R32 


M-2 in Baseline Config. 


20 


8.23 


7.26 


9.29 


34 i 


M-2 in Baseline Config. 


25 


10.28 


10.23 


12.97 


R34 


M-2 in Baseline Config. 


25 


10.28 


10.19 


13.03 


361 


M-2 in Baseline Config 


30 


12.35 ; 


14.54 


17.52 


R36 


M-2 in Baseline Config. 


30 


12.35 


1453 


17.48 


38 


M-2 in Baseline Config 


35 


14.4 


20.17 


22.33 


R38 


M-2 in Baseline Config. 


35 


14.4 


20.19 


22.3 


40 


M-2 in Baseline Config. 


40 


16.45 


26.42 


28.25 


R40 


M-2 in Baseline Config. 


40 


16.45 


26.36 


28.35 


42' 


M-2 in Baseline Config. 


11 


4.47 


9.81 


0.77 


43! 


M-2 in Baseline Config. 


11.5 


4.72 


10.78 


1.82 


44 


M-2 in Baseline Config 


13.5 


5.47 


6.68 


3.25 


45: 


M-2 in Baseline Config 


12.5 


5.22 


785 


2.86 


461 


M-2 in Baseline Config 


14.5 


5.97 


6.05 


4.27 
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48 


M-3 


6 : 


2.58 


0.801 


1.01 




M-3 


9! 


3.88 


4.22 


1.981 


49 


M-3 


121 


5.18 


6.08 


1.081 


53 


M-3 


16! 


6.91 


8.24 


4.481 


55 


M-3 


16 1 


6.92 


8.22 


4.63 


56 


M-3 


181 


7.78 


8.93 


6.41 


57 


M-3 


20! 


8.65 


9.74 


7.91 


58 


M-3 Towed Backwards 


7.5 : 


3.48 


3.17 


1.26 


59' 


M-3 


25 i 


10.81 


13.87' 


121 


R59 


M-3 


25: 


10.81 


13.86 


11.92 


61 


M-3 


30! 


13 


19.13 


16.68 


R61 


M-3 


30 i 


13 


19.03 


16.56 


63 


M-3 


35 i 


15.13 


26 


22.39 


R63 


M-3 


35' 


15.13 


25.75 


22.19 


66 


M-3 


40 


17.29 


33.27 


28.84 


R66 


M-3 


40' 


17.29 


33.1 


28.7 


68 


M-3 


11.5 


4.97 


6.35 


1.27 


69 


M-3 


14’ 


597 


6.74 


2.27 


70 


M-1 in Baseline Corrfig. 


10.25 


4 17 


6.01 


1.87 




M-1 in Baseline Corrfig. 


11.25 


4 58 


7.14 


1.42 




M-1 in Baseline Corrfig. 


12. 


4 78 


6.75 


1.27 


71 


M-1 in Baseline Config. 


10 


405 


6.37 


1.13 




M-1 in Baseline Config. 


11.5 


468 


9.43 


2.29 


72 


M-1 T+0.5" 


11.5 


4.71 


9.59 


-1.13 




M-1 T+1.0" 


11.5 


4.71 






74 


M-1 T+0.25” 


11.5 


467 


9.42 


0.89 


R74 


M-1 T+0.375” 


11.5 


AT 


9.45 


0.4' 


76 


M-1 Towed Backwards 


20 


8.07 


-9.91 


10.22 


77 


M-1 Draft Adjustment 


15 


5.97 


See Plots 


See Plots 




M-1 Draft Adjustment 


15 


5.97 


See Plots 


See Plots 


78 


M-1 Draft Adjustment 


15 


5.97 


See Plots 


See Plots 




M-1 Draft Adjustment 


15 


5.97 


See Plots 


See Plots 




M-1 Draft Adjustment 


15 


5.97 


See Plots 


See Plots 


79 


M-1 Draft Adjustment 


17.5 


6.97 


See Plots 


See Plots 


80 


M-1 Draft Adjustment 


17.5 


6.97 


See Plots 


See Plots 


81 


M-l Draft Adjustment 


17.5 


697 


See Plots 


See Plots 


82 


M-1 In Line Config 




4.05 


6.32 


1.27 




M-1 In Line Config 




4.85 


9.07 


2.67 


83 


M-1 In Line Config 




1.21 


0.158 


0.175 




M-1 In Line Config 




2.41 


0.693 


0.6 


84 


M-1 In Line Config 




364 


4.41 


1 ! 




M-1 In Line Config 




6.47 


5.74 


5.36 


85 


M-1 In Line Config 




8.08 


7.1 


7.54 


86 


M-1 In Line Config 




10.11 


11.48 


10.26 


87 


M-1 In Line Config 




12.13 


22.29 


14.22 


88 


M-1 In Line Config 




14.01 


35.17 


19.67 


89 


M-1 In Line Config 




16.19 


44.51 


26.45 
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THST: Model M-l Test I RIJN NC):00l 
05 / 30/95 01 : 44:43 
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DRAGOO OOI(LB)-AL.L T229/S0 
VliLOOO OOI(FPS)-ALL T229/S0 



TEST: Model M-l Test I RUN N():002 
05/30/95 01:48:09 
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TIME(SEC) 



TEST: Model M-l Test I RUN N0:003 




147 



TIME(SEC) 



TliS I Model M-l RUN N( );(><)(> 




148 



TIM H( SEC) 




TI-ST: Model M-l Test I RUN NO: 007 
05/30/95 01:59:58 
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ohsuni.i. 



ILiS I 'Model M-l lest I RUN N0.008 
05/30/95 02:03:35 




150 



TIME(SEC) 



TEST: Model M-l Test I RUN N0:009 
05/30/95 02:06:40 




151 



TIME(SEC) 



1 1 iS I’: Model M-l Tcsl I RUN NO R09 
05/30/95 02:10:44 




152 



TIME(SEC) 



TEST: Model M-l l est I RIJN NO:OI I 
05/30/95 02:17:02 



Csl 




153 



TIME(SEC) 



Tl iST: Model M-l lest I RUN NO 012 
05/30/95 02:20:27 




154 



TIME(SEC) 



TEST: Model M-l Test I RUN NO OI4 
05 / 30/95 02 : 24:06 




155 



TIME(SEC) 



TEST: Model M-l Test I RUN NO OI6 
05 / 30/95 02 : 28:26 




156 



TIME(SEC) 



TEST: Model M-l Test I RUN NO 01 7 



rs* 




157 



TIME(SEC) 



I'l'S'l : Model M-l lest I with Spray Rail RUN NO:OI8 
05 / 30/95 02 : 36:1 1 




158 



TIME(SEC) 



TEST:Model M-l Test I with Spray Rail RUN NO 019 






I 




159 



TIME(SEC) 



TEST: Model M-l 'l est I with Spray Rail RUN N0:020 
05/30/95 02:42:55 




-r c-n) rg 



160 



DRAG00 020(LB)-ALL T229/S0 
VELOOO 020(FPS)-ALLT229/S0 



11 iST: Model M -2 Test I RUN NO 025 
05/30/95 02:48:08 




161 



TIME(SEC) 



THST: Model M-2 Test 1 RUN NO:027 
05/30/95 02:51:38 




162 



TIME(SEC) 



TEST: Model M-2 Test I RUN NO 028 



^T) 

ri 

o 

o 

r^ ( 

in 

o 




163 



TIME(SEC) 



TEST: Model M-2 Test I RIJN NC):030 
05/30/95 03:27:50 







164 



TIMH(SEC) 



TF.ST: Model M-2 Test I RUN NO 03I 













165 



TIME(SEC) 



I liST: Model M -2 l est I RUN NO:032 
05/30/95 03:34:52 




166 



TIMt:(SE:C) 



TEST: Model M-2 Test I RUN NO:R32 



rn 

ri- 

ce 

rn 

rn 



O' 



WO 

O 




167 



TIME(SEC) 



TEST:Model M-2 Test I RUN N():034 
05/30/95 13:41:51 




C. yr, ~ ~ 



168 



TIME(SEC) 



THST: Model M-2 Test I RUN NO.R34 
05/30/95 13:43:50 




169 



TIMH(SEC) 




TliST: Model M -2 Test I RUN N():036 
05/30/95 13:45:37 







'r, ur. \r, 

rsj — • 
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TIME(SEC) 



TEST Model M-2 lest 1 RUN NO R36 




171 



TIME(SEC) 



TEST: Model M-2 Test I RIJN NO:038 
05/30/95 13:49:19 




172 



TIME(SEC) 



cc 

rO 

v 



04 

— 

- ^ 

04 0 s 
2 § 

“3 

-3 o 



y: 




173 



TIME(SEC) 




TEST: Model M-2 Test I RUN N0:040 
05/30/95 13:53:50 




174 



TIME(SEC) 



TKST: Model M-2 Test I RUN NO;R40 
05/30/95 13:55:24 




175 



TIME(SEC) 



TEST: Model M-2 Test I RUN NO:042 
05/30/95 13:57:01 




176 



TIME(SEC) 



TliST: Model M -2 Test I RIJN NO:043 
05/30/95 15:24:48 



< 
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177 



TIME(SEC) 



TEST: Model M-2 Test I RUN NO:044 
05/30/95 15:26:38 




178 



TIME(SEC) 



TEST: Model M -2 Test I RUN NO045 
05/30/95 15:28:28 




179 



TIMH(SEC) 



I IIS 1 Model M-2 Test I RUN Nl>:<>46 
05/30/95 15:30:15 




180 



TIME(SEC) 
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THST: Model M-3 Test I RUN NO 049 
05/30/95 15:34:28 




s r C 'r, c *0 



182 



TIME(SEC) 



THST : Model M-3 Test I RUN NO:053 
05/30/95 15:36:25 




yr, ~ C- 



183 



DRAG00 053(LB)-AI,L T229/S0 
VEI-OOO 053(FPS)-ALL T229/S0 




TliST: Model M-3 Test I RUN N( ) 055 
05/30/95 15:39 42 




184 



TIME(SEC) 



TUST: Model M-3 Test I RUN NO:056 







185 



DRAGOO 056(LB)-ALL T229/S0 
VELOOO 056(FPS)-ALL T229/S0 



THST; Model M-3 Test I RUN NO 057 
05/30'95 15:50:04 




186 



TIME(SEC) 



TliST: Model M-3 Reverse Direction Data Point RUN NO:058 

05/30/95 16:02:36 




187 



TIME(SEC) 



THST: Model M-3 Test I RUN NO059 
05/30/95 16:04:38 




188 



TIME(SEC) 



TEST:229 RUN NO:R59 




189 



DR AGOO. R59( LB)- A LL T229/SO 
VELOOO R59(FPS)-ALL T229/S0 



THST: Model M-3 lest I RUN NO 061 
05/30/95 16:09:20 



'r, 
r si 




© " -in: 

uuu=i^r“ 

KE 

::r- j _ 
twAVAV^H^ 

iu\\V^ ’**"**' 

GnmViyJVi — 
r 

£s#8T~~ 
mjyf 

• «---■’- 

t . .vlvv- ' 

1 

\ ^vSffe. 

v% s ... 

<* \ *=&s 

</ t \ =JS-" 

~"T— ^ 








190 



TIME(SEC) 



TEST: Model M-3 Test I RIJN NO:R61 
05/30/95 16:11:51 



rsi 
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TIME(SEC) 



I'liS I ' Model M-3 Test I RUN N():063 
05/30/95 16:13:45 




192 



TIME(SEC) 



TEST: Model M-3 Test I RIJN NOR63 
05/30/95 16:15 49 




CnI 



193 



TIME(SEC) 




TliST: Model M-3 Test I RUN NO 066 
05/30/95 16:17:49 




194 



TIME(SEC) 



TEST: Model M-3 Test I RUN NOR66 
05 / 30/05 16 : 10:38 



CN 




R 



195 



TIME(SEC) 



TEST: Model M -3 Test 1 RIJN N();068 
05/30/95 16:21:34 




tr, C 3 it, 



196 



TIME(SEC) 



TEST: Model M-3 Test I RUN NO:069 
05/30/95 16:23:51 




197 



TIME(SEC) 



TliST: Model M-l RUN NO. 7 1 




198 



TIME(SEC) 



mST: Model M-l Draft Variation T-t 0 375" NO:R74 




199 



TIME(SEC) 



TEST: Model M-l Draft Variation RIJN NO:072 
05/30/95 16:31:24 




200 



TIME(SEC) 



TEST: Model M-l Draft Variation T=.25in RUN NO:74 
05/30/95 16:33:46 




"" fN T? 



201 



TIME(SEC) 



TEST: Model M-l Test I RUN NO083 
Model Modified to In-Line Configuration 
05/30/95 16:36:08 
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TIME(SHC) 



Appendix G 



Data Scaling Worksheets 
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This Calculation Scales the Test Data obtained for Model M-1 to Full Scale. 



NOTE: The System of Subscripts in this Routine can become rather bulky... "R", "C", *M", "S" are Resistance, Coefficient 
of Resistance, Measured and Ship respectively... where measured refers to the measured model (Model M-2 Data 
unsealed) and Ship refers to the full scale calculation. Also, "s", "p", V, T, "t“ as subsenpts are strut, pod, residual, 
frictional and total respectively. Other uses of symbols are defined as they arise. 

Input Data Section 



i - 0.. 16 


Sets the data range 


RM t , =READ(mltlfx) 

i 


Total measured X component force for test #t: fixed in pitch and heave. 


VM. =READ(mltlv) 


Speed increments used in the test, (inputted in fps) 


X - 17.26 


Scaling Factor for Lengths 


VM ,. K = VM.- 1 
. 1 1 .688 


Measured Velocity in Knots (converts the inputs which are in fps) 


kts. _ ^^kts.'^ 

i i 


Scaled Velocity in Knots 


VS. =VM ,-a/x 

1 1 V 


Equivalent Scaled Velocity in FPS 


.= 1.2260- 10' 5 


ft*2/sec Kinematic Viscosity of Fresh Water at 59 F (Test Temp) 


u $w = 1 1794 10 5 


ft A 2/sec Kinematic Viscosity of Salt Water at 65 F (Assumed Operating Temp) 


Pfw " 1-9384 


lb sec A 2/ft A 4 Density of Fresh Water at 59 F (Test Temp) 



p SW = 1.9890 lbsec A 2/ft A 4 Density of Salt Water at 65 F (Assumed Operating Temp) 
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Geometry Section 

LM S = 1.39 

SM s =0.715 

LM = 2.14 
SM p = 2.01 

SM t =4 SM s +4 SM 



feet Length of one Strut 



ft A 2 Wetted surface of one Strut 



feet Length of one Pod 



ft A 2 Wetted Surface of one Pod 



ft A 2 Wetted surface of the entire underwater body 



Resistance Calculations 



NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 
Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 
for the model and the ship. 



VM.LM s 

RllM g operating Rn of the Strut 

1 v fw 



VM.LM n 
RnM = ! il 

1 u f\v 



operating Rn of the Pod 



CM 



fs. 



0.075 




Coefficient of Frictional Resistance based on Strut Length 
and using the 1957 ITTC line. 



CM 






0.075 

(logpnMpJ-2) 2 



Coefficient of Frictional Resistance based on Pod Length 
and using the 1957 ITTC line. 
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Frictional Resistance of four Struts 



RMfs. =CM fs.YPf\v( V M i ) 2 -4-SM s 



RM fp - CM j-p •— p ^VM.j^-4SM p Frictional Resistance of four Pods 



RM r =RM t , - RM f$ - RM fp 



RM r 

CM r =- 



Residual Resistance for the Model as measured. 



RnM s 

i 



137186 


273238 


412692 


459176 


462577 


504527 


530604 


549878 


595228 


620171 


641713 


733548 


916085 


1 1 46240 


1375261 


1603148 


1834437 



RnM 

Pi 



21 1207 


420669 


635367 


706933 


712170 


776754 


816900 


846574 


916395 


954796 


987961 


1129347 


1410375 


1764715 


21 17308 


2468157 


2824241 



RM ^ j RM 



0.151 


0.031 


0.689 


0.102 


4.24 


0.211 


6.29 


0.254 


6.46 


0.258 


8.91 


0.3 


9.43 


0.328 


9.08 


0.35 


6.86 


0.402 


6.16 


0.432 


5.65 


0.459 


5.61 


0.582 


7.18 


0.865 


10.645 


1.29 


19.465 


1.786 


28.82 


2.351 


35.35 


2.994 



RM f RM r 

i 



0.043 


0.328 


3.494 


5.389 


5.547 


7.845 


8.265 


7 839 


5.432 


4.624 


4.017 


3.538 


4.099 


6.043 


13.084 


20.413 


24.633 



0.077 


0.258 


0.535 


0.647 


0.655 


0.765 


0.836 


0.891 


1.026 


1.104 


1.173 


1.489 


2.216 


3.312 


4.595 


6.057 


7.723 



CM fs. CM fp. 

I r l 



0.008 


0.006 


0.006 


0.006 


0.006 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.004 


0.004 


0.004 



0.007 


0.006 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.004 


0.004 


0.004 


0.004 


0.004 



Note: these values for Rn are Low for ITTC 1957 use. Should be on order of 10 A 7 



CM r 

i 

0.003 

0.005 

0.025 

0.031 

0.032 

0.038 

0.036 

0.032 

0.019 

0.015 

0.012 

0.008 

0.006 

0.006 

0.008 

0.01 

0.009 
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Scaled Resistance Calculations 



CS j. - CM r Residual Resistance Coefficient remains the same in the saling 

i i 



RnS s . 



VS LM -X 

1 ^ 



U 



s\v 



operating Rn of the Start in the new scale 



RnS 



P, 




operating Rn of the Pod in the new scale 




i 



0.075 

tog, RnS s 'j - 2 



Coefficient of -Frictional Resistance based on scaled Start Length 
and using the 1957 ITTC line. 



CS 



fp i 



0.075 



log/ RnS p - 2\ 2 



Coefficient of Frictional Resistance based on scaled Pod Length 
and using the 1957 ITTC line. 



CS„ - 0.0004 

d 



Correlation Allowance 



RS is. 

i 


i z 


Frictional Resistance of the Struts 


RS f Pi 


= CS f Pi |p S w( v S,) 24 SM p ». 2 


Frictional Resistance of the Pods 


RS r 

i 


= r ^ p sw (VS.) 2 SM | X 2 


Scaled Residual Resistance 


RS a ( 


= CS a | Psw ( v; 5.) 2 SM, >. 2 


Correlation 


RS, 

i 


= RS Is. * RS fp. + RS r. + RS a 

i * l i i 
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Scaled 


Scaled 


Total M-2 


Speed 


Speed 


Resistance 


(fps) 


(knots) 


at Full 






Scale 


VS. 


vs kts. 

i 


RS t 

i 



5.03 


10.01 


15.12 


16.83 


16.95 


18.49 


19.44 


20.15 


21.81 


22.73 


23.51 


26.88 


33.57 


42 


50.39 


58.74 


67.22 



2.978 


5.932 


8.959 


9.968 


10.042 


10.952 


1 1.518 


11.937 


12.921 


13.463 


13.93 


15.924 


19887 


24.883 


29.854 


34.801 


39.822 



Residual 

Component 



RS, 



PodFrictional Strut Frictional 
Component Component 



488.89 


2680.46 


20480.51 


30929.77 


3 1798.85 


44370.38 


46882.96 


44860.57 


32717.51 


28777.95 


25868.95 


24672.03 


30738.84 


45766.63 


88587.27 


133802.96 


163618.18 



225.11 


1732.1 


18434.19 


28431.8 


29266.18 


41391.37 


43609.44 


41361.02 


28658.3 


24394.36 


21 195.81 


18667.49 


21627.81 


31883.98 


69030.12 


107699.01 


129965.39 



RS 



f P; 



167.14 


592.83 


1268.39 


1544.88 


1566.09 


1838.72 


2018.35 


2156.07 


2496.66 


2693.8 


2869.63 


3676.58 


5551.69 


8417.37 


1 1809.69 


15708.47 


20189.01 



RS 



fs. 



Reynolds # 
of the Strut 



RnS , 



Reynolds # 
of the Pod 



63.99 


226.03 


482.53 


587.39 


595.44 


698.78 


766.85 


819.03 


948.04 


1022.69 


1089.26 


1394.66 


2103.76 


3186.42 


4466.99 


5937.76 


7627.04 



10225869 


20367227 


307621 19 


34227083 


34480617 


37607536 


39551296 


40987988 


44368441 


46227690 


47833405 


54678821 


68285143 


85440940 


102512226 


119499001 


136739309 



RnS 



Pi 



15743425 


31356738 


47360385 


52694933 


53085266 


57899371 


60891923 


63103809 


68308247 


71170688 


73642796 


84181783 


10512964 


13154216 


15782457 


18397687 


21051951 



WRITE(mlfulfx) -RS^ WRITE(m lfulv) - VS j^ s Writes the scaled data to a file for future use. 
i i 

RS, 



R s fl =RSf p.fRS fs _ 

1 r l I 



FF. =- 



RS ft. 



' RS, 



RF. 



RS 



— Percent of total resistance 

[ attributable to friction/wave, 

i 



RS 



CA. = 

' RS, 



a. 



percent of total resistance that is attributable to correlation allowance 



WRITE(mlfulfp) =FF WRITE( m I fulrp ) = RF. WRITE(mlfulcap) = CA. 
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This Calculation Scales the Test Data obtained for Model M-2 to the scale that Model M-1 
is built to. This enables a direct comparison of the Data obtained from the two models. 



NOTE: The System of Subscripts in this Routine can become rather bulky... "R“, "C". *M“, "S" are Resistance, Coefficient 
of Resistance, Measured and Scaled respectively... where measured refers to the measured model (Model M-2 Data 
unsealed) and Scaled refers to the Model M-2 data scaled to the Model M-1 scale factor. Also, "s", “p\ V, T, "t" as 
subscripts are strut, pod, residual, frictional and total respectively. Other uses of symbols are defined as they arise. 

Input Data Section 



i :=0.. 15 


Sets the data range 


RM t , = READ(m2tlfx) 
i 


Total measured X component force for test #1: fixed in pitch and heave. 


VM. =READ(m2tlv) 


Speed increments used in the test, (inputted in fps) 


_ 16.68 
17.26 


Scaling Factor for Lengths M2/M1 scale factor. 


VM],., = VM.- ‘ 

. 1 1 .688 


Measured Velocity in Knots (converts the inputs which are in fps) 


VSkts. 

1 1 


Scaled Velocity in Knots 


VS. = VM.Jx 

1 1 V 


Equivalent Scaled Velocity in FPS 


u hv = 1.2260 10 5 


ft A 2/sec Kinematic Viscosity of Fresh Water at 59 F (Test Temp) 


u sw — 1.17 94 - 1 O' 5 


ft A 2/sec Kinematic Viscosity of Salt Water at 65 F (Assumed Operating Temp) 


Pfw = 1-9384 


lb sec A 2/ft A 4 Density of Fresh Water at 59 F (Test Temp) 


P sw •= 1.9890 


lb sec A 2/ft A 4 Density of Salt Water at 65 F (Assumed Operating Temp) 
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Geometry Section 

LM S = 1.4375 

SM s = 1 059 

LM p = 2.32 
SM p = 2.642 

SM, -4 SM s - 4 SM p 



feet Length of one Strut 



ft A 2 Wetted surface of one Strut 



feet Length of one Pod 



ft A 2 Wetted Surface of one Pod 



ft A 2 Wetted surface of the entire underwater body 



Resistance Calculations 



NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1 957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 
Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 
for the model and the ship. 



VM.-LM s 

RnM = operating Rn of the Strut 

' u f\v 



VM.-LM _ 
RnM = 1 E 

'■ u fw 



operating Rn of the Pod 



CM 



fs. 

I 



0.075 




Coefficient of Frictional Resistance based on Strut Length 
and using the 1957 ITTC line. 



CM fPi 



0.075 

(log(RnM pJ - 2 2 



Coefficient of Frictional Resistance based on Pod Length 
and using the 1957 ITTC line. 
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RMfs. 

l l z 



Frictional Resistance of four Struts 



RM 



RM r 



CM r. 



RnM s 

i 

143046 

288438 

431485 

481903 

524113 

553426 

578049 

612051 

641364 

699990 

771513 

964978 

1205343 

1448053 

1688418 

1928783 



>. = c Mrp.-P h ,(VM i ) 2 4 SMp 



Frictional Resistance of four Pods 



= RM j j - RM f s - RM 

i i * i 



Residual Resistance for the Model as measured. 



RM r 

i 

^Pfv(VM,! 2 SM t 



RnM 

* i 



230865 


465514 


696378 


777749 


845873 


893181 


932920 


987798 


1035106 


1129723 


1245155 


1557390 


1945318 


2337031 


2724959 


3112887 



RMu 

1 



0.172 


0.673 


4.52 


6.83 


9.81 


10.78 


9.84 


7.85 


6.68 


6.05 


6.29 


7.255 


10.21 


14.49 


20.18 


26.39 



RM fs RM f RM r 

i ‘ i i 



0.101 


0.346 


0.704 


0.857 


0.995 


1.096 


1.185 


1.312 


1.426 


1.667 


1.983 


2.96 


4.413 


6.139 


8.097 


10.297 



0046 


0.156 


0.316 


0.383 


0.445 


0.49 


0.529 


0.585 


0.636 


0.742 


0.882 


1.314 


1.954 


2.713 


3.573 


4 537 



0.025 


0.172 


3.5 


5.589 


8.37 


9.194 


8.126 


5.953 


4.619 


3.641 


3.425 


2.981 


3.843 


5.638 


8.51 


1 1.556 



CM j' s CM f CM r 

1 r i 1 



0.008 


0.006 


0.006 


0.006 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.004 


0.004 


0.004 



0.007 


0.006 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 



0.001 


0.002 


0.018 


0.023 


0.029 


0.029 


0.023 


0.015 


0.01 1 


0.007 


0.006 


0.003 


0.003 


0.003 


0.003 


0.003 





Note: these values for Rn are Low for ITTC 1957 use. Should be on order of 10 A 7 



Scaled Resistance Calculations 



CS r - CM r Residual Resistance Coefficient remains the same in the saling. 

i i 



RnS , 



VS.-LM ^k 



* u 



fw 



RnS 



VS.-LM n 
i P 



P; 



1 u 



fvv 



CS 



fs 



0.075 



1 (' log I RnS s 



t 



CS 



f p, 



0.075 



log/ RnS p - 2 



\ 2 



cs a =0 



operating Rn of the Strut in the new scale 



operating Rn of the Pod in the new scale 



Coefficient of Frictional Resistance based on scaled Strut Length 
and using the 1957 ITTC line. 



Coefficient of Frictional Resistance based on scaled Pod Length 
and using the 1957 ITTC line 



Correlation Allowance 



RS fs ^ fs.' . 'P fw 



VS. I 2 4 SM _.}. 2 

1 1 S 



Frictional Resistance of the Struts 



Rs fr, =cs fp. ^ Pf\v vs.) 2 j sm p ). 2 

RS , ^'Sr ^Pfc.VS^SM ,! 2 

i i 2 



Frictional Resistance of the Pods 



RS 



t. 



RS 



fs. 



RS 



f Pi 



RS 



r. 

i 



Scaled Residual Resistance 



Scaled 

Speed 

(fps) 


Scaled 

Speed 

(knots) 


Total M-2 
Resistance 
at 17.26:1 
Scale 


Residual 

Component 


PodFrictional 

Component 


Strut Frictional 
Component 


Reynolds # 
of the Strut 


Reynolds # 
of the Pod 


vs. 


VS kts 

1 


RS t 

i 


RS r 

i 


RS f Pi 


RS fs. 

1 


^ s. 

1 


RnSp 



135897 


274022 


409919 


457817 


497918 


525766 


549158 


581461 


609309 


665005 


732953 


916748 


1 145100 


1375679 


1604031 


1832383 



0.711 


1.433 


2.143 


2.394 


2.603 


2.749 


2.871 


3.04 


3.186 


3.477 


3.832 


4.793 


5.987 


7.192 


8.386 


9.58 



0.16 


0.61 


4.09 


6.18 


8.87 


9.75 


89 


7.1 


6.05 


5.48 


5.71 


6.59 


9.28 


13.16 


18.32 


23.95 



0.04 


0.14 


0.29 


0.35 


0.41 


0.45 


0.48 


0.53 


0.58 


0.68 


0.81 


1.2 


1.78 


2.47 


3.26 


4.14 



1.2 


2.42 


3.62 


4.04 


4.39 


4.64 


4.85 


5.13 


5.38 


5.87 


6.47 


8.09 


10.1 1 


12.14 


14.16 


16.17 



0.09 


0.32 


0.64 


0.78 


0.91 


1 

1.08 

1.2 


1.3 


1.52 


1.81 


2.7 


4.02 


5.6 


7.38 


9.39 



0.02 


0.16 


3.16 


5.04 


7.55 


8.3 


7.33 


5.37 


4.17 


3.29 


3.09 


2.69 


3.47 


5.09 


7.68 


10.43 



219326 


442248 


661574 


738877 


803596 


848540 


886293 


938428 


983372 


1073259 


1 182922 


1479552 


1 848092 


2220227 


2588767 


2957306 



WRITE(m2sclfx) -RS^ WRITE(m2sclv) = VS, Writes the scaled data to a file for future use. 



RS 



RSf| RSfp RS fs. 

i 1 i i 



FF. =- 



ft. 



1 RS, 



RF. = 1 - FF. Percent of total resistance 

i i 

attributable to friction/wave. 



WRITE(m2fT) = FF WRITE(m2rf) =RF. 




This Calculation Scales the Test Data obtained for Model M-2 to Full Scale . 



NOTE: The System of Subscripts in this Routine can become rather bulky... "R", "C". "M", "S’ are Resistance, Coefficient 
of Resistance, Measured and Ship respectively... where measured refers to the measured model (Model M-2 Data 
unsealed) and Ship refers to the full scale calculation. Also, "s’, "p", V, T, "t" as subscripts are strut, pod, residual, 
frictional and total respectively. Other uses of symbols are defined as they arise. 

Input Data Section 



II 

o 


Sets the data range 


RM t , = READ(m2tlf\) 
i 


Total measured X component force for test #1 : fixed in prtch and heave. 


VM. =READ(m2tlv) 


Speed increments used in the test, (inputted in fps) 


/. = 16.68 


Scaling Factor for Lengths M2/M1 scale factor. 


VM - VM.- ' 

i ' 1 .688 


Measured Velocity in Knots (converts the inputs which are in fps) 


VS kts. =VM ktt .-^ 

1 1 


Scaled Velocity in Knots 


VS. =VM.-v. 
1 1 M 


Equivalent Scaled Velocity in FPS 


Vfn - 1.2260- 10 5 


ft A 2/sec Kinematic Viscosity of Fresh Water at 59 F (Test Temp) 


Usw = 1 .1794 10 5 


ft A 2/sec Kinematic Viscosity of Salt Water at 65 F (Assumed Operating Temp) 


Pfw =1 9384 


lb sec A 2/ft A 4 Density of Fresh Water at 59 F (Test Temp) 



p gyy - 1 .9890 lb sec*2/ft*4 Density of Salt Water at 65 F (Assumed Operating Temp) 
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Geometry Section 



LM S 


= 1.4375 


feet Length of one Strut 


SM s 


= 1.059 


ft A 2 Wetted surface of one Strut 


LM p 


= 2.32 


feet Length of one Pod 


SM p 


= 2.642 


ft A 2 Wetted Surface of one Pod 



SM ( - 4 • SM <5+4 SM p ft A 2 Wetted surface of the entire underwater body 



Resistance Calculations 



NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1 957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 
Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 
for the model and the ship. 



VM.-LM s 

RnM s operating Rn of the Strut 

°f\v 



VM,LM n 
RnM = 1 l 

u fvv 



operating Rn of the Pod 



CM 



fs. 



0.075 




Coefficient of Frictional Resistance based on Strut Length 
and using the 1957 ITTC line. 



CM 






0.075 

log [ RnM p j - 2 2 



Coefficient of Frictional Resistance based on Pod Length 
and using the 1957 ITTC line. 
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RMfs. = CM fs.-|piW( VM i ) 2 - 4 - SM ! 



Frictional Resistance of four Struts 



RM 



<Pi 




Frictional Resistance of four Pods 



RM r =RM U - RM fs - RM f 

i i i *1 

RM r 

CM , ‘ 

jfpfw- vm/.smj 

L z J 



Residual Resistance for the Model as measured 



RnM s 

i 



143046 


288438 


431485 


481903 


524113 


553426 


578049 


612051 


641364 


699990 


771513 


964978 


1205343 


1448053 


1688418 


1928783 



RnM 

*1 



230865 


465514 


696378 


777749 


845873 


893181 


932920 


987798 


1035106 


1129723 


1245155 


1557390 


1945318 


2337031 


2724959 


3112887 



1 



0.172 


0.673 


4.52 


6.83 


9.81 


10.78 


9.84 


7.85 


6.68 


6.05 


6.29 


7.255 


10.21 


14.49 


20.18 


26.39 



RM fs RM f RM r 

i *i i 



0.046 


0.156 


0.316 


0.383 


0.445 


0.49 


0.529 


0.585 


0.636 


0.742 


0.882 


1.314 


1.954 


2.713 


3.573 


4.537 



0.101 


0.346 


0.704 


0.857 


0.995 


1.096 


1.185 


1.312 


1.426 


1.667 


1.983 


2.96 


4.413 


6.139 


8.097 


10.297 



0.025 


0.172 


3.5 


5.589 


8.37 


9.194 


8.126 


5.953 


4.619 


3.641 


3.425 


2.981 


3.843 


5.638 


8.51 


1 1.556 



CM fs. CM fp. 

I r ) 



0.008 


0.006 


0.006 


0.006 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.004 


0.004 


0.004 



0.007 


0.006 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 



CM f 

i 

0.001 

0.002 

0.018 

0.023 

0.0291 

0.029 

0.023 

0.015 

0.01 1 

0.007 

0.006 

0.003 

0.003 

0.003 

0.003 

0.003 
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Note: these values for Rn are Low for ITTC 1957 use. Should be on order of 10 A 7 



Scaled Resistance Calculations 



CS r CM r Residual Resistance Coefficient remains the same in the saling. 

i i 



FLn S s 



VS.-LM .-X 

I 



u 



s\v 



RnS 



Pi 



VS.-LM -J. 
1 P 



u 



s\v 



operating Rn of the Strut in the new scale 



operating Rn of the Pod in the new scale 



cs 



fs. 



0.075 

log|RnS s \- 2\ 2 



Coefficient of’Frictional Resistance based on scaled Strut Length 
and using the 1957 ITTC line. 



CS 




0.075 

log RnS > - 2\ 2 



Coefficient of Frictional Resistance based on scaled Pod Length 
and using the 1957 ITTC line. 



CS a - 0.0004 Correlation Allowance 



RS fs. = cs fs.-r-Pswi vs i ) 2 - 4 - SM s ^ 2 

1 I a- 


Frictional Resistance of the Struts 


RS fp. ^S {p l. Psw -{ys.f.4.SM p .X 2 


Frictional Resistance of the Pods 


RS Fj =CS r ip sw -(VS.) 2 SM t -L 2 


Scaled Residual Resistance 


RS a. =CS a ip sw .(VS i ) 2 .SM ( L 2 

i 2 


Correlation 


RS t. = RS fs.- RS fp. + RS r.- RS a. 

i i r i i i 
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Scaled 

Speed 

(fps) 


Scaled 

Speed 

(knots) 


Total M-2 
Resistance 
at Full 
Scale 


Residual 

Component 


PodFnctional 

Component 


StrutFnctional 

Component 


Reynolds # 
of the Strut 


Reynolds # 
of the Pod 


vs. 


VS kts. 

1 


RS t 

i 


RS r 


RS f P . 


RS fs. 

1 







200.31 


727.02 


1528.5 


1874.83 


2189.74 


2421.67 


2624.82 


2917.76 


3181.78 


3741.42 


4480.68 


6785.95 


10257.03 


14426.78 


19199.71 


24600.85 



446.65 


2025.26 


19224.96 


29761.5 


43537.84 


47855.66 


431 17.98 


33268.15 


27363.12 


23664.57 


23900 


25749.59 


35850.64 


51632.94 


73622.44 


97565.13 



2.952 


5.952 


8.904 


9.944 


10.815 


1 1.42 


1 1.928 


12.63 


13.235 


14.444 


15.92 


19.912 


24.872 


29.881 


34.841 


39.801 



118.56 


818.18 


16666.37 


26616 


39857.57 


43780.95 


38697.41 


28348.32 


21992.81 


17337.95 


16307.54 


14195.24 


18300.87 


26848. 1 8 


40524.35 


55028.49 



10129794 


2042565 1 


30555445 


34125783 


37114902 


39190680 


40934333 


43342235 


45418013 


49569568 


54634465 


68334597 


85355972 


102543410 


1 19564786 


136586162 



16348607 


32965225 


49313832 


55076046 


59900225 


63250349 


66064454 


69950598 


73300723 


80000972 


88175275 


110286097 


137757117 


165496147 


192967168 


220438188 



87.1 


314.68 


659.99 


809.02 


944.46 


1044.18 


1131.5 


1257.39 


1370.81 


1611.16 


1928.51 


2917.3 


4404.6 


6189.62 


8231.34 


10540.35 



4.98 


10.05 


15.03 


16.79 


18.26 


19.28 


20.13 


21.32 


22.34 


24.38 


26.87 


33.61 


41.98 


50.44 


58.81 


67.18 



WRITE(m2fulfx) =RS t 



WRITE(m2fulv) = VS kts 

i 



Writes the scaled data to a file for future use. 




i 




RS 



fs. 




I 



RF 

I 




1 



RS 



t. 



Percent of total resistance 
attributable to friction/wave. 



CA. 

1 



RS a 




percent of total resistance that is attributable to correlation allowance 



WRJTE( m2fulfp) =FF. WRITE(m2ftilrp) =RF. WRITE(m2fulcap) = CA. 



EHP 



ship. 



RS , VS. 

» • l 

I 

550 
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VS ^ EHP s hip 

l r I 



4.983 


10.047 


15.03 


16.786 


18.256 


19.277 


20.135 


21.319 


22.34 


24.382 


26.873 


33.612 


41.985 


50.439 


58.81 1 


67.184 



4.046 


36.996 


525.35 


908.305 


1445.139 


1677.3 


1578.486 


1289.539 


1 1 1 1.446 


1049.079 


1 167.775 


1573.641 


2736.688 


4735.1 


7872.416 


11917.794 



2.952 


5.952 


8.904 


9.944 


10.815 


11.42 


11.928 


12.63 


13.235 


14.444 


15.92 


19.912 


24.872 


29.881 


34.841 


39.801 



M-2 in full scale 
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This Calculation Scales the Test Data obtained for Model M-3 to the scale that Model M-1 
is built to. This enables a direct comparison of the Data obtained from the two models. 



NOTE: The System of Subscripts in this Routine can become rather bulky... "R", "C“, "M", "S" are Resistance, Coefficient 
of Resistance, Measured and Scaled respectively... where measured refers to the measured model (Model M-3 Data 
unsealed) and Scaled refers to the Model M-3 data scaled to the Model M-1 scale factor. Also, "s", "p", V, T, T as 
subscripts are strut, pod, residual, frictional and total respectively Other uses of symbols are defined as they arise. 

Input Data Section 



i r 0.. 16 


Sets the data range 


RM t] =READ(m3tlfx) 

i 


Total measured X component force for test #1: fixed in pitch and heave. 


VM. =READ(m3tlv) 


Speed increments used in the test, (inputted in fps) 


. _ 15.1 
17.26 


Scaling Factor for Lengths M3/M1 scale factor. 


VM kl , -VM.- 1 

i ' 1.688 


Measured Velocity in Knots (converts the inputs which are in fps) 


VS kts =VM kts.V^ 
1 1 


Scaled Velocity in Knots 


VS. =VM. a ,/. 


Equivalent Scaled Velocity in FPS 


= 1.2260- 10 5 


ft A 2/sec Kinematic Viscosity of Fresh Water at 59 F (Test Temp) 


u sw , = 1.1794 10' 5 


ft A 2/$ec Kinematic Viscosfty of Salt Water at 65 F (Assumed Operating Temp) 


Pfw -- 1.9384 


lb sec A 2/ft A 4 Densrty of Fresh Water at 59 F (Test Temp) 



P SW ~ 1 .9890 lb sec*2/ft*4 Density of Salt Water at 65 F (Assumed Operating Temp) 
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Geometry Section 



LM S = 1.589 

SM s - 1.716 

LM p =5.958 
SM p =6.45 

SM t =4 SM s 



2 SM p 



feet Length of one Strut 



ft A 2 Wetted surface of one Strut 



feet Length of one Demihull 



ft A 2 Wetted Surface of one Demihull 



ft A 2 Wetted surface of the entire underwater body 



Resistance Calculations 



NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 
Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 
for the model and the ship. 



VM,LM S 

RnM „ = operating Rn of the Strut 

u f\v 



RnM 



Pi 



VM. LM_ 
i P 



u hv 



operating Rn of the Demihull 



CM 



0.075 



fs. 



log[RnM \-2) 



Coefficient of Frictional Resistance based on Strut Length 
and using the 1957 ITTC line. 



CM 



*Pi 



0.075 



logi RnM - 2 
"i / 



Coefficient of Frictional Resistance based on Demihull Length 
and using the 1957 ITTC line. 
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RM 



fSj 



RM 



*1 



= CM fs .i p^, VM ) 2 -4 SM s 
* 2 7 

^ CM fp. jPfu- VM.) 2 2.SM p 



Frictional Resistance of four Struts 



Frictional Resistance of two Demihulls 



RM r =RM t , - RM fs - RM f 

i i i *i 

RM r 

CM r = 7 , 

' I Pftv . VM i 2 SM t( 



Residual Resistance for the Model as measured. 



RnM s 

i 



334390 


451038 


502881 


540467 


593607 


619529 


644154 


671372 


773763 


896243 


1008354 


11211 13 


1 289604 


1401068 


1684910 


1960976 


2240931 



RnM 

* i 



1253804 


1691178 


1885566 


2026498 


2225746 


2322940 


2415274 


2517328 


2901245 


3360487 


3780852 


4203646 


4835408 


5253343 


6317618 


7352736 


8402432 



RM t j RM 

i i 



0.801 


3.17 


4.22 


6.01 


7.14 


6.75 


6.35 


6.08 


6.74 


8.23 


8.93 


9.74 


12.2 


13.865 


19.08 


25.875 


33.185 



0.267 


0.453 


0.548 


0.623 


0.735 


0.793 


0.849 


0.914 


1.176 


1.527 


1.884 


2.276 


2.924 


3.391 


4.72 


6.199 


7.881 



RM fp RM r 



0.372 


0.635 


0.772 


0.879 


1.041 


1.124 


1.206 


1.299 


1.678 


2.188 


2.709 


3.282 


4.231 


4.918 


6.877 


9.064 


1 1.56 



0.162 


2.082 


2.899 


4.508 


5.364 


4.833 


4.295 


3.867 


3.886 


4.514 


4.337 


4.182 


5.046 


5.556 


7.483 


10.612 


13.744 



CM fs CM fp 

1 r I 



0.006 


0.006 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.004 


0.004 


0.004 


0.004 


0.004 



0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.003 


0.003 


0.003 


0.003 


0.003 



CM r 

i 

0.001 
0.009 
0.01 
0.014 
0.013 
0.01 1 
0.009 
0.008 
0.006 
0.005 
0.004 
0.003 
0.003 
0.002 
0.002 
0.002 
0.002 



Note: these values for Rn are Low for ITTC 1957 use. Should be on order of 10 A 7 




Scaled Resistance Calculations 



CS r = CM r 



Residual Resistance Coefficient remains the same in the sating. 



RnS 



RnS 



VS. LM 

_ 1 ^ 

s. 

' u f\v 

VS. LM „■/. 



Pi 



■ P 
u fw 



operating Rn of the Strut in the new scale 



operating Rn of the Pod in the new scale 



cs 



0.075 



fs. 



log RnS j - 2 



Coefficient of Frictional Resistance based on scaled Strut Length 
and using the 1957 ITTC line. 



CS 






0.075 



i log RnS j - 2 

\ \ 



2 



Coefficient of Frictional Resistance based on scaled Pod Length 
and using the 1957 ITTC line. 



cs a =0 



Correlation Allowance 



RS fs. 

1 


= cs f».-7 
1 . 


Rs fp, 


II 

n 

Zfi 

J? 


RS r. 

1 


= cs r -- 

r i 2 


RS t 


= RS fs. + ' 



cm 2 



/s.) 2 sM t ;. 2 



.2 



.2 



Frictional Resistance of the Struts 



■A. Frictional Resistance of the Demi Hulls 



Scaled Residual Resistance 



i i * i i 
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Scaled 


Scaled 


Total M-3 


Speed 


Speed 


Resistance 


(f ps) 


(knots) 


at 17.26:1 






Scale 


VS. 


VS kts. 

1 


RS, 

i 



1.43 


1.928 


2.15 


2.31 1 


2.538 


2.649 


2.754 


2.87 


3.308 


3.832 


4.31 1 


4.793 


5.513 


5.99 


7.203 


8.384 


9.581 



0.56 


2.16 


2.87 


4.07 


4.83 


4.58 


4.31 


4.14 


4.59 


5.62 


6.1 1 


6.67 


8.36 


9.51 


13.08 


17.72 


22.72 



2.41 


3.25 


3.63 


3.9 


4.28 


4.47 


4.65 


4.85 


5.58 


6.47 


7.28 


8.09 


9.31 


10.11 


12.16 


14.15 


16.17 



Residual PodFrictional 
Component Component 



0.11 




0.26 


1.39 




0.44 


1.94 




0.54 


3.02 




0.61 


3.59 




0.73 


3.24 




0.78 


2.88 




0 84 


2.59 




0.91 


2.6 




1 17 


3.02 




1.52 


2.9 




1.88 


2.8 




2.28 


3.38 




2.94 


3.72 




3.42 


5.01 




4.78 


7.11 




6.29 


9.2 
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WRITE(in3sclfx') = RS, WRITE(m3sclv) = VS. 
R - ^ ft. -^fp.~^fs. 

I f 1 I 

WRITE(m3ff) - FF. WRITE(m3rf) =RF. 



RS 



FF. =■ 



ft. 



RS, 



StrutFrictional Reynolds # Reynolds # 

Component of the Strut of the Pod 



RnS s 

1 

273626 

369077 

411500 

442256 

485739 

506951 

527101 

549373 

633158 

733381 

825120 

917390 

1055263 

1 146472 

1378736 

1604636 

1833718 



r l 



1025968 


1383864 


1542929 


1658251 


1821293 


1900825 


1976381 


2059890 


2374043 


2749834 


3093812 


3439778 


3956738 


4298728 


5169608 


6016628 


6875579 



RSfs 



0.19 


0.32 


0.38 


0.44 


0.52 


0.56 


0.6 


0.64 


0.82 


1.07 


1.32 


1.59 


2.04 


2.37 


3.29 


4.32 


5.49 



Writes the scaled data to a file for future use. 



RF. = 1 - FF. 

I I 



Percent of total resistance 
attributable to friction/wave. 
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This Calculation Scales the Test Data obtained for Model M-3 to full scale. 



NOTE: The System of Subscripts in this Routine can become rather bulky... "R", "C", "M", "S" are Resistance, Coefficient 
of Resistance, Measured and Scaled respectively... where measured refers to the measured model (Model M-3 Data 
unsealed) and Scaled refers to the Model M-3 data scaled to the Model M-1 scale factor. Also, "s", "p", V, T, "t" as 
subscripts are strut, pod, residual, frictional and total respectively Other uses of symbols are defined as they arise. 

Input Data Section 



II 

o 

ON 


Sets the data range 


RM t , =READ(m3tlfx) 
i 


Total measured X component force for test #1: fixed in pitch and heave. 


VM. =READ(m3llv) 


Speed increments used in the test (inputted in fps) 


X = 15.1 


Scaling Factor for Lengths M3/M1 scale factor. 


VM ,. tq - VM, ' 

i 1 1 688 


Measured Velocity in Knots (converts the inputs which are in fps) 


VS kts. =VMkts.-;». 

1 1 


Scaled Velocity in Knots 


f— 

VS = VM.a'/. 

1 1 


Equivalent Scaled Velocity in FPS 


u fw = 1.2260 10' 5 


ft A 2/sec Kinematic Viscosity of Fresh Water at 59 F (Test Temp) 


u sw = 1.1794 MO 5 


ft A 2/sec Kinematic Viscosity of Salt Water at 65 F (Assumed Operating Temp) 


Pfw =1.9384 


lb sec A 2/ft A 4 Density of Fresh Water at 59 F (Test Temp) 


P S w = 1-9890 


lb sec A 2/ft A 4 Density of Salt Water at 65 F (Assumed Operating Temp) 
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Geometry Section 

LM S = 1.589 

SM s = 1.716 

LM - 5.958 
SM p = 6.45 

SM t = 4 SM s - 2 • SM p 



feet Length of one Strut 



ft A 2 Wetted surface of one Strut 



feet Length of one Demihull 



ft A 2 Wetted Surface of one Demihull 



ft A 2 Wetted surface of the entire underwater body 



Resistance Calculations 



NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 
Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 
for the model and the ship. 



VM,LM S 

RnM _ '■ operating Rn of the Strut 

' j u fu 



VM.-LM n 
RnM = 1 £ 

' u t\v 



operating Rn of the Demihull 



CM 



fs. 



0.075 

log (RnM s - 2 



Coefficient of Frictional Resistance based on Strut Length 
and using the 1957 ITTC line. 



CM 



f Pi 



0.075 

1 logfRnM p ^ - 2 



Coefficient of Frictional Resistance based on Demihull Length 
and using the 1957 ITTC line. 
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RM 



RM 



RM r 



CM r 



RnM s 
" i 

334390 

451038 

502881 

540467 

593607 

619529 

644154 

671372 

773763 

896243 

1008354 

1 121 1 13 

1289604 

1401068 

1684910 

1960976 

2240931 



, = ^ fe.-J-P fW v VM.) 2 -4-SM s 



Frictional Resistance of four Struts 



= CM 



fp ' p fw 



• VM. 2 -2 SM, 



Frictional Resistance of two Demihulls 



: RM ^ j — RM — RM 



Residual Resistance for the Model as measured. 



RM, 



P f\ V ' ( VM. ) 2 SM t 



RnM p 



1253804 


1691178 


1885566 


2026498 


2225746 


2322940 


2415274 


2517328 


2901245 


3360487 


3780852 


4203646 


4835408 


5253343 


6317618 


7352736 


8402432 



RM U RM fs 

I I 



0.801 


0.267 


3.17 


0.453 


4.22 


0.548 


6.01 


0.623 


7.14 


0.735 


6.75 


0.793 


6.35 


0.849 


6.08 


0.914 


6.74 


1.176 


8.23 


1.527 


8.93 


1.884 


9.74 


2.276 


12.2 


2.924 


13.865 


3.391 


19.08 


4.72 


25.875 


6.199 


33.185 


7.881 



RM t -p RM r 



0.372 


0.635 


0.772 


0.879 


1.041 


1.124 


1.206 


1.299 


1.678 


2.188 


2.709 


3.282 


4.231 


4.918 


6.877 


9.064 


1 1.56 



0.162 


2.082 


2.899 


4.508 


5.364 


4.833 


4.295 


3.867 


3.886 


4.514 


4.337 


4.182 


5.046 


5.556 


7.483 


10.612 


13.744 



CM f s CM f CM r 

I I I 



0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.003 


0.003 


0.003 


0.003 


0.003 



0.001 


0.009 


0.01 


0.014 


0.013 


0.01 1 


0.009 


0.008 


0.006 


0.005 


0.004 


0.003 


0.003 


0.002 


0.002 


0.002 


0.002 



0.006 


0.006 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.005 


0.004 


0.004 


0.004 


0.004 


0.004 
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Note, these values for Rn are Low for ITTC 1957 use. Should be on order of 10 A 7 



Scaled Resistance Calculations 



CS r - CM r Residual Resistance Coefficient remains the same in the saling. 



RjlS 



RjlS 



VS.-LM <- /. 

I ^ 



s. 

1 u 



sw 

VS.-LM n /. 
i P 



P, 



u 



operating Rn of the Strut in the new scale 



operating Rn of the Pod in the new scale 



SW 



CS 



fs. 



0.075 



1 log/RnS s ( - 2 
i / 



Coefficient of Frictional Resistance based on scaled Strut Length 
and using the 1957 ITTC line. 



cs 



fPj 



0.075 



log(RnS - 2 



i / 



Coefficient of Fnctional Resistance based on scaled Pod Length 
and using the 1957 ITTC line. 



CS a =0.0004 



Correlation Allowance 



i i i 2 , . j 

RS r CSj' s — p s ^ v VS 4 SM s 7 Fnctional Resistance of the Struts 

i '12 1 

| , 2 2 

RS ^ r CS - p VS. • 2 SM p 7. Frictional Resistance of the Demi Hulls 



1 



RS r - CS r •— p sw - 1 VS. j 2 -SM t /. Scaled Residual Resistance 

i i 2 " ' X 

RS a = CS a — P sw '( VS/) 2 SM { K Correlation 

i 7 1 



RS t = RS f s t RS fp + RS r + RS a 

1 1 r I I I 
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Scaled 


Scaled 


Total M-3 


Speed 


Speed 


Resistance 


(fps) 


(knots) 


at Full 






Scale 


VS. 


VS kts. 

1 


RS t 

i 



Residual PodFrictional StrutFnctional Reynolds# 
Component Component Component of the Strut 



Reynolds # 
of the Pod 



10.03 


13.52 


15.08 


16.2 


17.8 


18.57 


19.31 


20.13 


23.2 


26.87 


30.23 


33.61 


38.66 


42.01 


50.52 


58.79 


67.19 



5.939 


8.01 1 


8.932 


9.6 


10.543 


1 1.004 


11.441 


1 1.925 


13.743 


15.919 


17.91 


19.913 


22.905 


24.885 


29.927 


34.83 


39.803 



1805.61 


9515.06 


12890.08 


18956.26 


22563.1 1 


20988.62 


19383.41 


1821 1.4 


19668.69 


23779.04 


25097.62 


26706.78 


33361.84 


37792.85 


52161.65 


71750.23 


92642.78 



RS, 



572.22 


7356.37 


10243.43 


15926.77 


18951.13 


17075.09 


15172.89 


13660.76 


13727.86 


15948.15 


15322.77 


14773.05 


17826.12 


19629.28 


26437.1 1 


37490.84 


48556.41 



RS 



f P ; 



RS 



fs. 



637.02 


1109.41 


1357.69 


1552.15 


1847.58 


2000.34 


2150.67 


2322.7 


3024.61 . 


3976.3 


4952.65 


6034.61 


7835.39 


9146.31 


12907.73 


17139.31 


21997.35 



416.2 


721.46 


881.47 


1006.64 


1196.59 


1294.71 


1391.23 


1501.61 


1951.46 


2560.23 


3183.76 


3873.76 


5020.44 


5854.09 


8242.17 


10923.55 


13996.93 



RnS 



20396136 


27511067 


30673259 


32965848 


36207094 


37788190 


39290231 


40950382 


47195710 


54666388 


61504627 


68382394 


78659517 


85458229 


102771228 


119609899 


136685734 



RnS 



P; 



76475883 


10315351 


11501024 


12360637 


13575951 


14168787 


14731982 


15354460 


17696163 


20497315 


23061332 


25640170 


29493606 


32042802 


38534359 


44848066 


51250698 



WRJTE(m3fulfx) - RS ^ WRITE(m3fulv) - VS Writes the scaled data to a file for future use. 



RS ^ RS fp .RS fs 



Rs ft. 

FF. . = 

' RS t 



RS, 



RF. = 

1 RS, 



Percent of total resistance 
attributable to friction/wave. 



RS, 



CA =- 



RS, 



percent of total resistance that is attributable to correlation allowance 
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WRITE( m3 fillip) =FF. WRITE(m3fulrp) =RF. 



EHP • 



S hip. 



RS . VS. 

I- 1 



550 



VS. 



VS 



lets. 



EHP 



shiPj 



10.026 


13.523 


15.077 


16.204 


17.797 


18.574 


19.313 


20.129 


23.199 


26.87 1 


30.232 


33.613 


38.664 


42.006 


50.516 


58.793 


67.187 



5.939 


8.01 1 


8.932 


9.6 


10.543 


1 1 .004 


1 1.441 


1 1.925 


13.743 


15.919 


17.91 


19.913 


22.905 


24.885 


29.927 


34.83 


39.803 



32.913 


233.946 


353.357 


558.489 


730.1 13 


708.823 


680.632 


666.498 


829.613 


1 161.749 


1379.551 


1632.163 


2345.303 


2886 43 


4790.937 


7669.871 


1 1317.027 



WRITE(m3fulcap) =CA 



M-3 in full scale 
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Appendix H 



Approximation of Impingement Force 
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This Worksheet approximates the impingement force created by 
the rooster tail on the towing frame. 

p =64 — density of water 

ft 3 



V 



xr 




sec 



estimate of the maximum relative x component 
of velocity between the frame and water before impact. 



V 



xa 




sec 



water x component of velocity after the impact. 



t = 25 sec 



an approximation for the time of impact 
of the water on the frame 



a 



x 




acceleration of the water 



a x =40* ft* sec 



vol = 1 .25- in 6- in V ^ t an approximation for the volume of water that is 

impacted during the impact period. 



vol =0.1 3* ft 3 



m %v = p' vol 



m w = 8. 333-lb 



mass of the water that is impacted during the impact 
period. 



^impact m \v' a x 



Approximated Impact Force on the frame. 



^impact 10 36-lbt 
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Appendix I 



Numerical Codes 
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********** 



Q m ********** 



% This routine solves for the total resistance of an 8 bodied SLICE type 
% vessel. Solution to the friction resistance is based on the ITTC 1957 
% friction line. Solution to the wave resistance portion is based on the 
% solution to Mitchell's Integral proposed by Lunde in 1951 and expanded 
% by Wilson. 



% GENERAL INPUTS SECTION 

ls= input(’Input the slowest integer speed as x.xx (f/s) — ): 
hs=input(’The number of increments- ); 
lani=input('The incrementation (decimal)-); 

g=32.2; 



water=input('For Fresh Water enter 1. Salt Water enter 2: '); 
if water==l; 
rho= 1.9367; 
nu= 1.0804E-5. 
elseif water==2; 
rho= 1.9912: 
nu= 1.3343E-5. 



else 



disp(' 



') 



dispC ') 

dispC** **************************************** f ^ 
disp('A valid entry was not made for water type') 
disp('The program may not provide valid results’) 

disp(’******** **********************************^ 

end 



end 



dispC ’) 
dispC ') 
dispC ') 

disp('The body numbering is: l=fwd stbd strut') 
dispC 2=fwd port strut') 
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dispC 

disp(' 

disp(' 

dispC 

disp(' 

dispC 

dispC ') 

dispC ') 

dispC ') 

disp(’The 

disp(' 

disp(' 

dispC 

dispC 



3=aft stbd strut') 
4=aft port strut') 
5=f\vd stbd pod') 
6=fwd port pod') 
7=aft stbd pod’) 
8=aft port pod') 



global coordinate system is:') 

x positive in the direction of motion') 
y positive to port') 
z positive upward') 

The origin is at the undisturbed free surface') 



TS= inputCThe height of the struts (ft)-): 

LS= inputCThe length of the struts (ft)-): 

°bWSS=input('The wetted Surface of one strut (ft A 2)-); 

cs6=input('The 6th power coefficient for the strut potential fn='): 
cs5=input('The 5th power coefficient for the strut potential fn='); 
cs4=input('The 4th power coefficient for the strut potential fn-); 
cs3=input('The 3rd power coefficient for the strut potential fn='): 
cs2=input('The 2nd power coefficient for the strut potential fn- ); 
cslHnputCThe 1st power coefficient for the strut potential fn-); 
csO=input('The Oth power coefficient for the strut potential fn='): 

TP= inputCThe max diameter of the pods (ft)-); 

LP= inputCThe length of the pods (ft)-); 

°bWSP=input('The wetted surface of one pod (ft A 2)-); 

cp6=input('The 6th power coefficient for the pod potential fn-); 
cp5=input('The 5th power coefficient for the pod potential fn-); 
cp4=input(’The 4th power coefficient for the pod potential fn='); 
cp3=input('The 3rd power coefficient for the pod potential fn='); 
cp2=input('The 2nd power coefficient for the pod potential fn- ); 
cpl=input('The 1st power coefficient for the pod potential fn='); 
cpO=input('The 0th power coefficient for the pod potential fn='); 

sbl=input('The x setback of body 1 l 4 - or - ft)-); 
osl =input('The v offset of body 1 (+ or - ft)='); 
subl=input('The z submergence of bodvl (-or - ft)-); 
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sb2=input('The x setback of body 2 (- or - ft)='); 
os2=input(’The y offset of body 2 (- or - ft)-); 
sub2=input('The z submergence of body 2 (+ or - ft)-); 



sb3=input('The x setback of body 3 (- or - ft)- ); 
os3=input('The v offset of body 3 (- or - ft)=’); 
sub3=input('The z submergence of body3 (+ or - ft)-); 



sb4=input('The x setback of body 4 (- or - ft)- ); 
os4=input('The y offset of body 4 (- or - ft)-); 
sub4=input('The z submergence of body4 (+ or - ft)-); 



sb5=input(’The \ setback of body 5 (- or - ft)- ); 
os5=input('The y offset of body 5 I - or - ft)- ); 
sub5=input('The z submergence of bod> 5 (i- or - ft)='); 



sb6=input('The x setback of body 6 (- or - ft)—); 
os6=input(’The y offset of body 6 (- or - ft)- ); 
sub6=input(’The z submergence of body6 (+ or - ft)=’); 



sb7=input('The x setback of body 7 (- or - ft)=’); 
os7=input('The v offset of body 7 (- or - ft)-); 
sub7=input(’The z submergence of body 7 (+ or - ft)='); 



sb8=input('The x setback of body 8 (- or - ft)-); 
os8=input('The y offset of body 8 (* or - ft)- ); 
sub8=input('The z submergence of body8 ( + or - ft)-); 



tseg=input('The number of Theta Segments- ); 
zseg= inputCThe number of segments in Z- ); 
xseg= input('The number of segments in X- ); 
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% This section calls the routines that computes the Rw value 
% for each of the body and interface terms and the Rf value for 
% the body terms only. 



R\vb=0: 

Rwl 1=0; 

R\v22=0: 

Rw33=0: 

R\v44=0. 

Rw55=0: 

Rw66=0; 

Rw77=0: 

Rw88=0; 

Rvvs=0; 

Rwl 5=0; 

Rw26=0; 

Rw37=0; 

R\v48=0; 

R\vl=0; 

Rwl 3=0. 

Rwl 7=0; 

Rw53=0; 

Rw57=0; 

Rw24=0; 

Rw28=0; 

Rw64=0; 

Rw68=0; 

Rvvd=0; 

Rwl 4=0; 

Rwl 8=0; 

Rw54=0; 

Rw58=0; 

Rw23=0; 

Rw27=0; 

Rw63=0; 

Rw67=0; 

Rwt=0; 

Rwl 2=0; 

Rwl 6=0; 

Rw52=0; 

Rw56=0; 

Rw34=0; 

Rw'38=0; 

Rw74=0; 
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Rw78=0; 

Rf=0; 

Rfs=0; 

Rfp=0; 

Rt=0: 



%*********** BODY TERMS********* 

on eon e 

twotwo 

thrthr 

forfor 

fivfiv 

sixsix 

sevsev 

ateate 

%***** **SAME CORNER INTERACTION TERMS 

onefiv 

twosix 

thrsev 

forate 

%******LONGITUDINAL INTERACTION TERMS' 

onethr 

onesev 

fivthr 

fivsev 

twofor 

twoate 

sixfor 

sixate 



%*****DIAGONAL INTERACTION TERMS****** 

onefor 

oneate 

fivfor 

fivate 

twothr 

twosev 

sixthr 

sixsev 



%******TRANSVERSE DIAGONAL TERMS***** 



onetwo 

onesix 

fivtwo 

fivsix 

thrfor 

thrate 

sevfor 

sevate 



Rwb=Rwl 1 ~Rw22-R\v33-Rw44+R\v55+Rw66+Rw77+Rw88; 
R\vs=Rwl 5 -tR\v26-Rvv37-!-Rw48: 

Rwl=Rw 1 3-Rw 1 7 -Rvv53~Rw57+Rw24+R\v28+Rw64+Rw68: 
Rwd=Rw]4-Rwl 8-R\\54-Rw58+R\v23-*-Rw27+Rw63- | -Ryv67; 
Rwt=R\vl2-R\vl6-R\v52-R\v56+Rw34+R\v38+Rw74+Rw78; 



R\v=Rwb-Rws-!-Rvvl-R\vd-R\\l 



Rf=Rfs-Rfp 

Rt=R\v-!-Rf 



plot(V,Rt,'-\V,Rw,’- ’) 
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******************* pqdgeo* ************************* 



% 

% This file creates the offsets for the pod by obtaining nose 
% tail and parallel midbodv information. The nose is formed 
% using an ellipse and the tail is formed using a parabola. 

% This is then fit to a seventh order polynomial and the results compared. 

seg=input('The min number of segments (must be even) for the pod = ’); 
D=input('The max. diameter(ft) of the pod = '); 

L2D=input(’The desired Legnth to Diameter ratio of the pod = '); 
Lpmb=input(' The amount of Parallel Midbody(ft) in the pod = '); 
Pn=input(’ The decimal amount of the pod that is nose section = '); 

Lpod= D*L2D; 

Ln= Pn*(Lpod-Lpmb); 

Lt= Lpod-Ln-Lpmb: 



npn=input(' The nose shape exponenUbigger is fuller) = '); 
npt=input(' The tail shape exponenUbigger is fuller) = '); 

ppod=zeros( 1 ,8); 

\n=zeros( l.seg*2); 
offn=zeros( 1 ,seg*2); 
xpmb=zeros( 1 .seg/2); 
offpmb=zeros( 1 ,seg/2); 
xt=zeros( 1 ,seg*2); 
offt=zerost 1 ,seg*2); 
xpod=zeros( 1 .(seg*4-*-seg/2)); 
offpod=zeros( 1 ,(seg*4-i-seg/2)); 

Xpod=zeros( 1 ,seg*6); 

OFFpod=zeros( 1 ,seg*6); 

xn= (linspace(0,-Ln,seg*2))'; 
offn= D/2*( 1 -((Ln-abs(xn ))./Ln)/ npn). A ( 1 /npn); 



xpmb= (linspace(-Ln,-(Ln+Lpmb),seg/2))’; 
offpmb= D/2*xpmb./xpmb; 



xt = (linspace(-(Ln+Lpmb),-Lpod.seg*2))'; 
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offt= D/2 * ( 1 *((abs( \t)-Ln-Lpm b)./Lt ) A npt ); 



\pod=[\n:xpmb:xt]: 

offpod=[offn;offprnb:offt]; 



ppod=polyflt(xpod.offpod,7)' 



Xpod=( linspace(0,-Lpod.seg*6)); 

OFFpod=(ppod( l).*Xpod. A 7)+(ppod(2).*Xpod A 6)+(ppod(3).*Xpod. A 5)+... 
(ppod(4).*Xpod. ' 4)-*-(ppod(5).*Xpod. 3)+(ppod(6).*Xpod. A 2)+... 
(ppod(7).*Xpod. A 1 )+ppod(8). 



plot(Xpod.OFFpod/-'.Xpod.-OFFpod.':') 
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0/ 0 ***************** **STRUTGEO**** ********** ****** ****** 

% 

% This file creates the offsets for the strut by obtaining nose 
% tail and parallel midbodv information. The nose and tail are formed 
% using a parabola. This is then fit to a seventh order polynomial 
% and the results compared. 

segs=input(’The min number of segments (must be even) for the strut = ’); 
B=input('The max. beam! ft) of the strut = '); 

L2Ds=input(’The desired Legnth to Diameter ratio of the strut = '); 
Lpmbs=inputf The amount of Parallel Midbody(ft) in the strut = ’); 
Pns=inputf The decimal amount of the strut that is nose section = '); 

Lstrut= B*L2Ds: 

Lns= Pns*(Lstrut-Lpmbs); 

Lts= Lstrut-Lns-Lpmbs; 



nsn=input(' The nose shape exponent,(bigger is fuller) = '): 
nst=inputf The tail shape exponent,(bigger is fuller) = '); 



xns=zeros( 1 ,segs*2): 
offns=zeros( 1 ,segs*2); 
xpmbs=zeros( 1 ,segs/2); 
offpmbs=zeros( 1 .segs/2); 
xts=zeros( l,segs*2); 
offts=zeros(T ,segs*2); 
xstrut=zeros( 1 ,(segs*4+segs/2)); 
offstrut=zeros(l.(segs*4+segs/2)); 
Xstrut=zeros( 1 ,segs*6); 

OFFstrut=zeros( 1 .segs*6); 

xns= (linspace(0,-Lns,segs*2))'; 
offns= B/2*( 1 -((Lns-abs(xns))./Lns). A nsn); 



xpmbs= (linspace(-Lns,-(Lns+Lpmbs),segs/2))'; 
offpmbs= B/2*xpmbs./xpmbs; 



xts = (linspace(-(Lns+Lpmbs),-Lstrut,segs*2))‘; 
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offts= B/2*( l-((abs(xts)-Lns-Lpmbs)./Lts). A nst); 



xstrut=[xns;xpmbs:xts]; 

offstrut=[offns:offpmbs:offts]; 



ps=polyfit(xstrut,offstrut.7)’ 



Xstrut=(linspace(0.-Lstnit.segs*6))'; 

OFFstrut=(ps( 1 ).*Xstnit.''7)+(ps(2').*Xstrut. A 6)+(ps(3).*Xstnit. A 5)+... 
(ps(4).*Xstrut. A 4) J -(ps(5).*Xstrut. A 3)+(ps(6).*Xstrut. A 2) j -... 
(ps(7).*Xstrut)-(ps(8)); 



plot( Xstrut.OFFstn.it. '--'.xstrut.offstrut.'-'); 
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O,^****#*»***********0J S J£QJ V J£ m ********************* 



theta l=zeros(tseg+l ,1 ); 
zl=zeros(zseg+l,l); 
x 1 =zeros( xsegf 1,1); 
plx=zeros(xseg+l ,1); 
q 1 x=zeros( xseg+ 1,1); 
p 1 z=zeros( zseg+ 1,1); 
q 1 z=zeros( zseg+ 1,1); 
plt=zeros<tseg+l,l); 
q 1 t=zeros( tseg-*- 1,1); 
res 1 1 =zeros( tseg+ 1,1); 
q 1 At=zeros( tseg-^ 1,1); 
p 1 At=zeros( tseg-*- 1.1); 
p 1 Az=zeros( zseg-*- 1,1); 
q 1 Az=zeros( zseg+ 1.1); 
zl A=zeros( zseg+ 1,1); 
p 1 Ax=zeros( xseg+ 1,1); 
q 1 Ax=zeros( xseg+ 1,1); 
x 1 A=zeros( xseg+ 1,1); 
sara 1 =zeros( xseg-*- 1,1); 
katep 1 =zeros( xseg+ 1,1); 
kateq 1 =zeros( xseg-t- 1,1); 
greg 1 =zeros< xseg- 1,1); 
sara 1 A=zeros( xseg+ 1.1); 
katep 1 A=zeros( xseg 4 - 1.1); 
kateq 1 A=zeros( xseg-"- 1.1); 
greg 1 A=zeros( xseg- 4 1,1); 
V=zeros(hs-l-ls,l), 

Rw 1 1 =zeros( hs+ 1 -Is, 1 ); 



for iv=l :hs; 
V(iv)=ls-(iv*lani); 
k=g/(V(iv)'2); 
w=-V(iv) (2*pi); 



for it=l ;tseg+l ; 
theta 1 (it)=(pi/2)*((it-l )/tseg); 
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for izl= 1 izseg-H . 
zl(izl) = -TS*((izl-l )/zseg); 

for ixl=l:xseg+l; 
x 1 (ix 1 ) =-LS*((ixl-l)/xsegk 



saraKixl )=w*((cs6*xl (ix 1 ) A 6)+(cs5*xl(ixl )'5)+(cs4*xl (ix 1 ) A 4)+(cs3*x 1 ( ix 1 ) A 3)+(cs2 
* \ 1 ( ix 1 ) ' 2)+(cs 1 *x 1 (ix 1 ))+csO); 

katepltixl )=cos(k*(sbl *cos(thetal(it))-osl *sin(thetal(it)))*(sec(thetal(it))) A 2): 
kateqHixl )=sin(k*(sbl *cos(thetal(it))-osl *sin( theta 1 (it)))*(sec(thetal (it))) A 2); 
greg 1 (ix 1 )=exp(k*subl *( sec(thetal(it))) N 2); 

plx(ixl)=saral(ixl )*katepl(ixl )*gregl(ixl ): 
ql \(ixl >=sara 1 ( ix 1 )*kateq 1 < ix 1 )*gregl(ixl ); 

end °othis end closes the **** for ixl = l:xseg-l loop **** 

plz(izl)=sum(plx)/xseg; 
qlzfizl )=sum(qlx)/xseg; 

end °othis end closes the **** for iz 1 = 1 :zseg— 1 loop **** 

plt(it)=sum(plz)/zseg; 
q 1 1( it)=sum( q 1 z)/zseg; 

o 0 *******This starts the next body calculations inside the theta loop***** 



for izI A= 1 :zseg-l ; 
zlA(izlA) = -TS*((izl A-l )/zseg); 



for ixl A=1 :xseg+l ; 
xlA(ixlA) =-LS*((ixl A-l )/xseg); 

saral A(ixl A)=w*((cs6*xlA(ixl A) A 6)-*-(cs5*xl A(ixl A) A 5)-(cs4*xl A(ixl A) A 4)+(cs3*x 
1 A(ix 1 A) A 3)+(cs2*x 1 A( ix 1 A) A 2)+(csl *x 1 A(ix 1 A))-csO); 

katepl A( ix 1 A)=cos(k*(sbl *cos(thetal (it))+osl *sin(thetal(it)))*(sec(thetal(it))) A 2); 
kateql A(ixl A)=sin(k*(sbl *cos(thetal(it))+osl *sin(thetal(it)))*(sec(thetal(it))) A 2): 
gregl Af ix 1 A)=exp(k*subl * (sec( theta 1 (it))) A 2). 

p 1 Ax( ix 1 A)=sara 1 A( ix 1 A)*katep 1 At ix 1 A)* greg 1 A(ix 1 A); 
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q 1 Ax(ix 1 A)=sara 1 A(ix 1 A)*kateq 1 A(ix 1 A)*greg 1 A(ix 1 A); 

end %this end closes the **** for ixl A=1 :xseg+l loop **** 

p 1 Az(izl A)=sum(p 1 Ax)/'xseg; 
ql Az(izl A)=sum(ql Ax)/xseg; 

end %this end closes the ***♦ for izl A=1 izseg-^l loop **** 

p 1 At(it)=sum(pl Az)/zseg; 
ql At(it)=sum(ql Az)/zseg; 

resl l(it)=16*pi*rho*k A 2*((plt(it)*pl At(it))+(qlt(it)*ql At(it)))*(sec(thetal(it)) A 3); 

end %this end closes the **** for it= l:tseg-*-l loop **** 

Rwll(iv)= sum(resl 1 )/tseg 
end % this closes the **** for iv = ... loop ***** 
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<^*********************gJJ,£gJJ( m * 



Rw66=zeros(hs- 1 -Is. 1 ); 
res66=zeros(tseg- 1.1): 
theta 1 =zeros(tseg~ 1,1): 
z^zerosCzseg 4 -!.!): 
xl=zeros(xseg"-1.1); 
p 1 x=zeros(xseg- 1.1); 
q 1 x=zeros(xseg-l . 1 ): 
p 1 z=zeros(zseg^ 1.1); 
q 1 z=zeros(zseg~ 1.1); 
plt=zeros(tseg+l .1): 
qlt=zeros(tseg+l.l): 
q 1 At=zeros( tseg- 1.1): 
pi At=zeros(tseg- 1.1): 
p 1 Az=zeros( zseg- 1.1): 
q 1 Az=zeros(zseg- 1.1); 
z 1 A=zeros( zseg- 1 . 1 ): 
p 1 Ax=zeros(xseg- 1.1): 
q 1 Ax=zeros(xseg- 1.1): 
x 1 A=zeros(xseg- 1.1); 
sara 1 =zeros(xseg- 1.1); 
katep 1 =zeros(xseg- 1,1); 
kateq 1 =zeros(xseg- 1.1): 
greg 1 =zeros(xseg- 1,1); 
sara 1 A=zeros(xseg- 1,1); 
katep 1 A=zeros(xseg j - 1.1): 
kateq 1 A=zeros(xseg-l .1 ): 
greg 1 A=zeros( xseg- 1,1): 
V=zeros(hs+l-ls.l ): 



for iv=l:hs; 
V(iv)=ls+(iv*lani): 

k=g/(V(iv) A 2); 

w=-V(iv)/(2*pi); 



for it=l:tseg+l: 
theta 1 (it)=(pi/2)*((it-l )/tseg); 



for izl= 1 :zseg-*-l : 
zl(izl) = -TP*((izl-l)/zseg); 



for ixl=l :xseg+l : 
xl(ixl) =-LP*((ixl-l)/xseg); 



saral(ixl )=\v*((cp6*xl(ixl ) A 6)+(cp5*xl(ix 1 ) A 5)+(cp4*xl(ixl) A 4)+(cp3*xl(ixl) A 3)+(cp 
2*x 1 ( ix 1 ) A 2)+(cpl *x 1 (ix 1 ))-cpO); 

katepl(ixl )=cos(k*(sb6*cos(thetal(it))+os6*sin(thetal(it)))*(sec(thetal(it))) A 2); 
kateql(ixl)=sin(k*(sb6*cos( theta l(it))+os6*sin(thetal (it)))*(sec( theta l(it))) A 2); 
greg 1 (ix 1 )=exp(k*sub6* (sect theta 1 1 it))) A 2); 

plx(ixl )=saral(ixl )*katepl(ix 1 )*gregl(ixl ). 
qlxtixl )=saral(ixl )*kateql(ixl )*gregl (ix 1 ); 

end %this end closes the **** for ixl = l:xseg-l loop **** 

plztizl l=sum(plx)/xseg; 
q 1 z(iz 1 )=sum(ql x)/xseg: 

end %this end closes the **** for izl = l:zseg~l loop **** 

pi t(it)=sum(pl z)/zseg: 
q 1 t(it)=sum(q 1 z)/zseg; 

% *******This starts the next body calculations inside the theta loop***** 



for tzl A= 1 :zseg+l : 
zlA(izlA) = -TP*f(izl A-l)/zseg); 

for ixlA=l:xseg~l. 
xlA(ixlA) =-LP*((ixlA-l)/xseg); 

saral A(ixl A)=w*((cp6*xl A(ixl A) A 6)+(cp5*xl A(ixl A) A 5)+(cp4*xl A( ix 1 A) A 4)+(cp3* 
x 1 A(ix 1 A) A 3) 4 -(cp2*x 1 A(ix 1 A) A 2) J -(cpl *x 1 A(ix 1 A))+cpO); 

katepl A(ixl A)=cos(k*(sb6*cos(thetal(it))+os6*sin(thetal(it)))*(sec(thetal(it))) A 2); 
kateql A(ixl A)=sin(k*(sb6*cos(thetal(it))+os6*sin(thetal(it)))*(sec(thetal(it))) A 2); 
greg 1 A(ix 1 A)=exp(k*sub6*(sec(thetal (it))) A 2); 

pi Ax(ix 1 A)=saral A(ix 1 A)*katep 1 A(ix 1 A)* greg 1 A(ix 1 A); 
ql Ax(ixl A)=saral A(ixl A)*kateql A(ixl A)*gregl A(ixlA); 

end %this end closes the **** for ixl A=1 :xseg*l loop **** 
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p 1 Az( iz 1 A)=sum(p 1 Ax)/xseg; 
ql Az(izl A)=sum(ql Ax)/xseg; 

end °othis end closes the **** for izl A=1 :zseg+l loop **** 

p 1 At( it)=sum(p 1 Az)/zseg; 
q 1 At( it)=sum(q 1 Az)/zseg; 

res66(it)=16*pi*rho*k A 2*((plt(it)*pl At(it))-(qlt(it)*ql At(it)))*(sec(thetal(it)) A 3); 

end %this end closes the **** for it= l:tseg->-l loop **** 

Rvv66(iv)= sum(res66)/tseg 
end °o this closes the **** for iv = ... loop ***** 
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%**♦* ************** ***ONESIX m ************************** 



Rwl 6=zeros(hs+l -Is, 1 ); 
res 1 6=zeros( tseg+ 1,1); 
theta 1 =zeros(tseg+ 1,1); 
z 1 =zeros(zseg 4 - 1,1); 
xl=zeros(xseg+l,l); 
p 1 x=zeros(xseg- 4 - 1,1); 
q 1 x=zeros(xseg+ 1,1); 
p 1 z=zeros(zseg+ 1,1); 
q 1 z=zeros(zseg+l ,1 ); 
plt=zeros(tseg J -l,l); 
qlt=zeros(tseg^l,l); 
q 1 At=zeros(tseg- t - 1,1); 
p 1 At=zeros(tseg+ 1,1); 
p 1 Az=zeros( zseg+ 1,1); 
q 1 Az=zeros< zseg-*- 1,1); 
z 1 A=zeros( zseg + 1.1); 
pi Ax=zeros(xseg-*-l .1); 
q 1 Ax=zeros( xsegi- 1,1); 
x 1 A=zeros( xseg 4 - 1 , 1 ); 
sara 1 =zeros( x seg+ 1 , 1 ) ; 
katep 1 =zeros( xseg+ 1,1); 
kateql=zeros(xseg-rl,l ); 
greg 1 =zeros( xseg+ 1,1); 
sara 1 A=zeros(xseg+ 1,1); 
katep 1 A=zeros(xseg 4 - 1,1 ); 
kateq 1 A=zeros( xseg+ 1,1); 
greg 1 A=zeros( xseg-*- 1,1); 
V=zeros(hs-'-l-ls.l); 



for iv=l :hs; 
V(iv)=ls~(iv*lani); 



k=g/(V(iv) A 2); 

w=-V(iv)/(2*pi); 



for it=l :tseg+l ; 
theta 1 (it)=(pi/2)*((it- 1 )/tseg); 



for izl= l:zseg+l; 
zl(izl) = -TP*((izl-l)/zseg); 
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for ixl = l:xseg+l ; 
xl(ixl) =-LP*((ixl-l)/xseg); 



saral(ixl)=w*((cs6*xl(ixl) A 6)+(cs5*xl(ixl) A 5)+(cs4*xl(ixl) A 4)+(cs3*xl(ixl) A 3)+(cs2 
*x 1 (ix 1 ) A 2)+(csl *x 1 (ix 1 ))+csO); 

katep 1 (ix 1 )=cos(k*(sbl *cos(theta 1 (it))+osl *sin(thetal (it)))*(sec(thetal (it))) A 2); 
kateq l(ixl)=sin(k*(sbl*cos(thetal(it))+osl*sin(thetal(it)))*(sec(thetal(it))) A 2); 
greg 1 (ix 1 )=exp(k*sub 1 *(sec(theta 1 (it))) A 2); 

plx(ixl)=saral(ixl)*katepl(ixl )* gregl (ixl ); 
q 1 x(ix 1 )=saral (ix 1 )*kateq 1 (ix 1 )*greg 1 (ix 1 ); 

end %this end closes the **** for ixl=l:xseg+l loop **** 

plz(izl)=sum(plx)/xseg; 

qlz(izl)=sum(qlx)/xseg; 

end %this end closes the **** for izl = l:zseg+l loop **** 

p 1 t(it)=sum(p 1 z)/zseg; 
qlt(it)=sum(qlz)/zseg; 

% ****** *This starts the next body calculations inside the theta loop***** 



for izl A= 1 izsegH ; 
zlA(izlA) = -TP*((izl A-l)/zseg); 

for ixl A=1 :xseg+l ; 
xlA(ixlA) =-LP*((ixl A-l)/xseg); 

sara 1 A(ix 1 A)=w*((cp6*x 1 A(ix 1 A) A 6)+(cp5 *x 1 A(ix 1 A) A 5)-^(cp4*x 1 A(ix 1 A) A 4)+(cp3* 
x 1 A(ix 1 A) A 3)+(cp2*x 1 A(ix 1 A) A 2)- | -(cpl *x 1 A(ix 1 A))+cpO); 

katep 1 A(ixl A)=cos(k*(sb6*cos(thetal(it))+os6*sin(thetal(it)))*(sec(thetal(it))) A 2); 
kateq 1 A(ix 1 A)=sin(k*(sb6*cos( theta 1 (it))+os6*sin(thetal (it)))*(sec(theta 1 (it))) A 2); 
gregl A(ixl A)=exp(k*sub6*(sec(thetal(it))) A 2); 

pi Ax(ixl A)=saral A(ixl A)*katepl A(ixl A)* greg 1 A(ixl A); 
ql Ax(ixl A)=saral A(ixl A)* kateq 1 A(ixl A)* gregl A(ixl A); 

end %this end closes the **** for ixl A=1 :xseg+l loop **** 
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plAz(izl A)=sum(pl Ax)/xseg; 
qlAz(izl A)=sum(qlAx)/xseg; 



end %this end closes the **** for izl A=1 :zseg+l loop **** 

p 1 At( it )=sum (p 1 Az)/zseg: 
q 1 At( it)=sum(q 1 AzVzseg: 

resl6(it)=32*pi*rho*k A 2*((plt(it)*pl At(it))+(qlt(it)*ql At(it)))*(sec( theta l(it)) A 3); 

end %this end closes the **** for it= l:tseg-i-l loop **** 

Rwl6(iv)= sum(resl6)/tseg 
end % this closes the **** for iv = ... loop ***** 
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